A procedure for defining virtual spaces, and the periodic one-electron and two-electron integrals, for plane-wave second quantized Hamiltonians has been developed, and it was validated using full configuration interaction (FCI) calculations, as well as executions of variational quantum eigensolver (VQE) circuits on Quantinuum’s ion trap quantum computers accessed through Microsoft’s Azure Quantum service. This work is an extension to periodic systems of a new class of algorithms in which the virtual spaces were generated by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. In this extension, the integration of the first Brillouin zone is automatically incorporated into the two-electron integrals. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, that were able to capture a significant amount of correlation. The focus in this manuscript is on comparing the simulations of small molecules calculated with plane-wave basis sets with large periodic unit cells at the Gamma-point, including images, to results for plane-wave basis sets with aperiodic unit cells. The results for this approach were promising, as we were able to obtain good agreement between periodic and aperiodic results for an LiH molecule. Calculations performed on the Quantinuum H1-1 quantum computer produced surprisingly good energies, in which the error mitigation played a small role in the quantum hardware calculations and the (noisy) quantum simulator results. Using a modest number of circuit runs (500 shots), we reproduced the FCI values for the 1 COVO Hamiltonian with an error of 11 milliHartree, which is expected to improve with a larger number of circuit runs.

Plasticity modelling has long relied on phenomenological models based on ad-hoc assumption of constitutive relations, which are then fitted to limited data. Other work is based on the consideration of physical mechanisms which seek to establish a physical foundation of the observed plastic deformation behavior through identification of isolated defect processes (’mechanisms’) which are observed either experimentally or in simulations and then serve to formulate so-called physically based models. Neither of these approaches is adequate to capture the complexity of plastic deformation which belongs into the realm of emergent collective phenomena, and to understand the complex interplay of multiple deformation pathways which is at the core of modern high performance structural materials. Data based approaches offer alternative pathways towards plasticity modelling whose strengths and limitations we explore here for a simple example, namely the interplay between rate and dislocation density dependent strengthening mechanisms in fcc metals.

The opportunities afforded by near-term quantum computers to calculate the ground-state properties of small molecules depend on the structure of the computational ansatz as well as the errors induced by device noise. Here we investigate the behavior of these noisy quantum circuits using numerical simulations to estimate the accuracy and fidelity of the prepared quantum states relative to the ground truth obtained by conventional means. We implement several different types of ansatz circuits derived from unitary coupled cluster theory for the purposes of estimating the ground-state energy of sodium hydride using the variational quantum eigensolver algorithm. We show how relative error in the energy and the fidelity scale with the levels of gate-based noise, the internuclear configuration, the ansatz circuit depth, and the parameter optimization methods.

The recently introduced coupled cluster (CC) downfolding techniques for reducing the dimensionality of quantum many-body problems recast the CC formalism in the form of the renormalization procedure allowing, for the construction of effective (or downfolded) Hamiltonians in small-dimensionality sub-space, usually identified with the so-called active space, of the entire Hilbert space. The resulting downfolded Hamiltonians integrate out the external (out-of-active-space) Fermionic degrees of freedom from the internal (in-the-active-space) parameters of the wave function, which can be determined as components of the eigenvectors of the downfolded Hamiltonians in the active space. This paper will discuss the extension of non-Hermitian (associated with standard CC formulations) and Hermitian (associated with the unitary CC approaches) downfolding formulations to composite quantum systems commonly encountered in materials science and chemistry. The non-Hermitian formulation can provide a platform for developing local CC approaches, while the Hermitian one can serve as an ideal foundation for developing various quantum computing applications based on the limited quantum resources. We also discuss the algorithm for extracting the semi-analytical form of the inter-electron interactions in the active spaces.

In the Portevin–Le Chatelier (PLC) effect sample plastic deformation takes place via localized bands. We present a model to account for band dynamics and the variability the bands exhibit. The approach is tuned to account for strain hardening and the strain-rate dependence for the case of so-called type A (propagating) bands. The main experimental features of the fluctuations are a reduction with strain and increase with the strain rate which is reproduced by a model of plastic deformation with Dynamic Strain Aging, including disorder as a key parameter. Extensions are discussed as are the short-comings in reproducing detailed avalanche statistics.

Ultra pure metals have various applications, e. g. as electrical conductors. Crystallization from the melt, e. g. via zone melting, using the segregation of impurities at the solidification front is the basic mechanism behind different technical processes for the refining of metals and semi-metals. In this paper, we focus on a crystallization methodology with a gas cooled tube (“cooled finger”) dipped into a metallic melt in a rotating crucible. The necessary requirement for purification in a solidification process is a morphologically stable solidification front. This is the only way to enable macroscopic separation of the impurities, e. g. by convection. For cellular or dendritic solidification morphologies, the segregated impurities are trapped into the interdendritic melt and remain as microsegregations in the solidified metal. Morphological stability depends on the temperature gradient G at the solidification front, the solidification front velocity V front and thermodynamic alloy properties like the segregation coefficients of the impurity elements. To quantify the impact of these parameters on the morphological evolution, especially on the planar/cellular transition and thus on microsegregation profiles, phase field simulations coupled to a thermodynamic database are performed for an aluminium melt with three impurities, Si, Mn and Fe. In particular, we have investigated the morphology evolution from the start of solidification at the cooled finger towards a stationary growth regime, because in the technical process a significant fraction of the melt solidifies along the initial transient. To solve the transient long range temperature evolution on an experimental length scale, the temperature field has been calculated using the homoenthalpic approach together with a 1D temperature field approximation. The simulations provide the process window for an energy efficient purification process, i. e. low thermal gradients, and elucidate the benefit of melt convection.

Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N-spin system, the constant-depth circuit contains only mathcal {O}(N^{2}) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Microstructure simulations for quaternary alloys are still a challenge, although it is of high importance for alloy development. This work presents a Phase field (PF) approach capable of resolving phase transformation in a multicomponent system with a simple and effective way to include the thermodynamic and kinetic information for such a complex system. The microstructure evolution during diffusional transformation between FCC and BCC phase at 700 °C for AlCrFeNi alloys was simulated, accounting for composition dependence and off-diagonal terms in the diffusion tensor. The reliability of the presented PF method is validated by comparing the 1-D simulation results with simulations by Diffusion Module (DICTRA) of Thermo-Calc Software. Additionally, 2-D PF simulations of precipitate growth and Ostwald ripening are performed for different alloy systems, and the coarsening behavior is compared. Results showed that thermodynamic and kinetic information is accurately described in the applied PF method. The simulation results show that the diffusion behavior is influenced evidently by variations in the amounts of the different elements in the system. These findings demonstrate the necessity of applying accurate thermodynamic and kinetic models to fully understand the complex interdiffusion behavior in high and medium entropy alloys.