Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SOC”

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 973 records · Page 54

Nearly massless Dirac fermions hosted by Sb square net in BaMnSb 2

Layered compounds AMnBi 2 (A = Ca, Sr, Ba, or rare earth element) have been established as Dirac materials. Dirac electrons generated by the two-dimensional (2D) Bi square net in these materials are normally massive due to the presence of a spin-orbital coupling (SOC) induced gap at Dirac nodes. Here we report that the Sb square net in an isostructural compound BaMnSb 2 can host nearly massless Dirac fermions. We observed strong Shubnikov-de Haas (SdH) oscillations in this material. From the analyses of the SdH oscillations, we find key signatures of Dirac fermions, including light effective mass (~0.052 m 0 ; m 0 , mass of free electron), high quantum mobility (1280 cm 2 V –1 S –1 ) and a π Berry phase accumulated along cyclotron orbit. Compared with AMnBi 2 , BaMnSb 2 also exhibits much more significant quasi two-dimensional (2D) electronic structure, with the out-of-plane transport showing nonmetallic conduction below 120 K and the ratio of the out-of-plane and in-plane resistivity reaching ~670. Additionally, BaMnSb 2 also exhibits a G-type antiferromagnetic order below 283 K. The combination of nearly massless Dirac fermions on quasi-2D planes with a magnetic order makes BaMnSb 2 an intriguing platform for seeking novel exotic phenomena of massless Dirac electrons.

36 MATERIALS SCIENCE↗

Dynamics of Energy Transfer and Soft-Landing in Collisions of Protonated Dialanine with Perfluorinated Self-Assembled Monolayer Surfaces

Chemical dynamics simulations are reported which provide atomistic details of collisions of protonated dialanine, ala2-H+, with a perfluorinateted octanethiolate self-assembled monolayer (F-SAM ) surface. The simulations are performed at collisions energy Ei of 5.0, 13.5, 22.5, 30.00, and 70 eV, and incident angles 0o 0 (normal) and grazing 45o. Excellent agreement with experiment (J. Am. Chem. Soc. 2000, 122, 9703-9714) is found for both the average fraction and distribution of the collision energy transferred to the ala2-H+ internal degrees of freedom. The dominant pathway for this energy transfer is to ala2-H+ vibration, but for Ei = 5.0 eV ~20% of the energy transfer is to ala2-H+ rotation. Energy transfer to ala2-H+ rotation decreases with increase in Ei and becomes negligible at high Ei. Three types of collisions are observed in the simulations: i.e. those for which ala2-H+ (1) directly scatters off the F-SAM surface; (2) sticks/physisorbs on//in the surface, but desorbs within the 10 ps numerical integration of the simulations; and (3) remains trapped (i.e. soft-landed) on/in the surface when the simulations are terminated. Penetration of the F-SAM by ala2-H+ is important for the latter two types of events. The trapped trajectories are expected to have relatively long residence times on the surface, since a previous molecular dynamics simulation (J. Phys. Chem. B 2014, 118, 5577-5588) shows that thermally accommodated ala2-H+ ions have an binding energy with the F-SAM surface of at least ~15 kcal/mol.

Pratihar, Subha↗

Gas phase vibrational spectroscopy of the protonated water pentamer: the role of isomers and nuclear quantum effects

