Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

$\mathrm{Perturbo}$: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

We report Perturbo is a software package for first-principles calculations of charge transport and ultrafast carrier dynamics in materials. The current version focuses on electron–phonon interactions and can compute phonon-limited transport properties such as the conductivity, carrier mobility and Seebeck coefficient. It can also simulate the ultrafast nonequilibrium electron dynamics in the presence of electron–phonon scattering. Perturbo uses results from density functional theory and density functional perturbation theory calculations as input, and employs Wannier interpolation to reduce the computational cost. It supports norm-conserving and ultrasoft pseudopotentials, spin–orbit coupling, and polar electron–phonon interactions for bulk and 2D materials. Hybrid MPI plus OpenMP parallelization is implemented to enable efficient calculations on large systems (up to at least 50 atoms) using high-performance computing. Taken together, Perturbo provides efficient and broadly applicable ab initio tools to investigate electron–phonon interactions and carrier dynamics quantitatively in metals, semiconductors, insulators, and 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Two-step annealing induced structural rejuvenation: A cause for memory effect in metallic glasses

The structural origin of the memory effect is still elusive in glassy materials. Here, in this letter, the memory effect in a series of metallic glasses (MGs) is observed by two step annealing using differential scanning calorimetry (DSC) and molecular dynamics simulations. It demonstrates that the Kovacs peak can be reflected from not only the system potential energy but also the atomic dynamics, showing the memory effect could originate from the formation of more loose packing configurations around mobile atoms by structural rejuvenation in the second step of high temperature annealing. In contrast, the high-low temperature annealing protocol can only induce the monotonic relaxation. Therefore, our results pave a pathway to better understand the origin of memory effect in MGs and other glassy materials from their dynamics and local atomic structure.

Memory effect↗
DOE OSTI2021

How Accurate Are Approximate Density Functionals for Noncovalent Interaction of Very Large Molecular Systems?

Noncovalent intermolecular interactions are very important in many research areas. Therefore, it is vital to understand the extent to which approximate density functionals give a proper description of noncovalent interactions. Previous research has demonstrated that some approximate density functionals can predict usefully accurate interaction energies for many noncovalent systems; however, most of that work is limited to small and moderate-sized molecules. Very recently though, accurate benchmarks have become available for some very large molecules. Here, the present work applies 21 approximate density functionals to compute the binding energies of seven large molecular systems that have a number of atoms ranging from 200 to 910. The results are judged by comparison to the recently published CIM-DLPNO-CCSD(T) results, which are assumed to provide a reliable benchmark. The five most accurate methods among those tested are found to be PW6B95-D4, PW6B95-D3(BJ), revM11, M06-L, and MN15.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

Here, we present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modeling Interactions within and between Peptide Amphiphile Supramolecular Filaments

Many peptides are able to self-assemble into one-dimensional (1D) nanostructures, such as cylindrical fibers or ribbons of variable widths, but the relationship between the morphology of 1D objects and their molecular structure is not well understood. Here, we use coarse-grained molecular dynamics (CG-MD) simulations to study the nanostructures formed by self-assembly of different peptide amphiphiles (PAs). The results show that ribbons are hierarchical superstructures formed by laterally assembled cylindrical fibers. Simulations starting from bilayer structures demonstrate the formation of filaments, whereas other simulations starting from filaments indicate varying degrees of interaction among them depending on chemical structure. These interactions are verified by observations using atomic force microscopy of the various systems. The interfilament interactions are predicted to be strongest in supramolecular assemblies that display hydrophilic groups on their surfaces, while those with hydrophobic ones are predicted to interact more weakly as confirmed by viscosity measurements. The simulations also suggest that peptide amphiphiles with hydrophobic termini bend to reduce their interfacial energy with water, which may explain why these systems do not collapse into superstructures of bundled filaments. The simulations suggest that future experiments will need to address mechanistic questions about the self-assembly of these systems into hierarchical structures, namely, the preformation of interactive filaments vs equilibration of large assemblies into superstructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Strong Relativistic Effects in Lanthanide-Based Single-Molecule Magnets

