Search NASA⌕ Search

SEARCH · Search NASA

Results for “classical”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Nonadiabatic nuclear dynamics of the ammonia cation studied by surface hopping classical trajectory calculations

The Landau–Zener (LZ) type classical-trajectory surface-hopping algorithm is applied to the nonadiabatic nuclear dynamics of the ammonia cation after photoionization of the ground-state neutral molecule to the excited states of the cation. The algorithm employs a recently proposed formula for nonadiabatic LZ transition probabilities derived from the adiabatic potential energy surfaces. The evolution of the populations of the ground state and the two lowest excited adiabatic states is calculated up to 200 fs. The results agree well with quantum simulations available for the first 100 fs based on the same potential energy surfaces. Three different time scales are detected for the nuclear dynamics: Ultrafast Jahn–Teller dynamics between the excited states on a 5 fs time scale; fast transitions between the excited state and the ground state within a time scale of 20 fs; and relatively slow partial conversion of a first-excited-state population to the ground state within a time scale of 100 fs. Beyond 100 fs, the adiabatic electronic populations are nearly constant due to a dynamic equilibrium between the three states. The ultrafast nonradiative decay of the excited-state populations provides a qualitative explanation of the experimental evidence that the ammonia cation is nonfluorescent.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Si nanowire growth on sapphire: Classical incubation, reverse reaction, and steady state supersaturation

Si nanowire growth on sapphire substrates by the vapor-liquid-solid (VLS) method using Au catalyst particles has been studied. Sapphire was chosen as the substrate to ensure that the vapor phase is the only source of Si. Three hitherto unreported observations are described. First, an incubation period of 120–480 s, which is shown to be the incubation period as defined in classical nucleation theory, is reported. This incubation period permits the determination of a desolvation energy of Si from Au-Si alloys of 15 kT. Two, transmission electron microscopy studies of incubation, point to Si loss by reverse reaction as an important part of the mechanism of Si nanowire growth by VLS. Three, calculations using these physico-chemical parameters determined from incubation and measured steady state growth rates of Si nanowires show that wire growth happens from a supersaturated catalyst droplet.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

The structure and IR signatures of the arginine-glutamate salt bridge. Insights from the classical MD simulations

Salt bridges and ionic interactions play an important role in protein stability, protein-protein interactions, and protein folding. Here, we provide the classical MD simulations of the structure and IR signatures of the arginine (Arg)–glutamate (Glu) salt bridge. The Arg-Glu model is based on the infinite polyalanine antiparallel two-stranded β-sheet structure. The 1 μs NPT simulations show that it preferably exists as a salt bridge (a contact ion pair). Bidentate (the end-on and side-on structures) and monodentate (the backside structure) configurations are localized [Donald et al., Proteins 79, 898–915 (2011)]. These structures are stabilized by the short {sup +}N–H⋯O{sup −} bonds. Their relative stability depends on a force field used in the MD simulations. The side-on structure is the most stable in terms of the OPLS-AA force field. If AMBER ff99SB-ILDN is used, the backside structure is the most stable. Compared with experimental data, simulations using the OPLS all-atom (OPLS-AA) force field describe the stability of the salt bridge structures quite realistically. It decreases in the following order: side-on > end-on > backside. The most stable side-on structure lives several nanoseconds. The less stable backside structure exists a few tenth of a nanosecond. Several short-living species (solvent shared, completely separately solvated ionic groups ion pairs, etc.) are also localized. Their lifetime is a few tens of picoseconds or less. Conformational flexibility of amino acids forming the salt bridge is investigated. The spectral signature of the Arg-Glu salt bridge is the IR-intensive band around 2200 cm{sup −1}. It is caused by the asymmetric stretching vibrations of the {sup +}N–H⋯O{sup −} fragment. Result of the present paper suggests that infrared spectroscopy in the 2000–2800 frequency region may be a rapid and quantitative method for the study of salt bridges in peptides and ionic interactions between proteins. This region is usually not considered in spectroscopic studies of peptides and proteins.

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