We use cryogenic ion trap vibrational spectroscopy to study the structure of the protonated water pentamer, H + (H 2 O) 5 , and its fully deuterated isotopologue, D + (D 2 O) 5 , over nearly the complete infrared spectral range (220-4000 cm-1) in combination with harmonic and anharmon ic electronic structure calculations as well as RRKM modelling. Isomer-selective IR–IR double-resonance measurements on the H + (H 2 O) 5 isotopologue establish that the spectrum is due to a single constitutional isomer, thus discounting the recent analysis of the band pattern in the context of two isomers based on AIMD simulations < W. Kulig and N. Agmon, Phys. Chem. Chem. Phys., 2014, 16, 4933–4941 >. The evolution of the persistent bands in the D + (D 2 O) 5 cluster allows the assignment of the fundamentals in the spectra of both isotopologues, and the simpler pattern displayed by the heavier isotopologue is consistent with the calculated spectrum for the branched, Eigen-based structure originally proposed < J.-C. Jiang, et al., J. Am. Chem. Soc., 2000, 122, 1398–1410 >. This pattern persists in the vibrational spectra of H + (H 2 O) 5 in the temperature range from 13 K up to 250 K. The present study also underscores the importance of considering nuclear quantum effects in predicting the kinetic stability of these isomers at low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Is there such a thing as a molecular organic alloy?

In inorganic materials, an alloy is a mixture of two or more substances that generally exhibits electronic and/or physical properties that differ from those of its constituents. In organic systems, the formation of a “molecular alloy” comprising mixtures of molecular organic materials has also been proposed. We test the validity of this concept via the study of the optoelectronic properties of a ternary system that has previously been identified to form a molecular acceptor alloy, namely a blend of a poly(3-hexylthiophene) (P3HT) donor, with two acceptors indene-C 60 bisadduct (ICBA) and phenyl-C 61 -butyric acid methyl ester (PC 61 BM) [R. A. Street, et al., J. Am. Chem. Soc., 2013, 135, 986–989]. Using photoelectron spectroscopy, we find that the ICBA:PC 61 BM blend shows the same highest occupied molecular orbital and exciton energies as that of ICBA, indicating the absence of a new exciton state in the blend. Furthermore, charge transfer state spectra of ternary blends are found to comprise a simple linear superposition of the corresponding binaries. From these results, no evidence of new, emergent electronic states is found to support the existence of a molecular alloy in this system. To our knowledge there is as yet no clear evidence of the existence of an alloy in any organic semiconductor system. Finally, we discuss the criteria that should be met by a molecular organic alloy and procedures needed for their unambiguous identification.

36 MATERIALS SCIENCE↗

High stability SrTi 1-x Fe x O 3-δ electrodes for oxygen reduction and oxygen evolution reactions

We report that Sr(Ti 1-x Fe x )O 3-δ (STF) has recently been explored as an oxygen electrode for solid oxide electrochemical cells (SOCs). Model thin film electrode studies show oxygen surface exchange rates that generally improve with increasing Fe content when x < 0.5, and are comparable to the best Co-containing perovskite electrode materials. Recent results on porous electrodes with the specific composition Sr(Ti 0.3 Fe 0.7 )O 3-δ show excellent electrode performance and stability, but other compositions have not been tested. Here we report results for porous electrodes with a range of compositions from x = 0.5 to 0.9. The polarization resistance decreases with increasing Fe content up to x = 0.7, but increases for further increases in x. This results from the interaction of two effects – the oxygen solid state diffusion coefficient increases with increasing x, but the electrode surface area and surface oxygen exchange rate decrease due to increased sinterability and Sr surface segregation for the Fe-rich compositions. Symmetric cells showed no degradation during 1000 h life tests at 700 °C even at a current density of 1.5 A cm -2 , showing that all the STF electrode compositions worked stably in both fuel cell mode and electrolysis modes. Lastly, the excellent stability may be explained by X-ray Photoelectron Spectroscopy (XPS) results showing that the amount of surface segregated Sr did not change during the long-term testing, and by relatively low polarization resistances that help avoid electrode delamination.

36 MATERIALS SCIENCE↗

Enhancement of Ni–(Y 2 O 3 ) 0.08 (ZrO 2 ) 0.92 fuel electrode performance by infiltration of Ce 0.8 Gd 0.2 O 2–δ nanoparticles

