Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “N”

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 883 records · Page 49

Development of (NO)Fe(N 2 S 2 ) as a Metallodithiolate Spin Probe Ligand: A Case Study Approach

The ubiquity of sulfur–metal connections in nature inspires the design of bi- and multimetallic systems in synthetic inorganic chemistry. Common motifs for biocatalysts developed in evolutionary biology include the placement of metals in close proximity with flexible sulfur bridges as well as the presence of π-acidic/delocalizing ligands. This Account will delve into the development of a (NO)Fe(N 2 S 2 ) metallodithiolate ligand that harnesses these principles. The Fe(NO) unit is the centroid of a N 2 S 2 donor field, which as a whole is capable of serving as a redox-active, bidentate S-donor ligand. Its paramagnetism as well as the ν(NO) vibrational monitor can be exploited in the development of new classes of heterobimetallic complexes. We offer four examples in which the unpaired electron on the {Fe(NO)} 7 unit is spin-paired with adjacent paramagnets in proximal and distal positions. First, the exceptional stability of the (NO)Fe(N 2 S 2 )-Fe(NO) 2 platform, which permits its isolation and structural characterization at three distinct redox levels, is linked to the charge delocalization occurring on both the Fe(NO) and the Fe(NO) 2 supports. This accommodates the formation of a rare nonheme {Fe(NO)} 8 triplet state, with a linear configuration. A subsequent FeNi complex, featuring redox-active ligands on both metals (NO on iron and dithiolene on nickel), displayed unexpected physical properties. Our research showed good reversibility in two redox processes, allowing isolation in reduced and oxidized forms. Various spectroscopic and crystallographic analyses confirmed these states, and Mössbauer data supported the redox change at the iron site upon reduction. Oxidation of the complex produced a dimeric dication, revealing an intriguing magnetic behavior. The monomer appears as a spin-coupled diradical between {Fe(NO)} 7 and the nickel dithiolene monoradical, while dimerization couples the latter radical units via a Ni 2 S 2 rhomb. Magnetic data (SQUID) on the dimer dication found a singlet ground state with a thermally accessible triplet state that is responsible for magnetism. A theoretical model built on an H 4 chain explains this unexpected ferromagnetic low-energy triplet state arising from the antiferromagnetic coupling of a four-radical molecular conglomerate. For comparison, two (NO)Fe(N 2 S 2 ) were connected through diamagnetic group 10 cations producing diradical trimetallic complexes. Antiferromagnetic coupling is observed between {Fe(NO)} 7 units, with exchange coupling constants (J) of -3, -23, and -124 cm –1 for Ni II , Pd II , and Pt II , respectively. This trend is explained by the enhanced covalency and polarizability of sulfur-dense metallodithiolate ligands. A central paramagnetic trans-Cr(NO)(MeCN) receiver unit core results in a cissoid structural topology, influenced by the stereoactivity of the lone pair(s) on the sulfur donors. This {Cr(NO)} 5 radical bridge, unlike all previous cases, finds the coupling between the distal Fe(NO) radicals to be ferromagnetic (J = 24 cm –1 ). The stability and predictability of this S = 1/2 moiety and the steric/electronic properties of the bridging thiolate sulfurs suggest it to be a likely candidate for the development of novel molecular (magnetic) compounds and possibly materials. The role of synthetic inorganic chemistry in designing synthons that permit connections of the (NO)Fe(N 2 S 2 ) metalloligand is highlighted as well as the properties of the heterobi- and polymetallic complexes derived therefrom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapid Automated Annotation and Analysis of N-Glycan Mass Spectrometry Imaging Data Sets Using NGlycDB in METASPACE