Thermoelectric properties of fully hydrogenated graphene: Semi-classical Boltzmann theory

Based on the calculated band structure, the electronic transport coefficients of chair-/boat-like graphane were evaluated by using the semi-classical Boltzmann theory and rigid band model. The maximum value of electrical conductivity for chair (boat)-like graphane of about 1.4 (0.6) × 10{sup 19} (Ωms){sup −1} is achieved at 600 K. The charge carrier concentration and the electrical conductivity linearly increase with increasing the temperature in agreement with the experimental work for graphene. The investigated materials exhibit the highest value of Seebeck coefficient at 300 K. We should emphasize that in the chemical potential between ∓0.125 μ(eV) the investigated materials exhibit minimum value of electronic thermal conductivity, therefore, maximum efficiency. As the temperature increases, the electronic thermal conductivity increases exponentially, in agreement with the experimental data of graphene. We also calculated the power factor of chair-/boat-like graphane at 300 and 600 K as a function of chemical potential between ∓0.25 μ(eV)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Decoherence and energy relaxation in the quantum-classical dynamics for charge transport in organic semiconducting crystals: An instantaneous decoherence correction approach

We explore an instantaneous decoherence correction (IDC) approach for the decoherence and energy relaxation in the quantum-classical dynamics of charge transport in organic semiconducting crystals. These effects, originating from environmental fluctuations, are essential ingredients of the carrier dynamics. The IDC is carried out by measurement-like operations in the adiabatic representation. While decoherence is inherent in the IDC, energy relaxation is taken into account by considering the detailed balance through the introduction of energy-dependent reweighing factors, which could be either Boltzmann (IDC-BM) or Miller-Abrahams (IDC-MA) type. For a non-diagonal electron-phonon coupling model, it is shown that IDC tends to enhance diffusion while energy relaxation weakens this enhancement. As expected, both the IDC-BM and IDC-MA achieve a near-equilibrium distribution at finite temperatures in the diffusion process, while in the Ehrenfest dynamics the electronic system tends to infinite temperature limit. The resulting energy relaxation times with the two kinds of factors lie in different regimes and exhibit different dependences on temperature, decoherence time, and electron-phonon coupling strength, due to different dominant relaxation processes.

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

Second order classical perturbation theory for the sticking probability of heavy atoms scattered on surfaces

A second order classical perturbation theory is developed to calculate the sticking probability of a particle scattered from an uncorrugated thermal surface. An analytic expression for the temperature dependent energy loss of the particle to the surface is derived by employing a one-dimensional generalized Langevin equation. The surface temperature reduces the energy loss, since the thermal surface transfers energy to the particle. Using a Gaussian energy loss kernel and the multiple collision theory of Fan and Manson [J. Chem. Phys. 130, 064703 (2009)], enables the determination of the fraction of particles trapped on the surface after subsequent momentum reversals of the colliding particle. This then leads to an estimate of the trapping probability. The theory is tested for the model scattering of Ar on a LiF(100) surface. Comparison with numerical simulations shows excellent agreement of the analytical theory with simulations, provided that the energy loss is determined by the second order perturbation theory.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Temperature-dependent kinetic measurements and quasi-classical trajectory studies for the OH{sup +} + H{sub 2}/D{sub 2} → H{sub 2}O{sup +}/HDO{sup +} + H/D reactions

We have measured the temperature-dependent kinetics for the reactions of OH{sup +} with H{sub 2} and D{sub 2} using a selected ion flow tube apparatus. Reaction occurs via atom abstraction to result in H{sub 2}O{sup +}/HDO{sup +} + H/D. Room temperature rate coefficients are in agreement with prior measurements and resulting temperature dependences are T{sup 0.11} for the hydrogen and T{sup 0.25} for the deuterated reactions. This work is prompted in part by recent theoretical work that mapped a full-dimensional global potential energy surface of H{sub 3}O{sup +} for the OH{sup +} + H{sub 2} → H + H{sub 2}O{sup +} reaction [A. Li and H. Guo, J. Phys. Chem. A 118, 11168 (2014)], and reported results of quasi-classical trajectory calculations, which are extended to a wider temperature range and initial rotational state specification here. Our experimental results are in excellent agreement with these calculations which accurately predict the isotope effect in addition to an enhancement of the reaction rate constant due to the molecular rotation of OH{sup +}. The title reaction is of high importance to astrophysical models, and the temperature dependence of the rate coefficients determined here should now allow for better understanding of this reaction at temperatures more relevant to the interstellar medium.

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