This paper addresses the use of Ce 0.8 Gd 0.2 O 2–δ (GDC) infiltration into the Ni–(Y 2 O 3 ) 0.08 (ZrO 2 ) 0.92 (YSZ) fuel electrode of solid oxide cells (SOCs) for improving their electrochemical performance in fuel cell and electrolysis operation. Although doped ceria infiltration into Ni–YSZ has recently been shown to improve the electrode performance and stability, the mechanisms defining how GDC impacts electrochemical characteristics are not fully delineated. Furthermore, the electrochemical characteristics have not yet been determined over the full range of conditions normally encountered in fuel cell and electrolysis operation. Here we present a study of both symmetric and full cells aimed at understanding the electrochemical mechanisms of GDC-modified Ni–YSZ over a wide range of fuel compositions and temperatures. Single-step GDC infiltration at an appropriate loading substantially reduced the polarization resistance of Ni–YSZ electrodes in electrolyte-supported cells, as measured using electrochemical impedance spectroscopy (EIS) at various temperatures (600–800 °C) in a range of H 2 O–H 2 mixtures (3–90 vol% H 2 O). Fuel-electrode-supported cells had significant concentration polarization due to the thick Ni–YSZ supports. Here, a distribution of relaxation times approach is used to develop a physically-based electrochemical model; the results show that GDC reduces the reaction resistance associated with three-phase boundaries, but also appears to improve oxygen transport in the electrode. Increasing the H 2 O fraction in the H 2 –H 2 O fuel mixture reduced both the three-phase boundary resistance and the gas diffusion resistance for Ni–YSZ; with GDC infiltration, the electrode resistance showed less variation with fuel composition. GDC infiltration improved the performance of fuel-electrode-supported full cells, which yielded a maximum power density of 2.28 W cm –2 in fuel cell mode and an electrolysis current density at 1.3 V of 2.22 A cm –2 , both at 800°C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A redox-active organic cation for safer high energy density Li-ion batteries

Ni-rich layered cathode materials are at the forefront to be deployed in high energy density Li-ion batteries for the automotive market. However, the intrinsic poor structural and interfacial stability during overcharging could trigger violent thermal failure, which severely limits their wide application. To protect the Ni-rich cathode from overcharging, we firstly report a redox-active cation, thioether-substituted diaminocyclopropenium, as an electrolyte additive to limit the cell voltage within the safe value during overcharging. The organic cation demonstrates a record-breaking electrochemical reversibility at ~4.55 V versus Li + /Li and solubility (0.5 M) in carbonate-based electrolyte. The protection capability of the additive was explored in two cell chemistries: a LiNi 0.8 Co 0.15 Al 0.05 O 2 /graphite cell and a LiNi 0.8 Co 0.15 Al 0.05 O 2 /silicon–graphene cell with areal capacities of ~2.2 mA h cm -2 and ~3 mA h cm -2 , respectively. With 0.2 M addition, the LiNi 0.8 Co 0.15 Al 0.05 O 2 /graphite cell survived 54 cycles at 0.2C with 100% overcharge. Moreover, the cell can carry an utmost 4.4 mA cm -2 (2C) with 100% overcharge and a maximum capacity of 7540% SOC at 0.2C.

25 ENERGY STORAGE↗

Simultaneous enhancement of thermally activated delayed fluorescence and photoluminescence quantum yield via homoconjugation

A critical challenge facing thermally activated delayed fluorescence (TADF) is to facilitate rapid and efficient electronic transitions while ensuring a narrow singlet–triplet energy gap (ΔE ST ) in a single luminophore. We present a TADF-active iptycene that clearly demonstrates that homoconjugation can be harnessed as a viable design strategy towards answering this challenge. A homoconjugated analogue of an established quinoxaline-based TADF luminophore has been produced by fusing three of these luminophores together across a shared triptycene core. Homoconjugation was confirmed by electrochemistry, and as a direct consequence of this phenomenon we observed synergistic improvements to photoluminescence quantum yield (Φ PL ), radiative rate of singlet decay ($k^{S}_{r}$), delayed fluorescence lifetime (τ TADF ), and rate of reverse intersystem crossing (k rISC ), all while narrowing the ΔE ST . The enhancement is rationalised with TD-DFT calculations including spin–orbit coupling (SOC). A facile synthesis, the ubiquity of the pyrazine motif in state-of-the-art TADF materials of all colours, and the extent of the overall performance enhancement leads to a great potential for generality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-polarization anisotropy controlled by bending in tungsten diselenide nanoribbons and tunable excitonic states

