Quantum computation of stopping power for inertial fusion target design.

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it-one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [PRX Quant. 2, 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born-Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

Nonequilibrium Dynamics of Electron Emission from Cold and Hot Graphene under Proton Irradiation.

Characteristic properties of secondary electrons emitted from irradiated two-dimensional materials arise from multi-length and multi-time-scale relaxation processes that connect the initial nonequilibrium excited electron distribution with their eventual emission. To understand these processes, which are critical for using secondary electrons as high-resolution thermalization probes, we combine first-principles real-time electron dynamics with irradiation experiments. Our data for cold and hot proton-irradiated graphene show signatures of kinetic and potential emission and generally good agreement for electron yields between experiment and theory. The duration of the emission pulse is about 1.5 fs, which indicates high time resolution when used as a probe. Our newly developed method to predict kinetic energy spectra shows good agreement with electron and ion irradiation experiments and prior models. We find that the lattice temperature significantly increases secondary electron emission, whereas electron temperature has a negligible effect.

Identifying native point defect configurations inα-alumina.

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case ofα-Al2O3by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. These new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

Real-Time Exciton Dynamics with Time-Dependent Density-Functional Theory.

Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as -1/k^{2} to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate, for various materials, that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.

Anomalous Stopping and Charge Transfer in Proton-Irradiated Graphene.

We use first-principles calculations to uncover and explain a new type of anomalous low-velocity stopping effect in proton-irradiated graphene. We attribute a shoulder feature that occurs exclusively for channeling protons to enhanced electron capture from σ- and π-orbitals. Our analysis of electron emission indicates that backward emission is more sensitive to proton trajectory than forward emission and could thus produce higher contrast images in ion microscopy. For slow protons, we observe a steep drop in emission, consistent with predictions from analytical models.