Quantum and classical dynamics of water dissociation on Ni(111): A test of the site-averaging model in dissociative chemisorption of polyatomic molecules

Recently, we reported the first highly accurate nine-dimensional global potential energy surface (PES) for water interacting with a rigid Ni(111) surface, built on a large number of density functional theory points [B. Jiang and H. Guo, Phys. Rev. Lett. 114, 166101 (2015)]. Here, we investigate site-specific reaction probabilities on this PES using a quasi-seven-dimensional quantum dynamical model. It is shown that the site-specific reactivity is largely controlled by the topography of the PES instead of the barrier height alone, underscoring the importance of multidimensional dynamics. In addition, the full-dimensional dissociation probability is estimated by averaging fixed-site reaction probabilities with appropriate weights. To validate this model and gain insights into the dynamics, additional quasi-classical trajectory calculations in both full and reduced dimensions have also been performed and important dynamical factors such as the steering effect are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Ion-Stockmayer clusters: Minima, classical thermodynamics, and variational ground state estimates of Li{sup +}(CH{sub 3}NO{sub 2}){sub n} (n = 1–20)

Structural optimizations, classical NVT ensemble, and variational Monte Carlo simulations of ion Stockmayer clusters parameterized to approximate the Li{sup +}(CH{sub 3}NO{sub 2}){sub n} (n = 1–20) systems are performed. The Metropolis algorithm enhanced by the parallel tempering strategy is used to measure internal energies and heat capacities, and a parallel version of the genetic algorithm is employed to obtain the most important minima. The first solvation sheath is octahedral and this feature remains the dominant theme in the structure of clusters with n ≥ 6. The first “magic number” is identified using the adiabatic solvent dissociation energy, and it marks the completion of the second solvation layer for the lithium ion-nitromethane clusters. It corresponds to the n = 18 system, a solvated ion with the first sheath having octahedral symmetry, weakly bound to an eight-membered and a four-membered ring crowning a vertex of the octahedron. Variational Monte Carlo estimates of the adiabatic solvent dissociation energy reveal that quantum effects further enhance the stability of the n = 18 system relative to its neighbors.

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

Direct measurements of classical and enhanced gradient-aligned cross-field ion flows in a helicon plasma source using laser-induced fluorescence

Direct laser induced fluorescence measurements are shown of cross-field ion flows normal to an absorbing boundary that is aligned parallel to the axial magnetic field in a helicon plasma. We show Langmuir and emissive probe measurements of local density and plasma potential in the same region, as well as floating probe spectra near the boundary. With these measurements, we investigate the influence of ion-neutral collisionality on radial ion transport by varying the ratio of the ion gyro-radius, ρ{sub i}, to the ion-neutral collision length, λ, over the range 0.34 ≤ ρ{sub i}λ{sup −1} ≤ 1.60. Classical drift-diffusion transport along density and potential gradients is sufficient to describe flow profiles for most cases. For two parameter regimes (ρ{sub i}λ{sup −1} = 0.65 and 0.44), low-frequency electrostatic fluctuations (f < 10 kHz) and enhanced cross-field bulk ion flow to the boundary are observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Classical calculation of the equilibrium constants for true bound dimers using complete potential energy surface

The present paper aims at deriving classical expressions which permit calculation of the equilibrium constant for weakly interacting molecular pairs using a complete multidimensional potential energy surface. The latter is often available nowadays as a result of the more and more sophisticated and accurate ab initio calculations. The water dimer formation is considered as an example. It is shown that even in case of a rather strongly bound dimer the suggested expression permits obtaining quite reliable estimate for the equilibrium constant. The reliability of our obtained water dimer equilibrium constant is briefly discussed by comparison with the available data based on experimental observations, quantum calculations, and the use of RRHO approximation, provided the latter is restricted to formation of true bound states only.

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