A WSe 2 monolayer shows many interesting properties due to its spin-orbit coupling induced spin splitting in bands around the Fermi level and the spin-valley configuration. The orientation of the spin polarization in the relevant bands is crucial for the nature of exciton states and the optical valley selectivity. In this work, we studied the WSe 2 nanoribbons under different mechanical bending curvatures and electron/hole doping with density functional theory and their optical absorption and excitonic states with many-body perturbation GW and BSE (Bethe-Salpeter equation) methods. We found that the WSe 2 nanoribbons can exhibit an enhanced SOC effect and a spatially varying spin polarization in bands around the Fermi level under bending. The spin-polarization can show an anisotropy (or asymmetry) in those nearly degenerate bands, leading to a controllable magnetism via bending and electron/hole doping to the nanoribbons, suggesting a potential application in compact and controllable magnetic nanodevices and spintronics. The optical absorption spectrum of the nanoribbon presents a large tunability with bending within the near infrared region of about 0.4 to 1.5 eV, showing an enhanced absorption at a large bending condition. Finally, the exciton states generally show mixed or various spin configuration in the electron and hole pairs that are controlled by bending, potentially useful for applications in spin-based quantum information processes.

36 MATERIALS SCIENCE↗

The role of oxygen in automotive grade lithium-ion battery cathodes: an atomistic survey of ageing

The rising demand for high-performance lithium-ion batteries, pivotal to electric transportation, hinges on key materials like the Ni-rich layered oxide LiNi x Co y Al z O 2 (NCA) used in cathodes. The present study investigates the redox mechanisms, with particular focus on the role of oxygen in commercial NCA electrodes, both fresh and aged under various conditions (aged cells have performed >900 cycles until a cathode capacity retention of ~80%). Our findings reveal that oxygen participates in charge compensation during NCA delithiation, both through changes in transition metal (TM)–O bond hybridization and formation of partially reversible O 2 , the latter occurs already below 3.8 V vs. Li/Li + . Aged NCA material undergoes more significant changes in TM–O bond hybridization when cycling above 50% SoC, while reversible O 2 formation is maintained. Nickel is found to be redox active throughout the entire delithiation and shows a more classical oxidation state change during cycling with smaller changes in the Ni–O hybridization. By contrast, Co redox activity relies on a stronger change in Co–O hybridization, with only smaller Co oxidation state changes. The Ni–O bond displays an almost twice as large change in its bond length on cycling as the Co–O bond. The Ni–O 6 octahedra are similar in size to the Co–O 6 octahedra in the delithiated state, but are larger in the lithiated state, a size difference that increases with battery ageing. These contrasting redox activities are reflected directly in structural changes. The NCA material exhibits the formation of nanopores upon ageing, and a possible connection to oxygen redox activity is discussed. The difference in interaction of Ni and Co with oxygen provides a key understanding of the mechanism and the electrochemical instability of Ni-rich layered transition metal oxide electrodes. Our research specifically highlights the significance of the role of oxygen in the electrochemical performance of electric-vehicle-grade NCA electrodes, offering important insights for the creation of next-generation long-lived lithium-ion batteries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxygen redox activities governing high-voltage charging reversibility of Ni-rich layered cathodes

