SEARCH · Engineering Papers
Results for “ab initio phasing”
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.
Compressive strain turns 𝑠 ± - into 𝑑-wave pairing in a one-unit-cell La 3 Ni 2 O 7 thin film via substrate-induced hole doping
Motivated by recent reports of ambient-pressure superconductivity in La 3 Ni 2 O 7 films grown on LaSrAlO 4 , we investigate the superconducting instability in a one-unit-cell (1UC) thin film using ab initio and random-phase approximation techniques. Compared to the high-pressure bulk system, the ratio of interlayer 𝑑 3𝑧 2 −𝑟 2 hopping to intralayer 𝑑 𝑥 2 −𝑦 2 hopping is suppressed in the 1UC thin film, and the crystal-field splitting of the 𝑒 𝑔 orbitals is increased. Here, our calculation indicates that spin-fluctuation-driven pairing correlations are weak for the stoichiometric case at ambient pressure, but increase significantly under hole doping. The leading pairing symmetry is also found to change by hole doping. Specifically, we obtain a leading 𝑑 𝑥 2 −𝑦 2 pairing state at moderate hole doping, followed by a 𝑑 𝑥𝑦 state at higher doping. These states are driven by intraband spin-fluctuation scattering within the 𝛾 hole pocket centered around the 𝑀 point, and arise primarily from states in the Ni layer farther from the substrate. These results strongly suggest that the thin-film superconducting samples are hole-doped and that pairing in this system predominantly arises in the layer, as opposed to the interlayer pairing in the pressurized bulk system.
Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La 3 Ni 2 O 7
Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 Ni 2 O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3 Ni 2 O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 Ni 2 O 7 always has lower energy than the single-layer trilayer stacking La 3 Ni 2 O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.
Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure
Motivated by the recently reported signatures of superconductivity in trilayer La 4 Ni 3 O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4 Ni 3 O 10 , analogous to the dimers of La 3 Ni 2 O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 Ni 2 O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3 Ni 2 O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 Ni 3 O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.
Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase
Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.
Integrated design of aluminum-enriched high-entropy refractory B2 alloys with synergy of high strength and ductility
Refractory high-entropy alloys (RHEAs) are promising high-temperature structural materials. Their large compositional space poses great design challenges for phase control and high strength-ductility synergy. The present research pioneers using integrated high-throughput machine learning with Monte Carlo simulations supplemented by ab initio calculations to effectively navigate phase selection and mechanical property predictions, developing single-phase ordered B2 aluminum-enriched RHEAs (Al-RHEAs) demonstrating high strength and ductility. These Al-RHEAs achieve remarkable mechanical properties, including compressive yield strengths up to 1.7 gigapascals, fracture strains exceeding 50%, and notable high-temperature strength retention. They also demonstrate a tensile yield strength of 1.0 gigapascals with a ductility of 9%, albeit with B2 ordering. Furthermore, we identify valence electron count domains for alloy ductility and brittleness with the explanation from density functional theory and provide crucial insights into elemental influence on atomic ordering and mechanical performance. The work sets forth a strategic blueprint for high-throughput alloy design and reveals fundamental principles governing the mechanical properties of advanced structural alloys.
Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]
Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.
Tuning charge carrier dynamics through spacer cation functionalization in layered halide perovskites: an ab initio quantum dynamics study
Dion–Jacobson (DJ) phase two-dimensional (2D) hybrid halide perovskites are promising for environmentally stable optoelectronic device applications due to their attractive photophysical properties, including long charge carrier lifetime and partially suppressed non-radiative losses. However, the atomistic details of the structure–property relationship are severely limited, which substantially restricts the strategic materials designed for these layered halide perovskites (LHPs). Here, we combine nonadiabatic molecular dynamics and time-domain density functional theory to understand the effects of spacer cation functionalization on the ground and excited-state charge carrier dynamics in DJ phase LHPs. Our in-depth study reveals that the fluorination of spacer cations considerably restricts the thermal motions of LHPs at ambient conditions. The compact structure weakens the overall electron–phonon interactions and reduces the thermal fluctuations of the band edges in real-time. These dynamic modifications partially mitigate the non-radiative recombination, prolonging the lifetime of photogenerated charge carriers. The suppressed carrier loss mechanism strongly suggests that the fluorinated spacer cation-based LHPs would exhibit enhanced performance as optoelectronic materials. These systematic simulations of excited state carrier dynamics elaborate that chemically viable functionalization of the spacer cations is a robust approach to improve the photophysical properties of 2D halide perovskites strategically. These insights also guide us to propose a few potential strategies to design highly beneficial spacer cations of LHPs that can be introduced in next-generation optoelectronic devices.
Crystallographic map: A general lattice and basis formalism enabling efficient and discretized exploration of crystallographic phase space
Three-dimensional lattices are fundamental to solid-state physics. The description of a lattice with an atomic basis constitutes the necessary information to predict solid phase properties and evolution. Here, we present an algorithm for systematically exploring crystallographic phase space. Further, coupled with ab initio techniques, such as density functional theory, this algorithm offers an approach for exploring and tuning materials behavior, with a broad range of potential applications: particularly martensitic phase transformations and materials stability.
Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study
The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.
Floquet band engineering and topological phase transitions in 1T’ transition metal dichalcogenides
Using ab initio tight-binding approaches, we investigate Floquet band engineering of the 1T' phase of transition metal dichalcogenides (MX 2 , M = W, Mo and X = Te, Se, S) monolayers under the irradiation with circularly polarized light. Furthermore, our first principles calculations demonstrate that light can induce important transitions in the topological phases of this emerging materials family. For example, upon irradiation, Te-based MX2 undergoes a phase transition from quantum spin Hall (QSH) semimetal to time-reversal symmetry broken QSH insulator with a nontrivial band gap of up to 92.5 meV. On the other hand, Se- and S-based MX 2 undergoes the topological phase transition from the QSH effect to the quantum anomalous Hall effect and into trivial phases with increasing light intensity. From a general perspective, this theoretical work brings further insight into non-equilibrium topological systems.
Phase transitions and potential ferroelectricity in noncentrosymmetric KNaNbOF 5
Increasing the number of known ferroelectrics requires expanding the chemical design space and mechanisms producing ferroelectricity. To that end, we examine the displacive, noncentrosymmetric-to-centrosymmetric phase transition in the oxyfluoride KNaNbOF 5 using ab initio calculations and Landau theory. Here we predict an intermediate Pnma phase occurs in the transition from the known low- and high-temperature phases, Pna2 1 and Cmcm, respectively, which we subsequently confirm using high-resolution and in situ 19 F magic-angle spinning nuclear magnetic resonance. Using the intermediate structure, we then calculate the monodomain polarization switching barrier in KNaNbOF 5 to be ~ 93 meV per formula unit, comparable to existing ferroelectrics. The reversal of the polarization is due to uncompensated antipolar displacements of oxygen and fluorine and does not require oxygen and fluorine site exchange as in other nonswitchable heteroanionic materials, which makes KNaNbOF 5 a viable oxyfluoride for a switchable electric polarization.
Structural and magnetic transitions of CoFeMnNiAl high-entropy alloys caused by composition and annealing
Tailoring magnetic behavior by composition and annealing is an effective way. The phase transition from the face-centered-cubic (FCC) to B2 structure influenced by the Al content on the CoFeMnNiAl x alloys leads to the enhanced magnetization. However, the excessive Al addition in the B2-structured alloys impairs the magnetization. The phase transformation caused by annealing also indicates that the saturation magnetization of the alloy depends on the volume fraction of the B2 phase, and a higher content of the B2 phase is beneficial to the large saturation magnetization. Ab initio calculations are used to explain the magnetic behavior of the present HEAs and the effects of phase structures on magnetic characteristics. The phase stability of the CoFeMnNiAl x alloy is thoroughly studied by both the existing phase-formation empirical criteria and annealing. Here, the existing empirical criteria used to predict the phase formation are not universal. However, combining the ΔH mix -δ relation and the φ value, whether the single disordered solid-solution phase is formed can be preliminarily predicted.
Liquid-Vapor Coexistence and Critical Point of Mg 2 SiO 4 From Ab Initio Simulations
Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.
Ab Initio Simulations of Tritium Diffusion in Al-Rich Iron Aluminide Coating Phases
Surface coatings of steels used in extreme conditions and corrosive environments generally aim to provide protection and increased durability. In the case of tritium-producing burnable absorber rods (TPBARs) used in nuclear reactors, a 316 stainless steel has been coated with Al. Scanning transmission electron microscopy (STEM) characterization of the coating found three Al-rich (>60 atom % Al) iron aluminide alloys identified as hexagonal FeNiAl 5 , monoclinic Fe 4 Al 13 , and orthorhombic Fe 2 Al 5 . Density functional theory simulations using nudged elastic band have been performed to investigate the diffusion of interstitial tritium in each Al-rich iron aluminide phase. While FeNiAl 5 and Fe 4 Al 13 can be viewed as the stacking of two layers, the structural peculiarity of Fe 2 Al 5 is that channels of variable Al vacancy content are present along the c-axis. Therefore, three stoichiometries for Fe 2 Al x phase, namely, Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 , have been considered to evaluate the impact of Al vacancy concentration on tritium diffusion behavior. Altogether, we found that at 600 K, tritium diffusion decreases from a faster rate in the channels of Fe 2 Al x phases (D T ≤ 10 –11 m 2 ·s –1 ) to Fe 4 Al 13 (D T ≈ 10–12 m2·s–1), and finally in FeNiAl5 (D T ≈ 10 –13 m 2 ·s –1 ). Here, we also find that interstitial tritium generally diffuses faster in Fe–Al coating phases than in the tritium breeding material γ-LiAlO 2 (D T ≈ 10 –14 m 2 ·s –1 ) but slightly slower than in 316 stainless steel (D T ≈ 10 –10 m 2 ·s –1 ).
Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data
Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.
Structures of elemental potassium at terapascal pressures
Here we investigate the pressure-temperature phase diagram of elemental potassium (K) up to multiterapascal (TPa) pressures using ab initio random structure searching (AIRSS), discovering eight structural phase transitions beyond the already known double hexagonal close-packed (dhcp) structure. Starting at 1.15 TPa, K transitions from the dhcp structure and passes through a variety of close-packed hexagonal and trigonal phases (dhcp → P6 3 /mmc → P3m1 → P$\bar{6}m2$), which differ only in their stacking sequence along the c axis. At 2.55 TPa, K adopts the hexagonal close-packed (hcp) structure, and at 7.21 TPa transitions into a complex orthorhombic phase of Fddd symmetry with 64 atoms in its (conventional) unit cell, before assuming an Ibam structure at 8.60 TPa, and, eventually, the face centered cubic (fcc) structure at 23.6 TPa, which persists well into the petapascal (PPa) regime. Further, we calculate the full pressure-temperature dependence of the melting line of K through extensive molecular dynamics simulations. We study the evolution of the bonding topology of K with pressure, finding that K passes through two quasi-molecular phases featuring diatomic pairs (the Fddd and Ibam structures), before ultimately becoming a high pressure electride (HPE) in the fcc phase. The electron-phonon and superconducting properties of K at these extreme pressures are also investigated, where we find the critical temperature T c rises to a maximum of between 7.65 K and 15.70 K in the hcp phase, before eventually decreasing to essentially zero in the fcc phase. Our results fully elucidate the structural and electronic behavior of K under the most extreme conditions, and provide new case studies for multi-TPa dynamic compression experiments.
Ab Initio Calculations on the Ground and Excited Electronic States of Thorium–Ammonia, Thorium–Aza-Crown, and Thorium–Crown Ether Complexes
Positively charged metal–ammonia complexes are known to host peripheral, diffuse electrons around their molecular skeleton. The resulting neutral species form materials known as expanded or liquid metals. Alkali, alkaline earth, and transition metals have been investigated previously in experimental and theoretical studies of both the gas and condensed phase. This work is the first ab initio exploration of an f-block metal–ammonia complex. The ground and excited states are calculated for Th0–3+ complexes with ammonia, crown ethers, and aza-crown ethers. For Th3+ complexes, the one valence electron Th populates the metal’s 6d or 7f orbitals. For Th0–2+, the additional electrons prefer occupation of the outer s- and p-type orbitals of the complex, except Th(NH3)10, which uniquely places all four electrons in outer orbitals of the complex. Although thorium coordinates up to ten ammonia ligands, octa-coordinated complexes are more stable. Crown ether complexes have a similar electronic spectrum to ammonia complexes, but excitations of electrons in the outer orbitals of the complex are higher in energy. Aza-crown ethers disfavor the orbitals perpendicular to the crowns, attributed to the N-H bonds pointing along the plane of the crowns.