New insights into the nonadiabatic state population dynamics of model proton-coupled electron transfer reactions from the mixed quantum-classical Liouville approach

In a previous study [F. A. Shakib and G. Hanna, J. Chem. Phys. 141, 044122 (2014)], we investigated a model proton-coupled electron transfer (PCET) reaction via the mixed quantum-classical Liouville (MQCL) approach and found that the trajectories spend the majority of their time on the mean of two coherently coupled adiabatic potential energy surfaces. This suggested a need for mean surface evolution to accurately simulate observables related to ultrafast PCET processes. In this study, we simulate the time-dependent populations of the three lowest adiabatic states in the ET-PT (i.e., electron transfer preceding proton transfer) version of the same PCET model via the MQCL approach and compare them to the exact quantum results and those obtained via the fewest switches surface hopping (FSSH) approach. We find that the MQCL population profiles are in good agreement with the exact quantum results and show a significant improvement over the FSSH results. All of the mean surfaces are shown to play a direct role in the dynamics of the state populations. Interestingly, our results indicate that the population transfer to the second-excited state can be mediated by dynamics on the mean of the ground and second-excited state surfaces, as part of a sequence of nonadiabatic transitions that bypasses the first-excited state surface altogether. This is made possible through nonadiabatic transitions between different mean surfaces, which is the manifestation of coherence transfer in MQCL dynamics. We also investigate the effect of the strength of the coupling between the proton/electron and the solvent coordinate on the state population dynamics. Drastic changes in the population dynamics are observed, which can be understood in terms of the changes in the potential energy surfaces and the nonadiabatic couplings. Finally, we investigate the state population dynamics in the PT-ET (i.e., proton transfer preceding electron transfer) and concerted versions of the model. The PT-ET results confirm the participation of all of the mean surfaces, albeit in different proportions compared to the ET-PT case, while the concerted results indicate that the mean of the ground- and first-excited state surfaces only plays a role, due to the large energy gaps between the ground- and second-excited state surfaces.

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

State-specific tunneling lifetimes from classical trajectories: H-atom dissociation in electronically excited pyrrole

A trajectory method of calculating tunneling probabilities from phase integrals along straight line tunneling paths, originally suggested by Makri and Miller [J. Chem. Phys. 91, 4026 (1989)] and recently implemented by Truhlar and co-workers [Chem. Sci. 5, 2091 (2014)], is tested for one- and two-dimensional ab initio based potentials describing hydrogen dissociation in the {sup 1}B{sub 1} excited electronic state of pyrrole. The primary observables are the tunneling rates in a progression of bending vibrational states lying below the dissociation barrier and their isotope dependences. Several initial ensembles of classical trajectories have been considered, corresponding to the quasiclassical and the quantum mechanical samplings of the initial conditions. It is found that the sampling based on the fixed energy Wigner density gives the best agreement with the quantum mechanical dissociation rates.

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

Insight into the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} eutectic mixture from classical molecular dynamics: Thermodynamics, structure, and dynamics

We use molecular dynamics simulations to study the thermodynamics, structure, and dynamics of the Li{sub 2}CO{sub 3}–K{sub 2}CO{sub 3} (62:38 mol. %) eutectic mixture. We present a new classical non-polarizable force field for this molten salt mixture, optimized using experimental and first principles molecular dynamics simulations data as reference. This simple force field allows efficient molecular simulations of phenomena at long time scales. We use this optimized force field to describe the behavior of the eutectic mixture in the 900–1100 K temperature range, at pressures between 0 and 5 GPa. After studying the equation of state in these thermodynamic conditions, we present molecular insight into the structure and dynamics of the melt. In particular, we present an analysis of the temperature and pressure dependence of the eutectic mixture’s self-diffusion coefficients, viscosity, and ionic conductivity.

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