A schematic illustration indicative of the improved structural reversibility through stabilizing oxidized oxygen by Al doping. The chemical reactions and phase transitions at high voltages are generally considered to determine the electrochemical properties of high-voltage layered cathodes such as Ni-rich rhombohedral oxides. Even if significantly higher SOCs (states-of-charge) are utilized above the capability of transition metal redox (primarily Ni and Co), the effect of oxygen redox on Ni-rich rhombohedral oxides still looks mysterious thereby necessitating research that can clarify the relationship between redox reactions and phase transitions. Here, we performed a comprehensive and comparative study of the cationic and anionic redox reactions, as well as the structural evolution of a series of commercial Ni-rich layered oxides with and without Al doping. We combined the results from X-ray spectroscopy, operando electrochemical mass spectrometry, and neutron diffraction with electrochemical properties and thereby revealed the different oxygen redox activities associated with structural and electrochemical degradations. We reveal that Al doping suppresses the irreversible oxygen release but enhances the lattice oxygen oxidization. With this modulated oxygen redox activity, the Ni-rich layered oxides’ notorious H2–H3 structural phase transition becomes highly reversible. Our findings disentangle the different oxygen redox activities during high-voltage cycling and clarify the role of dopants in the Ni-rich layered oxides in terms of structural and electrochemical stability, shedding light on the future direction of optimizing layered cathode materials for safer high energy-density secondary batteries.

Lee, Gi-Hyeok↗

Aging matrix visualizes complexity of battery aging across hundreds of cycling protocols

To reliably deploy lithium-ion batteries, a fundamental understanding of cycling aging behavior is critical. Battery aging consists of complex and highly coupled phenomena, making it challenging to develop a holistic interpretation. In this work, we generate a diverse battery cycling dataset with a broad range of degradation trajectories, consisting of 359 high energy density commercial Li(Ni,Co,Al)O 2 /graphite + SiO x cylindrical 21 700 cells cycled across 207 unique cycling protocols. We consolidate aging via 16 mechanistic state-of-health (SOH) metrics, including cell-level performance metrics, electrode-specific capacities/state-of-charges (SOCs), and aging trajectory metrics. We develop a framework using interpretable machine learning and explainable features to generate an aging matrix that visually deconvolutes the complex battery degradation behavior. This generalizable data-driven mechanistic framework simplifies the complex interplay between cycling conditions, degradation modes, and SOH, acting as a hypothesis-generation tool to aid battery users in identifying key degradation regimes for further study and experimentation.

25 ENERGY STORAGE↗

Effects of catholyte aging on high-nickel NMC cathodes in sulfide all-solid-state batteries

Sulfide solid-state electrolytes (SSEs) in all-solid-state batteries (SSBs) are recognized for their high ionic conductivity and inherent safety. The LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) cathode offers a high thermodynamic potential of approximately 3.8 V vs. Li/Li + and a theoretical specific capacity of 200 mA h g −1 . However, the practical utilization of NMC811 in sulfide SSBs faces significant interfacial challenges. The oxidation instability of sulfide solid electrolytes against NMC811 and the formation of the cathode electrolyte interphase (CEI) during cycling lead to degradation and reduced cell performance. Volumetric changes in NMC during lithiation and de-lithiation can also cause detachment from sulfide electrolytes or internal particle cracking. Despite extensive galvanostatic cycling studies to address the issues, the calendar life of sulfide SSBs remains poorly understood. Here, we systematically studied the effects of four different catholytes on the calendar aging of LiNbO 3 (LNO)-coated NMC811, including Li 6 PS 5 Cl (LPSCl), Li 3 InCl 6 –Li 6 PS 5 Cl (LIC–LPSCl), Li 3 YCl 6 –Li 6 PS 5 Cl (LYC–LPSCl), and Li 10 GeP 2 S 12 (LGPS). Our results indicate that LPSCl provides optimal capacity retention when stored at high state-of-charge (SOC) at room temperature, but the LIC–LPSCl cathode shows significant capacity degradation and chemical incompatibility. We also established an effective electrochemical calendar aging testing protocol to simulate daily usage, enabling quick inference of the calendar life of SSBs. In conclusion, this new testing approach accelerates materials selection strategies for high-nickel NMC composite cathodes in sulfide SSBs.

25 ENERGY STORAGE↗

