SEARCH · Engineering Papers
Results for “Dynamics”
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.
Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
Dynamics of Phospholipid Membranes beyond Thermal Undulations
Here, we investigated the molecular dynamics of unilamellar liposomes by neutron spin echo spectroscopy. We report the first experimental evidence of a short-range motion at the length scale of the size of the headgroup of a lipid. The associated mean squared displacement shows a t 0.26 dependence in the pico- to nanosecond region that indicates another process beyond the predictions of the Zilman–Granek (ZG) model (t 0.66 ) and translational diffusion (t 1 ). A comparison with theory shows that the observed low exponent is associated with a non-Gaussian transient trapping of lipid molecules in a local area and supports the continuous time random walk model. The analysis of the mean squared displacement leads to the important conclusion that the friction at the interface between water and liposomes plays a minor role. Center of mass diffusion of liposomes and transient trapping of lipids define the range in which the ZG model can be applied to analyze membrane fluctuations.
Dynamic Control of Optical Response in Layered Metal Chalcogenide Nanoplates
Explore the source record for details and available documents.
Endogenous Dynamic Nuclear Polarization for Natural Abundance 17O and Lithium NMR in the Bulk of Inorganic Solids
Explore the source record for details and available documents.
Photoinduced Dynamics of Charge Separation: From Photosynthesis to Polymer–Fullerene Bulk Heterojunctions
Understanding charge separation and charge transport is crucial for improving the efficiency of organic solar cells. Their active media are based on organic molecules and polymers, serving as both light-absorbing and transport layers. The charge-transfer (CT) states play an important role, being intermediate for free carrier generation and charge recombination. Here, we use light-induced electron paramagnetic resonance spectroscopy to study the CT dynamics in blends of the polymers P3HT, PCDTBT, and PTB7 with the fullerene derivative C-60-PCBM. Time-resolved EPR measurements show strong spin-polarization patterns for all polymer-fullerene blends, confirming predominant generation of singlet CT states and partial orientation ordering near the donor-acceptor interface. These observations allow a comparison with charge separation processes in molecular donor-acceptor systems and in natural and artificial photosynthetic assemblies, and thus the elucidation of the initial steps of sequential CT in organic photovoltaic materials.
Molecular dynamics simulations of a hydrophilic MIL-160-based membrane demonstrate pressure-dependent selective uptake of industrially relevant greenhouse
Oral presentation for ACS Spring Meeting
Structure, solvation, and dynamics of Mg{sup 2+}, Ca{sup 2+}, Sr{sup 2+}, and Ba{sup 2+} complexes with 3-hydroxyflavone and perchlorate anion in acetonitrile medium: A molecular dynamics simulation study
Molecular dynamics simulations of complexes of Mg{sup 2+}, Ca{sup 2+}, Sr{sup 2+}, and Ba{sup 2+} with 3-hydroxyflavone (flavonol, 3HF) and ClO {sub 4}{sup −} in acetonitrile were performed. The united atoms force field model was proposed for the 3HF molecule using the results of DFT quantum chemical calculations. 3HF was interpreted as a rigid molecule with two internal degrees of freedom, i.e., rotation of the phenyl ring and of the OH group with respect to the chromone moiety. The interatomic radial distribution functions showed that interaction of the cations with flavonol occurs via the carbonyl group of 3HF and it is accompanied with substitution of one of the acetonitrile molecules in the cations’ first solvation shells. Formation of the cation–3HF complexes does not have significant impact on the rotation of the phenyl ring with respect to the chromone moiety. However, the orientation of the flavonol's OH-group is more sensitive to the interaction with doubly charged cations. When complex with Mg{sup 2+} is formed, the OH-group turns out of the plane of the chromone moiety that leads to rupture of intramolecular H-bond in the ligand molecule. Complexation of Ca{sup 2+}, Sr{sup 2+}, and BaClO {sub 4}{sup +} with 3HF produces two structures with different OH-positions, as in the free flavonol with the intramolecular H-bond and as in the complex with Mg{sup 2+} with disrupted H-bonding. It was shown that additional stabilization of the [MgClO{sub 4}(3HF)]{sup +} and [BaClO{sub 4}(3HF)]{sup +} complexes is determined by strong affinity of perchlorate anion to interact with flavonol via intracomplex hydrogen bond between an oxygen atom of the anion and the hydrogen atom of the 3-hydroxyl group. Noticeable difference in the values of the self-diffusion coefficients for Kt{sup 2+} from one side and ClO {sub 4}{sup −}, 3HF, and AN in the cations’ coordination shell from another side implies quite weak interaction between cation, anion, and ligands in the investigated complexes.
Vibrational dynamics of zero-field-splitting hamiltonian in gadolinium-based MRI contrast agents from ab initio molecular dynamics
The electronic relaxation of gadolinium complexes used as MRI contrast agents was studied theoretically by following the short time evolution of zero-field-splitting parameters. The statistical analysis of ab initio molecular dynamics trajectories provided a clear separation between static and transient contributions to the zero-field-splitting. For the latter, the correlation time was estimated at approximately 0.1 ps. The influence of the ligand was also probed by replacing one pendant arm of our reference macrocyclic complex by a bulkier phosphonate arm. In contrast to the transient contribution, the static zero-field-splitting was significantly influenced by this substitution.
Dissociative chemisorption of methane on metal surfaces: Tests of dynamical assumptions using quantum models and ab initio molecular dynamics
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Nonlinear space charge dynamics in mixed ionic-electronic conductors: Resistive switching and ferroelectric-like hysteresis of electromechanical response
We performed self-consistent modelling of nonlinear electrotransport and electromechanical response of thin films of mixed ionic-electronic conductors (MIEC) allowing for steric effects of mobile charged defects (ions, protons, or vacancies), electron degeneration, and Vegard stresses. We establish correlations between the features of the nonlinear space-charge dynamics, current-voltage, and bending-voltage curves for different types of the film electrodes. A pronounced ferroelectric-like hysteresis of the bending-voltage loops and current maxima on the double hysteresis current-voltage loops appear for the electron-transport electrodes. The double hysteresis loop with pronounced humps indicates a memristor-type resistive switching. The switching occurs due to the strong nonlinear coupling between the electronic and ionic subsystems. A sharp meta-stable maximum of the electron density appears near one open electrode and moves to another one during the periodic change of applied voltage. Our results can explain the nonlinear nature and correlation of electrical and mechanical memory effects in thin MIEC films. The analytical expression proving that the electrically induced bending of MIEC films can be detected by interferometric methods is derived.
Mean ionic activity coefficients in aqueous NaCl solutions from molecular dynamics simulations
The mean ionic activity coefficients of aqueous NaCl solutions of varying concentrations at 298.15 K and 1 bar have been obtained from molecular dynamics simulations by gradually turning on the interactions of an ion pair inserted into the solution. Several common non-polarizable water and ion models have been used in the simulations. Gibbs-Duhem equation calculations of the thermodynamic activity of water are used to confirm the thermodynamic consistency of the mean ionic activity coefficients. While the majority of model combinations predict the correct trends in mean ionic activity coefficients, they overestimate their values at high salt concentrations. The solubility predictions also suffer from inaccuracies, with all models underpredicting the experimental values, some by large factors. These results point to the need for further ion and water model development.
Development of a model electronic Hamiltonian for understanding electronic relaxation dynamics of [Fe(bpy){sub 3}]{sup 2+} through molecular dynamics simulations
A model electronic Hamiltonian of [Fe(bpy){sub 3}]{sup 2+}, which was recently refined for use in molecular dynamics simulations, is reviewed with some additional results. In particular, the quality of the refined model Hamiltonian is examined in terms of the vibrational frequencies and solvation structures of the lowest singlet and quintet states.
Understanding misfit strain releasing mechanisms via molecular dynamics simulations of CdTe growth on {112}zinc-blende CdS
Molecular dynamics simulations have been used to analyse microstructures of CdTe films grown on {112} surfaces of zinc-blende CdS. Interestingly, CdTe films grow in <331> orientations as opposed to <112> epitaxial orientations. At the CdTe-{331}/CdS-{112} interface, however, there exists an axis that is parallel to the <110> orientation of both CdS and CdTe. It is the direction orthogonal to this <110> that becomes different, being <116> for CdTe and <111> for CdS, respectively. Missing CdTe-{110} planes are found along the <110> axis, suggesting that the misfit strain is released by the conventional misfit dislocation mechanism along this axis. In the orthogonal axis, the misfit strain is found to be more effectively released by the new grain orientation mechanism. Our finding is supported by literature experimental observations of the change of growth direction when Cd 0.96 Zn 0.04 Te films are deposited on GaAs. Lastly the analyses of energetics clearly demonstrate the cause for the formation of the new orientation, and the insights gained from our studies can help understand the grain structures experimentally observed in lattice mismatched systems.
System-size corrections for self-diffusion coefficients calculated from molecular dynamics simulations: The case of CO{sub 2}, n-alkanes, and poly(ethylene glycol) dimethyl ethers
Molecular dynamics simulations were carried out to study the self-diffusion coefficients of CO{sub 2}, methane, propane, n-hexane, n-hexadecane, and various poly(ethylene glycol) dimethyl ethers (glymes in short, CH{sub 3}O–(CH{sub 2}CH{sub 2}O){sub n}–CH{sub 3} with n = 1, 2, 3, and 4, labeled as G1, G2, G3, and G4, respectively) at different conditions. Various system sizes were examined. The widely used Yeh and Hummer [J. Phys. Chem. B 108, 15873 (2004)] correction for the prediction of diffusion coefficient at the thermodynamic limit was applied and shown to be accurate in all cases compared to extrapolated values at infinite system size. The magnitude of correction, in all cases examined, is significant, with the smallest systems examined giving for some cases a self-diffusion coefficient approximately 15% lower than the infinite system-size extrapolated value. The results suggest that finite size corrections to computed self-diffusivities must be used in order to obtain accurate results.
Dynamic mapping of conical intersection seams: A general method for incorporating the geometric phase in adiabatic dynamics in polyatomic systems
The incorporation of the geometric phase in single-state adiabatic dynamics near a conical intersection (CI) seam has so far been restricted to molecular systems with high symmetry or simple model Hamiltonians. This is due to the fact that the ab initio determined derivative coupling (DC) in a multi-dimensional space is not curl-free, thus making its line integral path dependent. In a recent work [C. L. Malbon et al. , J. Chem. Phys. 145, 234111 (2016)], we proposed a new and general approach based on an ab initio determined diabatic representation consisting of only two electronic states, in which the DC is completely removable, so that its line integral is path independent in the simply connected domains that exclude the CI seam. Then with the CIs included, the line integral of the single-valued DC can be used to construct the complex geometry-dependent phase needed to exactly eliminate the double-valued character of the real-valued adiabatic electronic wavefunction. This geometry-dependent phase gives rise to a vector potential which, when included in the adiabatic representation, rigorously accounts for the geometric phase in a system with an arbitrary locus of the CI seam and an arbitrary number of internal coordinates. Here, we demonstrate this approach in a three-dimensional treatment of the tunneling facilitated dissociation of the S 1 state of phenol, which is affected by a C s symmetry allowed but otherwise accidental seam of CI. Here, since the space is three-dimensional rather than two-dimensional, the seam is a curve rather than a point. The nodal structure of the ground state vibronic wavefunction is shown to map out the seam of CI.
The role of sub-surface hydrogen on CO 2 reduction and dynamics on Ni(110): An ab initio molecular dynamics study
The catalytic reduction of carbon dioxide is a crucial step in many chemical industrial reactions, such as methanol synthesis, the reverse water-gas shift reaction (rWGS) and formic acid synthesis. Here, we investigate the role of bulk hydrogen, where hydrogen atoms are found deep inside a metal surface as opposed to subsurface ones, upon CO 2 reduction over a Ni(110) surface using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. While it has previously been shown that subsurface hydrogen stabilizes CO 2 and can aid in overcoming reaction barriers, the role of bulk hydrogen is less studied and thus unknown with regard to CO 2 reduction. Furthermore, we find that the presence of bulk hydrogen can significantly alter the electronic structure of the Ni(110) surface, particularly the work function and d-band center, such that CO 2 adsorbs more strongly to the surface and is more easily reduced. Our results show an enhanced CO 2 dissociation in the presence of bulk hydrogen, shedding light on a hitherto underappreciated mechanistic pathway for CO 2 reduction on metal surfaces.