Lanthanide-based single-molecule magnets (SMMs) are promising building blocks for quantum memory and spintronic devices. Designing lanthanide-based SMMs with long spin relaxation time requires a detailed understanding of their electronic structure, including the crucial role of the spin-orbit coupling (SOC). While traditional calculations of SOC using the perturbation theory applied to a solution of the non-relativistic Schrödinger equation are valid for light atoms, this approach is questionable for systems containing heavy elements such as lanthanides. We investigate the accuracy of the perturbation estimates of SOC by variationally solving the Dirac equation for the [DyO] + molecule, a prototype of a lanthanide-based SMM. We show that the energy splittings between the M J states involved in spin relaxation depend on the interplay between strong SOC and dynamic electron correlations. Here, we demonstrate that this interplay affects the resonances between the spin and vibrational transitions, and, therefore, the spin relaxation time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Mechanistic Advantages of Organotin Molecular EUV Photoresists

Here, extreme ultraviolet (EUV)-induced radiation exposure chemistry in organotin-oxo systems, represented by the archetypal [(R-Sn) 12 O 14 (OH) 6 ](A) 2 cage, has been investigated with density functional theory. Upholding existing experimental evidence of Sn-C cleavage-dominant chemistry, computations have revealed that either electron attachment or ionization can single-handedly trigger tin-carbon bond cleavage, partially explaining the current EUV sensitivity advantage of metal oxide systems. We have revealed that tin atoms at different parts of the molecule react differently to ionization and electron attachment and have identified such selectivity as a result of local coordination chemistry instead of the macro geometry of the molecule. An ionization-deprotonation pathway has also been identified to explain the observed evolution of an anion conjugate acid upon exposure and anion mass dependence in resist sensitivity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Origin of giant electric-field-induced strain in faulted alkali niobate films

A large electromechanical response in ferroelectrics is highly desirable for developing high-performance sensors and actuators. Enhanced electromechanical coupling in ferroelectrics is usually obtained at morphotropic phase boundaries requiring stoichiometric control of complex compositions. Recently it was shown that giant piezoelectricity can be obtained in films with nanopillar structures. Here, we elucidate its origin in terms of atomic structure and demonstrate a different system with a greatly enhanced response. This is in non-stoichiometric potassium sodium niobate epitaxial thin films with a high density of self-assembled planar faults. A giant piezoelectric coefficient of ~1900 picometer per volt is demonstrated at 1 kHz, which is almost double the highest ever reported effective piezoelectric response in any existing thin films. The large oxygen octahedral distortions and the coupling between the structural distortion and polarization orientation mediated by charge redistribution at the planar faults enable the giant electric-field-induced strain. Our findings demonstrate an important mechanism for realizing the unprecedentedly giant electromechanical coupling and can be extended to many other material functions by engineering lattice faults in non-stoichiometric compositions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A graph neural network-state predictive information bottleneck (GNN-SPIB) approach for learning molecular thermodynamics and kinetics

Molecular dynamics simulations offer detailed insights into atomic motions but face timescale limitations. Enhanced sampling methods have addressed these challenges but even with machine learning, they often rely on pre-selected expert-based features. Here, in this work, we present a Graph Neural Network-State Predictive Information Bottleneck (GNN-SPIB) framework, which combines graph neural networks and the state predictive information bottleneck to automatically learn low-dimensional representations directly from atomic coordinates. Tested on three benchmark systems, our approach predicts essential structural, thermodynamic and kinetic information for slow processes, demonstrating robustness across diverse systems. The method shows promise for complex systems, enabling effective enhanced sampling without requiring pre-defined reaction coordinates or input features.

Zou, Ziyue↗
DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Ultrafast state-selective tunneling in two-dimensional semiconductors with a phase- and amplitude-controlled THz-scanning tunneling microscope

THz-pulse driven scanning tunneling microscopy (THz-STM) enables access to the ultrafast quantum dynamics of low-dimensional material systems at simultaneous ultrafast temporal and atomic spatial resolution. State-selective tunneling requires precise amplitude and phase control of the THz pulses combined with quantitative near-field waveform characterization. Here, we employ our state-of-the-art THz-STM with multi-MHz repetition rates, efficient THz generation, and precisely tunable THz waveforms to investigate a single sulfur vacancy in monolayer MoS2. We demonstrate that 2D transition metal dichalcogenides (TMDs) are an ideal platform for near-field waveform sampling by THz cross-correlation. Furthermore, we determine the THz voltage via QEV scans, which measure the THz rectified charge Q as a function of THz field amplitude E and dc bias Vdc. Mapping the complex energy landscape of localized states with a resolution down to 0.01 electrons per pulse facilitates state-selective tunneling to the HOMO and LUMO orbitals of a charged sulfur vacancy.

