Testing the r 2 SCAN Density Functional for the Thermodynamic Stability of Solids with and without a van der Waals Correction
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “energy density functional”
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.
Nonintuitive Surface Self-Assembly of Functionalized Molecules on Ag(111)
The fabrication of nanomaterials involves selfordering processes of functional molecules on inorganic surfaces. To obtain specific molecular arrangements, a common strategy is to equip molecules with functional groups. However, focusing on the functional groups alone does not provide a comprehensive picture. Especially at interfaces, processes that govern self-ordering are complex and involve various physical and chemical effects, often leading to unexpected structures, as we showcase here on the example of a homologous series of quinones on Ag(111). Naively, one could expect that such quinones, which all bear the same functionalization, form similar motifs. In salient contrast, our joint theoretical and experimental study shows that profoundly different structures are formed. Using a machine-learning-based structure search algorithm, we find that this is due to a shift of the balance of three antagonizing driving forces: adsorbate-substrate interactions governing adsorption sites, adsorbate-adsorbate interactions favoring close packing, and steric hindrance inhibiting certain otherwise energetically beneficial molecular arrangements. The theoretical structures show excellent agreement with our experimental characterizations of the organic/ inorganic interfaces, both for the unit cell sizes and the orientations of the molecules within. The nonintuitive interplay of similarly important interaction mechanisms will continue to be a challenging aspect for the design of functional interfaces. With a detailed examination of all driving forces, we are, however, still able to devise a design principle for self-assembly of functionalized molecules.
structure prediction↗
Proton Relay for the Rate Enhancement of Electrochemical Hydrogen Reactions at Heterogeneous Interfaces
Proton transfer is critically important to many electrocatalytic reactions, and directed proton delivery could open new avenues for the design of electrocatalysts. However, although this approach has been successful in molecular electrocatalysis, proton transfer has not received the same attention in heterogeneous electrocatalyst design. Here, in this work, we report that a metal oxide proton relay can be built within heterogeneous electrocatalyst architectures and improves the kinetics of electrochemical hydrogen evolution and oxidation reactions. The volcano-type relationship between activity enhancement and pK a of amine additives confirms this improvement; we observe maximum rate enhancement when the pK a of a proton relay matches the pH of the electrolyte solution. Density-functional-theory-based reactivity studies reveal a decreased proton transfer energy barrier with a metal oxide proton relay. These findings demonstrate the possibility of controlling the proton delivery and enhancing the reaction kinetics by tuning the chemical properties and structures at heterogeneous interfaces.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Electronically-coupled redox centers in trimetallic cobalt complexes
Synthesis and isolation of molecular building blocks of metal–organic frameworks (MOFs) can provide unique opportunities for characterization that would otherwise be inaccessible due to the heterogeneous nature of MOFs. Herein, we report a series of trinuclear cobalt complexes incorporating dithiolene ligands, triphenylene-2,3,6,7,10,11-hexathiolate (THT) (1 3+ ), and benzene hexathiolate (BHT) (2 3+ ), with 1,1,1,-tris(diphenylphosphinomethyl)ethane (triphos) employed as the capping ligand. Single crystal X-ray analyses of 1 3+ and 2 3+ display three five-coordinate cobalt centers bound to the triphos and dithiolene ligands in a distorted square pyramidal geometry. Cyclic voltammetry studies of 1 3+ and 2 3+ reveal three redox features associated with the formation of mixed valence states due to the sequential reduction of the redox-active metal centers (Co III/II ). Using this electrochemical data, the comproportionality values were determined for 1 and 2 (log K c = 1.4 and 1.5 for 1, and 4.7 and 5.8 for 2), suggesting strong resonance-stabilized coupling of the metal centers, with stronger electronic coupling observed for complex 2 compared to that for complex 1. Cyclic voltammetry studies were also performed in solvents of varying polarity, whereupon the difference in the standard potentials (ΔE 1/2 ) for 1 and 2 was found to shift as a function of the polarity of the solvent, indicating a negative correlation between the dielectric constant of the electrochemical medium and the stability of the mixed valence species. Spectroelectrochemical studies of in situ generated multi-valent (MV) states of complexes 1 and 2 display characteristic NIR intervalence charge transfer (IVCT) bands, and analysis of the IVCT transitions for complex 2 suggests a weakly coupled class II multi-valent species and relatively large electronic coupling factors (1700 cm –1 for the first multi-valent state of 2 2+ , and 1400 and 4000 cm –1 for the second multi-valent state of 2 + ). Here, density functional theory (DFT) calculations indicate a significant deviation in relative energies of the frontier orbitals of complexes 1 3+ , 2 3+ , and 3 + that contrasts those calculated for the analogous trinuclear cobalt dithiolene complexes employing pentamethylcyclopentadienyl (Cp*) as the capping ligand (Co 3 Cp* 3 THT and Co 3 Cp* 3 BHT, respectively), and may be a result of the cationic nature of complexes 1 3+ , 2 3+ , and 3 + .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
The striking influence of oxophilicity differences in heterometallic Mo–Mn oxide cluster reactions with water
Mixed-metal oxides have proven to be effective catalysts for the hydrogen evolution reaction, often outperforming either of the binary metal oxides. The reactivity of Mn x MoO y - (x = 1, 2; y = 3, 4) clusters toward H 2 O was investigated via time-of-flight mass spectrometry with clear evidence of cluster oxidation and corresponding H 2 production, specifically for Mn x MoO 3 - (x = 1, 2) clusters. Unlike previously studied MoxOy- clusters, which assumed a broad distribution of stoichiometries (typically x ≤ y ≤ 3x), both MnMoOy- and Mn 2 MoO y - preferentially formed y = 3 and 4 compositions in significant quantities under our source conditions. The electronic and molecular structures of the Mn x MoO y (x = 1, 2; y = 3, 4) anion and neutral clusters were probed with anion photoelectron spectroscopy and analyzed with supporting density functional theory calculations. Our studies suggest that both metal centers are involved in initial cluster–water complex formation, while Mo is the center that undergoes oxidation; hence, reactivity terminates when Mo is saturated in its highest oxidation state of +6. Across these four clusters, Mn remains relatively reduced and is stable in a high-spin electronic configuration. The preferential reactivity of water molecules toward the Mo center rather than Mn is rationalized by the much lower relative oxophilicity of Mn.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation
Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point
The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NiGa 2 O 4 interfacial layers in NiO/Ga 2 O 3 heterojunction diodes at high temperature
NiO/Ga 2 O 3 heterojunction diodes have attracted attention for high-power applications, but their high temperature performance and reliability remain underexplored. Here, we report the time evolution of the electrical properties in the widely studied p-NiO/n-Ga 2 O 3 heterojunction diodes and formation of NiGa 2 O 4 interfacial layers at high temperatures. Results of our thermal cycling experiment show an initial leakage current increase which stabilizes after sustained thermal load, due to reactions at the NiO–Ga 2 O 3 interface. High-resolution TEM microstructure analysis of the devices after thermal cycling indicates that the NiO–Ga 2 O 3 interface forms a ternary compound at high temperatures, and thermodynamic calculations suggest the formation of the spinel NiGa 2 O 4 layer between NiO and Ga 2 O 3 . First-principles defect calculations find that NiGa 2 O 4 shows low p-type intrinsic doping and hence can serve to limit electric field crowding at the interface. Vertical NiO/ Ga 2 O 3 diodes with intentionally grown 5 nm thin spinel-type NiGa 2 O 4 interfacial layers show an excellent device ON/OFF ratio of >10 10 (± 3 V), V ON of ~1.9 V, and increased breakdown voltage of ~1.2 kV for an initial unoptimized 300 lm diameter device. These p–n heterojunction diodes are promising for high-voltage, high temperature applications.
30 DIRECT ENERGY CONVERSION↗
Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics
In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.
Sharma, Vidushi [Princeton Plasma Physics Laborato↗
Tailoring electrocatalytic activity of in situ crafted perovskite oxide nanocrystals via size and dopant control
Significance The ability to scrutinize the correlation of dimension, composition, and dopant to electrocatalytic performance renders the development of highly active electrocatalysts. This work reports a general and robust strategy for crafting uniform perovskite oxide nanoparticles (i.e., BaTiO 3 and La- and Co-doped BaTiO 3 ) with controlled sizes and dopant compositions by employing amphiphilic star-like diblock copolymers as nanoreactors. Their size- and dopant-dependent oxygen reduction reaction (ORR) performances are unveiled. Particularly, La- and Co-doped BaTiO 3 nanoparticles exhibit markedly improved ORR activities over the pristine counterparts due to reduced free energy barrier of ORR and increased conductivity as substantiated by density functional theory. The nanoreactor strategy could be conveniently extended to yield various functional nanomaterials of interest with tailored sizes and compositions for electrocatalysis.
36 MATERIALS SCIENCE↗
Molecular surface programming of rectifying junctions between InAs colloidal quantum dot solids
Heavy-metal-free III–V colloidal quantum dots (CQDs) show promise in optoelectronics: Recent advancements in the synthesis of large-diameter indium arsenide (InAs) CQDs provide access to short-wave infrared (IR) wavelengths for three-dimensional ranging and imaging. In early studies, however, we were unable to achieve a rectifying photodiode using CQDs and molybdenum oxide/polymer hole transport layers, as the shallow valence bandedge (5.0 eV) was misaligned with the ionization potentials of the widely used transport layers. This occurred when increasing CQD diameter to decrease the bandgap below 1.1 eV. Here, we develop a rectifying junction among InAs CQD layers, where we use molecular surface modifiers to tune the energy levels of InAs CQDs electrostatically. Previously developed bifunctional dithiol ligands, established for II-VI and IV-VI CQDs, exhibit slow reaction kinetics with III-V surfaces, causing the exchange to fail. Here, we study carboxylate and thiolate binding groups, united with electron-donating free end groups, that shift upward the valence bandedge of InAs CQDs, producing valence band energies as shallow as 4.8 eV. Photophysical studies combined with density functional theory show that carboxylate-based passivants participate in strong bidentate bridging with both In and As on the CQD surface. The tuned CQD layer incorporated into a photodiode structure achieves improved performance with EQE (external quantum efficiency) of 35% (>1 μm) and dark current density < 400 nA cm -2 , a >25% increase in EQE and >90% reduced dark current density compared to the reference device. This work represents an advance over previous III-V CQD short-wavelength IR photodetectors (EQE < 5%, dark current > 10,000 nA cm -2 ).
36 MATERIALS SCIENCE↗
Nucleonic shells and nuclear masses
The binding energy of an isotope is a sensitive indicator of the underlying shell structure as it reflects the net energy content of a nucleus. Since magic nuclei are significantly lighter, or more bound, compared to their neighbors, the presence of nucleonic shell structure makes an imprint on nuclear masses. In this study, using a carefully designed binding-energy indicator, we catalog the appearance of spherical and deformed shell and subshell closures throughout the nuclear landscape. After presenting experimental evidence for shell and subshell closures as seen through the lens of nuclear masses, we study the ability of global nuclear mass models to predict local binding-energy variations related to shell effects.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Boron nitride on SiC(0001)
Here, in the field of van der Waals heterostructures, the twist angle between stacked two-dimensional layers has been identified to be of utmost importance for the properties of the heterostructures. In this context, we previously reported the growth of a single layer of unconventionally oriented epitaxial graphene that forms in a surfactant atmosphere [F. C. Bocquet et al., Phys. Rev. Lett. 125, 106102 (2020)]. The resulting G-R0° layer is aligned with the SiC lattice, and hence represents an important milestone towards high-quality twisted bilayer graphene, a frequently investigated model system in this field. Here, we focus on the surface structures obtained in the same surfactant atmosphere, but at lower preparation temperatures at which a boron nitride template layer forms on SiC(0001). In a comprehensive study based on complementary experimental and theoretical techniques, we find—in contrast to the literature—that this template layer is a hexagonal B x N y layer, but not high-quality hBN. It is aligned with the SiC lattice and gradually replaced by low-quality graphene in the 0° orientation of the B x N y template layer upon annealing.
36 MATERIALS SCIENCE↗
First-principles investigation of phase stability in layered Na x CrO 2
Layered oxide intercalation compounds continue to attract interest as electrode materials for Na-ion batteries. However, many of these materials undergo complex phase transitions during cycling that influence battery performance but are still not completely understood. We have conducted a first-principles study of layered Na x CrO 2 (0 ≤ x ≤ 1) to assess phase stability between various Na-vacancy ordered phases in the O3 and P3 host structures. We predict that many of the low-energy phases belong to families of Na orderings that follow specific patterns. At high x, we identify families of vacancy row orderings in O3, which may also couple to magnetic orderings of the Cr spins. We predict similar orderings at intermediate x in P3 that contain antiphase boundaries between regions of the x = 1/2 ordering. In both cases, the average spacing between rows/boundaries is set by the overall composition. At x = 0, we find a strong preference for charge disproportionation and migration of Cr to tetrahedral sites in the intercalation layers. Finally, we obtain generally good agreement with experimental observations and rationalize key discrepancies.
25 ENERGY STORAGE↗
Coupled Process/Device Modeling and Point Defect Engineering of Cu(In,Ga)Se 2 Solar Cells
Point defects directly impact solar cell device performance by limiting the carrier lifetime. In this work, density functional theory calculations are first used to determine the formation energy and diffusion energy barriers of dominant defects in Cu(In,Ga)Se 2 . Next, continuum reaction–diffusion models are developed to analyze the redistribution of defects during manufacturing processes. We estimate defect capture cross sections using a first-principles-based approach. Here, these cross sections are combined with our calculated defect profiles and trap energy levels to parameterize a Shockley–Read–Hall recombination model, which we implement into a device simulator to predict carrier lifetimes and device performance. In that way, a predictive technology computer aided design model is built to predict and optimize the performance of Cu(In,Ga)Se 2 solar cells.
14 SOLAR ENERGY↗
Van der Waals density functional study of hydrocarbon adsorption and separation in metal–organic frameworks without open metal sites
Metal–organic frameworks (MOFs) have received significant attention thanks to their promising features in the storage and separation of guest molecules. MOFs without open metal sites are emerging as they are often less susceptible to poisoning compared to those with open metal sites. However, a complete understanding of the binding and gas separation mechanisms in such materials is still missing. In this work, we perform a comparative study of two classes of vanadium-based MOFs without open metal sites: MFM-300-V (III) and MFM-300-V (IV) , as well as MIL-47-V (III) and MIL-47-V (IV) . We employ first-principles van der Waals density functional theory to find the optimal binding conformations and binding energies of a series of small hydrocarbons within the pores of the aforementioned MOFs. Our study provides insight into the host–guest interactions in such MOFs without open metal sites, especially the role played by the bridging hydroxyl group (μ 2 –OH). We conclude that the bridging –OH group acts as a pseudo open metal site in these MOFs.
36 MATERIALS SCIENCE↗
Elucidating Hydrogen Reaction-Induced Water Desorption from Oxide-Passivated Metal Surfaces for Plasma Applications
Elucidating the mechanisms responsible for sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for understanding pulsed power behavior. Ionization of desorbed impurities in the vacuum regions causes power or current loss; devising methods to limit such desorption during the short time scale of pulsed power is needed to improve corresponding applications. Previous molecular modeling studies have strongly suggested that, under high vacuum conditions, the amount of water impurity desorbing from oxide surfaces on metal electrodes is at a sub-monolayer level at room temperature, which appears insufficient to explain observed pulsed power energy losses at high current densities. In this work, we apply Density Functional Theory (DFT) techniques to show that hydrogen trapped inside iron metal can diffuse into hematite (α-Fe 2 O 3 ) on the metal surface, ultimately reacting with the oxide to form Fe(II) and H 2 O. The latter desorbs at elevated temperature and may explain the anomalous amount of desorbed impurity inferred from pulsed-power experiments. We also apply a suite of characterization techniques to demonstrate that when iron metal is heated to 650 °C, the dominant surface oxide component becomes α-Fe 2 O 3 . The oxide facets exposed are found to be a mixture of (0001), (10-10), and others, in agreement with the DFT models used.
74 ATOMIC AND MOLECULAR PHYSICS↗