Imaging N-glycan spatial distribution in tissues using mass spectrometry imaging (MSI) is emerging as a promising tool in biological and clinical applications. However, there is currently no high throughput tool for visualization and molecular annotation of N-glycans in MSI data, which significantly slows down data processing and hampers the applicability of this approach. In this work, we present how METASPACE, an open-source cloud engine for molecular annotation of MSI data, can be used to automatically annotate, visualize, analyze, and interpret high-resolution mass spectrometry-based spatial N-glycomics data. METASPACE is an emerging tool in spatial metabolomics, but the lack of compatible glycan databases has limited its application for comprehensive N-glycan annotations from MSI datasets. We created NGlycDB, a public database of N-glycans, by adapting available glycan databases. We demonstrate the applicability of NGlycDB in METASPACE by analyzing MALDI-MSI data from formalin-fixed paraffin-embedded (FFPE) human kidney and mouse lung tissue sections. We added NGlycDB to METASPACE for public use, thus facilitating applications of MSI in glycobiology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Situ Neutron Scattering Study of the Structure Dynamics of the Ru/Ca 2 N:e – Catalyst in Ammonia Synthesis

NH 3 synthesis is one of the most critical industrial processes. Compared to commercial iron catalysts, Ru catalysts show high intrinsic activity in this reaction but suffer from hydrogen poisoning. By loading Ru onto supports such as electrides and hydrides, the hydrogen poisoning problem can be significantly alleviated. However, relevant studies on the structural dynamics of the Ru/electride catalysts under reaction conditions are very scarce. Taking advantage of the high sensitivity to hydrogen species, it is possible to obtain insights into the structural changes during the reaction using in situ neutron techniques. In this study, we have investigated the structural evolution of the Ru/Ca 2 N:e – catalyst during the ammonia synthesis reaction by in situ neutron scattering (inelastic neutron scattering, INS) technique. In situ INS experiments suggest that Ca 2 N:e – is likely converted to the Ca 2 NH phase during the reaction. Unlike the previously known structure where H and N atoms are intermixed, the formed Ca 2 NH exhibits a segregated structure where the H and N atoms are located in different layers separated by the Ca layer. In conclusion, density functional theory calculations of the reaction energetics reveal that there are minor changes in the barriers and thermodynamics of the first N hydrogenation step between the two phases (Ca 2 NH phase with segregated H/N layers and intermixed Ca 2 NH phase), suggesting the impact of the phases on the reaction kinetics to be relatively minimal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗

Magnetic and EPR Spectroscopic Studies of Thiolate Bridged Divalent Ni, Pd, and Pt Ions Capped with VO(N 2 S 2 ) Metalloligands

Reactions of the metallodithiolate complex VO(bme-dach) (hereafter abbreviated as V, where bme-dach = N,N ′-bis(2-mercaptoethyl)-1,4-diazacycloheptane) with [Pd II (CH 3 CN) 4 ](BF 4 ) 2 and [Pt II (CH 3 CN) 4 ](BF 4 ) 2 yield the V–M–V trimetallic compounds [VPdV](BF 4 ) 2 (2) and [VPtV](BF 4 ) 2 (3). Reaction of a similar metalloligand, VO(bme-daco) (hereafter abbreviated as V′ where bme-daco = N,N′-bis(2-mercaptoethyl)-1,5-diazacyclooctane) with [Ni II (CH 3 CN) 6 ](BF 4 ) 2 afforded the related salt [V′NiV′](BF 4 ) 2 (1). X-ray structural analyses revealed that cations in 1, 2, and 3 adopt a stairstep C 2h structure consisting of two terminal VO(N 2 S 2 ) moieties bridged via thiolate sulfur to the group 10 metal ions. Weak ferromagnetic superexchange coupling (J = 0.282 cm –1 for 1, 0.954 cm –1 for 2, and 1.372 cm –1 for 3) was observed between the two S = 1/2 V IV centers separated by distances in the range of 5.9–6.3 Å, with J values varying following the order Ni < Pd < Pt. Frozen-solution EPR spectra measured on the more soluble [BArF 24 ] − (BArF 24 – = tetrakis((3,5-trifluoromethyl)phenyl)borate) analogues revealed that the [VPtV] 2+ cation exhibits a 15-line hyperfine splitting of 225 MHz at g = 4 in parallel mode, confirming exchange coupling between the two 51 V, I = 7/2 nuclei. Density-functional theory (DFT) calculations indicate an S = 1 ground state for 1–3. These results demonstrate that the choice of paramagnetic metallodithiolate ligand and diamagnetic bridge in such trimetallic species influences the sign and magnitude of magnetic interactions.