Comparison of Ce( iv )/Th( iv )-alkynyl complexes and observation of a trans -influence ligand series for Ce( iv )

Organometallic cerium(iv) complexes have been challenging to isolate and characterize due to the strongly oxidizing nature of the cerium(iv) cation. Herein, we report two cerium(iv) alkynyl complexes, [Ce(TriNOx)(C[triple bond, length as m-dash]C-SiMe3)] (1-CeTMS) and [Ce(TriNOx)(C[triple bond, length as m-dash]C-Ph)] (1-CePh) (TriNOx3- = tris(2-tert-butylhydroxylaminato)benzylamine), that include terminal alkyne moieties. The isostructural thorium analogue [Th(TriNOx)(C[triple bond, length as m-dash]C-SiMe3)] (1-ThTMS) was also synthesized and compared with 1-CeTMS in bond distance, 13C-NMR spectra, vibrational spectra and electronic structure. The Ce-C bond distances were 2.501(3) Å for 1-CePh and 2.513(5) Å for 1-CeTMS on the shorter end of the few reported CeIV-C single bonds (2.478(3)-2.705(2) Å), possibly indicating significant Ce 5d- and 4f-orbital involvement. 13C-NMR spectroscopy was also consistent with Ce-C covalency, with significantly deshielded resonances ranging from 185-213 ppm. Such 13C-NMR shifts demonstrate a strong influence from spin-orbit coupling (SOC) effects, corroborated by computational studies. Raman analysis showed ν C[triple bond, length as m-dash]C stretching frequencies of 2000 cm-1 (1-CeTMS) and 2052 cm-1 (1-CePh), indicating the cerium(iv)-alkynyl interaction, compared to the parent HC[triple bond, length as m-dash]CPh (IR = 2105 cm-1 and Raman = 2104 cm-1). L3-edge X-ray absorption measurements revealed a predominant Ce(iv) electronic configuration, and magnetic measurements revealed temperature-independent paramagnetism. Electrochemical studies similarly revealed the electron donating ability of the alkynyl ligands, stronger than either fluoride or imido ligands for the Ce(iv)(TriNOx)-framework, with a cerium(iv/iii) reduction potential of E pc = -1.58 to -1.66 V vs. Fc/Fc+. Evidence for a trans-influence has been observed by evaluating a series including previously reported [CeIV(TriNOx)X]+/0 complexes with axial ligands X = THF, I-, Br-, Cl-, F-, -C[triple bond, length as m-dash]C-Ph, -C[triple bond, length as m-dash]C-SiMe3, -NH(3,5-(CF3)2-Ar), -OSiPh3, -N(M(L))(3,5-(CF3)2-Ar) [M(L) = Li(TMEDA), K(DME)2 or Cs(2,2,2-crypt)]. These data stand in contrast with previous reports of an inverse trans-influence at cerium(iv) and point to differences in involvement of cerium 4f- versus 5d-orbitals in the electronic structures of the complexes.

Yang, Qiaomu↗

Presymplectic current and the inverse problem of the calculus of variations

The inverse problem of the calculus of variations asks whether a given system of partial differential equations (PDEs) admits a variational formulation. We show that the existence of a presymplectic form in the variational bicomplex, when horizontally closed on solutions, allows us to construct a variational formulation for a subsystem of the given PDE. No constraints on the differential order or number of dependent or independent variables are assumed. The proof follows a recent observation of Bridges, Hydon, and Lawson [Math. Proc. Cambridge Philos. Soc. 148(01), 159–178 (2010)] and generalizes an older result of Henneaux [Ann. Phys. 140(1), 45–64 (1982)] from ordinary differential equations (ODEs) to PDEs. Uniqueness of the variational formulation is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theoretical study of Al{sub n}V{sup +} clusters and their interaction with Ar