Quantum and classical contributions to linear magnetoresistance in topological insulator thin films

Three dimensional topological insulators possess backscattering immune relativistic Dirac fermions on their surface due to nontrivial topology of the bulk band structure. Both metallic and bulk insulating topological insulators exhibit weak-antilocalization in the low magnetic field and linear like magnetoresistance in higher fields. We explore the linear magnetoresistance in bulk insulating topological insulator Bi{sub 2-x}Sb{sub x}Te{sub 3-y}Se{sub y} thin films grown by pulsed laser deposition technique. Thin films of Bi{sub 2-x}Sb{sub x}Te{sub 3-y}Se{sub y} were found to be insulating in nature, which conclusively establishes the origin of linear magnetoresistance from surface Dirac states. The films were thoroughly characterized for their crystallinity and composition and then subjected to transport measurements. We present a careful analysis taking into considerations all the existing models of linear magnetoresistance. We comprehend that the competition between classical and quantum contributions to magnetoresistance results in linear magnetoresistance in high fields. We observe that the cross-over field decreases with increasing temperature and the physical argument for this behavior is explained.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Application of classical thermodynamics to the conductivity in non-polar media

Electrical conductivity in non-polar media is a subject which recently regained interest. If most of experiments and theoretical developments were done more than 50 years ago, new experiments and theories have been recently published. As the electrical conductivity describes, at low field, the equilibrium state of a system, it is natural to apply theories based on equilibrium thermodynamics. In this article, well-established classical thermodynamics and solvations models are applied to recently published data. This enables to get a new insight in intriguing phenomena, such as the linear dependence of the conductivity on the concentration of ionic surfactant and the evaluation of conductivity for the mixture of two miscible fluids, such as alcohol and alcane, which have very different conductivities.

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

Effects of guest atomic species on the lattice thermal conductivity of type-I silicon clathrate studied via classical molecular dynamics

The effects of guest atomic species in Si clathrates on the lattice thermal conductivity were studied using classical molecular dynamics calculations. The interaction between a host atom and a guest atom was described by the Morse potential function while that between host atoms was described by the Tersoff potential. The parameters of the potentials were newly determined for this study such that the potential curves obtained from first-principles calculations for the insertion of a guest atom into a Si cage were successfully reproduced. The lattice thermal conductivities were calculated by using the Green-Kubo method. The experimental lattice thermal conductivity of Ba{sub 8}Ga{sub 16}Si{sub 30} can be successfully reproduced using the method. As a result, the lattice thermal conductivities of type-I Si clathrates, M{sub 8}Si{sub 46} (M = Na, Mg, K, Ca Rb, Sr, Cs, or Ba), were obtained. It is found that the lattice thermal conductivities of M{sub 8}Si{sub 46}, where M is IIA elements (i.e., M = Mg, Ca, Sr, or Ba) tend to be lower than those of M{sub 8}Si{sub 46}, where M is IA elements (i.e., M = Na, K, Rb, or Cs). Those of {sup m}M{sub 8}Si{sub 46}, where m was artificially modified atomic weight were also obtained. The obtained lattice thermal conductivity can be regarded as a function of a characteristic frequency, f{sub c}. That indicates minimum values around f{sub c}=2-4 THz, which corresponds to the center of the frequencies of the transverse acoustic phonon modes associated with Si cages.

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

Fuzzy logic controller versus classical logic controller for residential hybrid solar-wind-storage energy system

Several researches for management of diverse hybrid energy systems and many techniques have been proposed for robustness, savings and environmental purpose. In this work we aim to make a comparative study between two supervision and control techniques: fuzzy and classic logics to manage the hybrid energy system applied for typical housing fed by solar and wind power, with rack of batteries for storage. The system is assisted by the electric grid during energy drop moments. A hydrogen production device is integrated into the system to retrieve surplus energy production from renewable sources for the household purposes, intending the maximum exploitation of these sources over years. The models have been achieved and generated signals for electronic switches command of proposed both techniques are presented and discussed in this paper.

14 SOLAR ENERGY↗