anions↗

One-Step Synthesis and Operando Electrochemical Impedance Spectroscopic Characterization of Heterostructured MoP–Mo 2 N Electrocatalysts for Stable Hydrogen Evolution Reaction

Here this study presents a novel synthesis of self-standing MoP and Mo 2 N heterostructured electrocatalysts with enhanced stability and catalytic performance. Facilitated by the controlled phase and interfacial microstructure, the seamless structures of these catalysts minimize internal resistivity and prevent local corrosion, contributing to increased stability. The chemical synthesis proceeds with etching step to activate the surface, followed by phosphor-nitriding in a chemical vapor deposition chamber to produce MoP-Mo 2 N@Mo heterostructured electrocatalysts. X-ray diffraction analyses confirm the presence of MoP, Mo 2 N, and Mo phases in the electrocatalyst. Morphology studies using scanning electron microscopy characterizes the hierarchical growth of structures, indicating successful formation of the heterostructure. X-ray photoelectron spectroscopy (XPS) analyses of the as-synthesized and post-catalytic activity samples reveal the chemical shift in terms of binding energy (BE) of Mo 3d XPS peak, especially after catalytic activity. The XPS BE shifts attributed to changes in oxidation state, electron transfer, and surface reconstruction during catalysis. Electrochemical evaluation of the catalysts demonstrates the superior performance of the MoP-Mo 2 N@Mo heterostructured catalyst in hydrogen evolution reactions (HER), with lower overpotentials and enhanced Tafel slopes. Stability tests reveal changes in double layer capacitance over time, suggesting surface reconstruction and increased active surface area during catalysis. Operando electrochemical impedance spectroscopy (EIS) further elucidates the dynamic changes in resistance and charge transfer during HER. Overall, comprehensive understanding of the synthesis, characterization, and electrochemical behavior of the developed MoP-Mo 2 N@Mo heterostructured electrocatalyst as presented in this paper highlights their potential utilization in sustainable energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Salts as Additives: A Route to Improve Performance and Stability of n-Type Organic Electrochemical Transistors

Organic electrochemical transistors (OECTs) are becoming increasingly ubiquitous in various applications at the interface with biological systems. However, their widespread use is hampered by the scarcity of electron-conducting (n-type) backbones and the poor performance and stability of the existing n-OECTs. Here, we introduce organic salts as a solution additive to improve the transduction capability, shelf life, and operational stability of n-OECTs. We demonstrate that the salt-cast devices present a 10-fold increase in transconductance and achieve at least one year-long stability, while the pristine devices degrade within four months of storage. The salt-added films show improved backbone planarity and greater charge delocalization, leading to higher electronic charge carrier mobility. These films show a distinctly porous morphology where the interconnectivity is affected by the salt type, responsible for OECT speed. The salt-based films display limited changes in morphology and show lower water uptake upon electrochemical doping, a possible reason for the improved device cycling stability. Our work provides a new and easy route to improve n-type OECT performance and stability, which can be adapted for other electrochemical devices with n-type films operating at the aqueous electrolyte interface.

36 MATERIALS SCIENCE↗

Synthetically Reversible, Proton-Mediated Nitrite N–O Bond Cleavage at a Dicopper Site