Recently, it has been experimentally elucidated whether a V impurity in Al{sub n}V{sup +} clusters occupies an external or an internal site by studying their interaction with argon as a function of cluster size [S. M. Lang, P. Claes, S. Neukermans, and E. Janssens, J. Am. Soc. Mass Spectrom. 22, 1508 (2011)]. In the work presented here we studied, by means of density functional theoretic calculations, the structural and electronic properties of Al{sub n}V{sup +} clusters with n = 14–21 atoms, as well as the adsorption of a single Ar atom on them. For n < 17 the lowest energy structure of Al{sub n}V{sup +} is related to that of the pure Al {sub n+1}{sup +} cluster with the V atom substituting a surface Al atom. For n ⩾ 17 the V impurity becomes embedded in the cluster, in agreement with the experimental results, and the clusters adopt a fcc-like structure instead of the icosahedral-like skeleton of pure Al {sub n+1}{sup +}. We have studied the binding energy per atom, the second energy difference, and the V and Al atom separation energies, in comparison with those of pure Al {sub n+1}{sup +}. We also studied the adsorption of atomic Ar on endohedral and exohedral V doped clusters. The optimized Ar adsorption geometries are formed with Ar on top of a surface atom (V for n < 17, and Al for n ⩾ 17) without noticeable structural distortion of the host cluster. At the critical size (n = 17) of the exohedral-endohedral transition, the calculated Ar adsorption energy exhibits a drop and the Ar-cluster distance increases drastically, indicating that Ar becomes physisorbed rather than chemisorbed. All these results confirm the assumptions made by the experimentalists when interpreting their measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density functional theory and evolution algorithm calculations of elastic properties of AlON

Different models for aluminum oxynitride (AlON) were calculated using density functional theory and optimized using an evolutionary algorithm. Evolutionary algorithm and density functional theory (DFT) calculations starting from several models of AlON with different Al or O vacancy locations and different positions for the N atoms relative to the vacancy were carried out. The results show that the constant anion model [McCauley et al., J. Eur. Ceram. Soc. 29(2), 223 (2009)] with a random distribution of N atoms not adjacent to the Al vacancy has the lowest energy configuration. The lowest energy structure is in a reasonable agreement with experimental X-ray diffraction spectra. The optimized structure of a 55 atom unit cell was used to construct 220 and 440 atom models for simulation cells using DFT with a Gaussian basis set. Cubic elastic constant predictions were found to approach the experimentally determined AlON single crystal elastic constants as the model size increased from 55 to 440 atoms. The pressure dependence of the elastic constants found from simulated stress-strain relations were in overall agreement with experimental measurements of polycrystalline and single crystal AlON. Calculated IR intensity and Raman spectra are compared with available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Adaptive two-regime method: Application to front propagation

The Adaptive Two-Regime Method (ATRM) is developed for hybrid (multiscale) stochastic simulation of reaction-diffusion problems. It efficiently couples detailed Brownian dynamics simulations with coarser lattice-based models. The ATRM is a generalization of the previously developed Two-Regime Method [Flegg et al., J. R. Soc., Interface 9, 859 (2012)] to multiscale problems which require a dynamic selection of regions where detailed Brownian dynamics simulation is used. Typical applications include a front propagation or spatio-temporal oscillations. In this paper, the ATRM is used for an in-depth study of front propagation in a stochastic reaction-diffusion system which has its mean-field model given in terms of the Fisher equation [R. Fisher, Ann. Eugen. 7, 355 (1937)]. It exhibits a travelling reaction front which is sensitive to stochastic fluctuations at the leading edge of the wavefront. Previous studies into stochastic effects on the Fisher wave propagation speed have focused on lattice-based models, but there has been limited progress using off-lattice (Brownian dynamics) models, which suffer due to their high computational cost, particularly at the high molecular numbers that are necessary to approach the Fisher mean-field model. By modelling only the wavefront itself with the off-lattice model, it is shown that the ATRM leads to the same Fisher wave results as purely off-lattice models, but at a fraction of the computational cost. The error analysis of the ATRM is also presented for a morphogen gradient model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