Materials Science↗
DOE OSTI2025

Keldysh tuning of photoluminescence in a lead halide perovskite crystal

In 1964, Keldysh laid the groundwork for strong-field physics in atomic, molecular, and solid-state systems by delineating a ubiquitous transition from multiphoton absorption to quantum electron tunneling under intense AC driving forces. While both processes in semiconductors can generate carriers and result in photon emission through electron-hole recombination, the low quantum yields in most materials have hindered direct observation of the Keldysh crossover. Leveraging the large quantum yields of photoluminescence in lead halide perovskites, we show that we can not only induce bright light emission from extreme sub-bandgap light excitation but also distinguish between photon-induced and electric-field-induced processes. Our results are rationalized by the Landau-Dykhne formalism, providing insights into the non-equilibrium dynamics of strong-field light-matter interactions. These findings open new avenues for light upconversion and sub-bandgap photon detection, highlighting the potential of lead halide perovskites in advanced optoelectronic applications.

FOS: Physical sciences↗
DOE OSTI2025

Review: moiré-of-moiré superlattice in twisted trilayer graphene

When a third layer of graphene is transferred on top of twisted bilayer graphene with a second twist, the atomic and electronic structure of the system can be significantly enriched. Generally, the two coexisting moiré superlattices give rise to a higher-order moiré-of-moiré (MoM) superlattice, with a plethora of new length scales and associated novel quantum phenomena. This article reviews the current theoretical understanding and experimental observations of twisted trilayer graphene MoM (or supermoiré) superlattices, and the theoretical predictions and ongoing experimental efforts for unraveling rich exotic quantum phenomena and underlying physics.

electron correlation↗
DOE OSTI2026

Probing Excited-State Dynamics of Transmon Ionization

The fidelity and quantum nondemolition character of the dispersive readout in circuit QED are limited by unwanted transitions to highly excited states at specific photon numbers in the readout resonator. This observation can be explained by multiphoton resonances between computational states and highly excited states in strongly driven nonlinear systems, analogous to multiphoton ionization in atoms and molecules. In this work, we utilize the multilevel nature of high-𝐸 𝐽 /𝐸 𝐶 transmons to probe the excited-state dynamics induced by strong drives during readout. With up to ten resolvable states, we quantify the critical photon number of ionization, the resulting state after ionization, and the fraction of the population transferred to highly excited states. Moreover, using pulse shaping to control the photon number in the readout resonator in the high-power regime, we tune the adiabaticity of the transition and verify that transmon ionization is a Landau-Zener-type transition. We further extend these methods to a typical transmon with 𝐸 𝐽 /𝐸 𝐶 ≈ 55 and probe the offset-charge dependence of ionization dynamics in a timed-resolved manner. Our experimental results agree well with the theoretical prediction from a semiclassical driven transmon model and may guide future exploration of strongly driven nonlinear oscillators.

cavity quantum electrodynamics↗
DOE OSTI2020

Combining the synergistic control capabilities of modeling and experiments: Illustration of finding a minimum-time quantum objective

A common way to manipulate a quantum system, for example spins or artificial atoms, is to use properly tailored control pulses. In order to accomplish quantum information tasks before coherence is lost, it is crucial to implement the control in the shortest possible time. Here we report the near time-optimal preparation of a Bell state with fidelity higher than 99% in an NMR experiment, which is feasible by combining the synergistic capabilities of modeling and experiments operating in tandem. In this work, the pulses preparing the Bell state are found by experiments that are recursively assisted with a gradient-based optimization algorithm working with a model. Thus, we exploit the interplay between model-based numerical optimal design and experimental-based learning control. Utilizing the balanced synergism between the dual approaches, as dictated by the case specific capabilities of each approach, should have broad applications for accelerating the search for optimal quantum controls.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