A monocationic dicopper(I,I) nitrite complex [Cu 2 (μ-κ 1 :κ 1 -O 2 N)DPFN][NTf 2 ] (2) (DPFN = 2,7-bis(fluoro-di(2-pyridyl)methyl)-1,8-naphthyridine, NTf 2 – = N(SO 2 CF 3 ) 2 – ), was synthesized by treatment of a dicopper acetonitrile complex, [Cu 2 (μ-MeCN)DPFN][NTf 2 ] 2 (1), with tetrabutylammonium nitrite ([nBu 4 N][NO 2 ]). DFT calculations indicate that 2 is one of three linkage isomers that are close in energy and presumably accessible in solution. Reaction of the μ-κ 1 :κ 1 -O 2 N complex with p-TolSH produces nitrous acid (HONO) and the corresponding dicopper thiolate species via an acid–base exchange reaction. Notably, treatment of 2 with HNTf 2 results in N–O bond cleavage in the putative, HONO-ligated complex to form the more thermodynamically favorable nitrosyl-bridged dicopper complex [Cu 2 (μ-NO)(μ-OH)DPFN][NTf 2 ] 2 (4). This scission can be reversed via deprotonation of the hydroxy ligand with KO t Bu. X-ray diffraction studies confirmed the solid-state molecular structures of 2 and 4. DFT calculations were used to construct a reaction coordinate diagram detailing formation of the μ-NO complex and to describe its electronic structure. The nitrosyl ligand in 4 is chemically labile, as demonstrated by its ready displacement in reactions with CO or NO 2 – .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unravelling the Stability Stressors of Atomically Dispersed Fe–N–C Oxygen Reduction Catalysts

Enhancing the catalytic stability of Fe–N–C catalysts for cathodic oxygen reduction in proton-exchange membrane fuel cells (PEMFCs) necessitates an in-depth understanding of their degradation mechanisms. Here, this study identifies key stressors affecting the stability of Fe–N–C catalysts, specifically acidic environment, oxygen (O 2 ), and reactive oxygen species (ROS). Through ex situ/operando experiments, we show that the oxidation of local carbon by acidic environment + O 2 + ROS, along with the demetalation of catalytic FeN x C y sites by O 2 or O 2 + ROS, is the primary factor responsible for the initial fast degradation of Fe–N–C catalysts. The demetalation of FeN x C y sites, influenced by O 2 , in particular by O 2 + ROS, leads to the subsequent gradual degradation of Fe–N–C. Notably, FeN 4 C 12 -type active sites are more susceptible to demetalation than FeN 4 C 10 -type sites in O 2 or O 2 + ROS. Our findings indicate that, besides constructing more stable FeN x C y sites, preventing local carbon oxidation and scavenging of ROS are all critical for maintaining the stability of Fe–N–C catalysts.

Xie, Xiaohong [Pacific Northwest National Laborato↗

Chemical manipulation of hydrogen induced high p-type and n-type conductivity in Ga 2 O 3

Advancement of optoelectronic and high-power devices is tied to the development of wide band gap materials with excellent transport properties. However, bipolar doping (n-type and p-type doping) and realizing high carrier density while maintaining good mobility have been big challenges in wide band gap materials. Here P-type and n-type conductivity was introduced in β-Ga 2 O 3 , an ultra-wide band gap oxide, by controlling hydrogen incorporation in the lattice without further doping. Hydrogen induced a 9-order of magnitude increase of n-type conductivity with donor ionization energy of 20 meV and resistivity of 10 -4 Ω.cm. The conductivity was switched to p-type with acceptor ionization energy of 42 meV by altering hydrogen incorporation in the lattice. Density functional theory calculations were used to examine hydrogen location in the Ga 2 O 3 lattice and identified a new donor type as the source of this remarkable n-type conductivity. Positron annihilation spectroscopy measurements confirm this finding and the interpretation of the experimental results. This work illustrates a new approach that allows a tunable and reversible way of modifying the conductivity of semiconductors and it is expected to have profound implications on semiconductor field. At the same time, it demonstrates for the first time p-type and remarkable n-type conductivity in Ga 2 O 3 which should usher in the development of Ga 2 O 3 devices and advance optoelectronics and high-power devices.

36 MATERIALS SCIENCE↗

Effect of metallosupramolecular polymer concentration on the synthesis of poly[n]catenanes

Poly[n]catenanes are a class of polymers that are composed entirely of interlocked rings. One synthetic route to these polymers involves the formation of a metallosupramolecular polymer (MSP) that consists of alternating units of macrocyclic and linear thread components. Ring closure of the thread components has been shown to yield a mixture of cyclic, linear, and branched poly[n]catenanes. Reported herein are investigations into this synthetic methodology, with a focus on a more detailed understanding of the crude product distribution and investigations into how the concentration of the MSP during the ring closing reaction impacts the resulting poly[n]catenanes. In addition to a better understanding of the molecular products obtained in these reactions, the results show that the concentration of the reaction can be used to tune the size and type of poly[n]catenanes accessed. At low concentrations the interlocked product distribution is limited to primarily oligomeric and small cyclic catenanes (number average molar mass $\overline{Mn}$ = 6.2 kg mol -1 , number of rings = 4). However, the same reaction at increased concentration can yield branched poly[n]catenanes with an $\overline{Mn}$ ca. 21 kg mol -1 and evidence of structures with as many as 640 interlocked rings (1000 kg mol -1 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis and symmetry of perovskite oxynitride CaW(O,N) 3

Perovskite oxynitrides, in addition to being promising electrocatalysts and photoabsorbers, present an interesting case study in crystal symmetry. Full or partial ordering of the O and N anions affects global symmetry and influences material performance and functionality; however, anion ordering is challenging to detect experimentally. In this work, we synthesize a novel perovskite oxynitride CaW(O,N) 3 and characterize its crystal structure using both X-ray and neutron diffraction. Through co-refinement of the diffraction patterns with a range of literature and theory-derived model structures, we demonstrate that CaW(O,N) 3 adopts an orthorhombic Pnma average structure and exhibits octahedral distortion with evidence for preferred anion site occupancy. However, through comparison with a large, low-symmetry unit cell, we identify the presence of disorder that is not fully accounted for by the high-symmetry model. We compare CaW(O,N) 3 with SrW(O,N) 3 to demonstrate the broader presence of such disorder and identify contrasting features in the electronic structures. This work signifies an updated perspective on the inherent crystal symmetry present in perovskite oxynitrides.

36 MATERIALS SCIENCE↗

Unveiling the covalency of versatile Pu(iii)-N bonds in a unique plutonium(iii) complex

A trivalent plutonium–pyrazinyl–tetrazolate complex Na 2 [Pu(Hdtp)(dtp) 2 (H 2 O) 4 ]·9H 2 O (Pu_dtp, H 2 dtp = 2,3-di-1H-tetrazol-5-ylpyrazine) was synthesized through metathesis reaction of plutonium bromide and Na 2 (dtp)·2H 2 O in water. This structure is particularly notable among complexes formed by trivalent f-elements and the dtp 2− ligand in aqueous media. In contrast to other trivalent f-elements, including all Ln 3+ (with the exception of Pm 3+ ) and Cm 3+ , preferentially coordinated with eight water molecules rather than the nitrogen donors of the dtp 2− ligand, Pu 3+ exhibits a distinct affinity for nitrogen coordination. This observation provides strong evidence that the 5f electrons in Pu 3+ are more delocalized than other studied trivalent f-elements. In Pu_dtp, three distinct Pu(III)–N bonds are present: one Pu(III)–N5 from pyrazinyl, one Pu(III)–N4 from the least electronegative nitrogen in the tetrazolate, and three Pu(III)–N1/N2/N3 from the most electronegative nitrogens in the tetrazolate. Experimental Pu(III)–N bond lengths, Wiberg bond indices (WBI), natural localized molecular orbitals (NLMO), quantum theory of atoms in molecules (QTAIM), and energy decomposition analysis (EDA), reveal a covalency trend: Pu(III)–N from the most electronegative nitrogen in tetrazolate > Pu(III)–N from the least electronegative nitrogen in tetrazolate > Pu(III)–N from pyrazinyl. This trend arises from the increased negative charge on the most electronegative nitrogen atoms in the tetrazolate ring, enhancing electrostatic Pu–N1/N2/N3 interactions. These stronger electrostatic interactions lead to shorter bond lengths, thereby enhancing orbital overlap and greater covalency, compared to the less electronegative nitrogen in tetrazolate (Pu–N4) and the neutral pyrazinyl nitrogen (Pu–N5).

Bai, Zhuanling [Colorado School of Mines, Golden, ↗

Galactic ArchaeoLogIcaL ExcavatiOns (GALILEO): I. An updated census of APOGEE N-rich giants across the Milky Way

We use the 17th data release of the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) to provide a homogenous census of N-rich red giant stars across the Milky Way (MW). We report a total of 149 newly identified N-rich field giants toward the bulge, metal-poor disk, and halo of our Galaxy. They exhibit significant enrichment in their nitrogen abundance ratios ([N/Fe] ≳ +0.5), along with simultaneous depletions in their [C/Fe] abundance ratios ([C/Fe] < +0.15), and they cover a wide range of metallicities (-1.8 < [Fe/H] < -0.7). The final sample of candidate N-rich red giant stars with globular-cluster-like (GC-like) abundance patterns from the APOGEE survey includes a grand total of ~412 unique objects. These strongly N-enhanced stars are speculated to have been stripped from GCs based on their chemical similarities with these systems. Even though we have not found any strong evidence for binary companions or signatures of pulsating variability yet, we cannot rule out the possibility that some of these objects were members of binary systems in the past and/or are currently part of a variable system. In particular, the fact that we identify such stars among the field stars in our Galaxy provides strong evidence that the nucleosynthetic process(es) producing the anomalous [N/Fe] abundance ratios occurs over a wide range of metallicities. This may provide evidence either for or against the uniqueness of the progenitor stars to GCs and/or the existence of chemical anomalies associated with likely tidally shredded clusters in massive dwarf galaxies such as “Kraken/Koala”, Gaia-Enceladus-Sausage, among others, before or during their accretion by the MW. In conclusion, a dynamical analysis reveals that the newly identified N-rich stars exhibit a wide range of dynamical characteristics throughout the MW, indicating that they were produced in a variety of Galactic environments.

59 BASIC BIOLOGICAL SCIENCES↗

Probing electronic dead layers in homoepitaxial n-SrTiO 3 (001) films

We combine state-of-the-art oxide epitaxial growth by hybrid molecular beam epitaxy with transport, x-ray photoemission, and surface diffraction, along with classical and first-principles quantum mechanical modeling to investigate the nuances of insulating layer formation in otherwise high-mobility homoepitaxial n-SrTiO 3 (001) films. Our analysis points to charge immobilization at the buried n-SrTiO 3 /undoped SrTiO 3 (001) interface as well as within the surface contamination layer resulting from air exposure as the drivers of electronic dead-layer formation. As Fermi level equilibration occurs at the surface and the buried interface, charge trapping reduces the sheet carrier density (n 2D ) and renders the n-STO film insulating if n 2D falls below the critical value for the metal-to-insulator transition.

36 MATERIALS SCIENCE↗

N -representability violations in truncated equation-of-motion coupled-cluster methods

One-electron reduced density matrices (1RDMs) from equation-of-motion (EOM) coupled-cluster with single and double excitations (CCSD) calculations are analyzed to assess their N-representability (i.e., whether they are derivable from a physical N-electron state). We identify EOM-CCSD stationary states whose 1RDMs violate either ensemble-state N-representability conditions or pure-state conditions known as generalized Pauli constraints. As such, these 1RDMs do not correspond to any physical N-electron state. Unphysical states are also encountered in the course of time-dependent EOM-CC simulations; when an external field drives transitions between a pair of stationary states with pure-state N-representable 1RDMs, the 1RDM of the time-dependent state can violate ensemble-state conditions. Furthermore, these observations point to potential challenges in interpreting the results of time-dependent EOM-CCSD simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pitfalls in the n -mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Chemistry↗

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