Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cluster chemistry”

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 253 records · Page 14

Porting Fragmentation Methods to Graphical Processing Units Using an OpenMP Application Programming Interface: Offloading the Fock Build for Low Angular Momentum Functions

Here, a framework to offload four-index two-electron repulsion integrals to graphical processing units (GPUs) using OpenMP is discussed. The method has been applied to the Fock build for low angular momentum s and p functions in both the restricted Hartree–Fock (RHF) and in the effective fragment molecular orbital (EFMO) framework. Benchmark calculations for the GPU code for the pure RHF method show an increasing speedup relative to the existing OpenMP CPU code in GAMESS from 1.04 to 52× for clusters of 70–569 water molecules. The parallel efficiency on 24 NVIDIA V100 GPU boards also increases when increasing the system size: from 75 to 94% for water clusters that contain 303–1120 molecules. In the EFMO framework, the GPU Fock build shows a high linear scalability up to 4608 V100s with a parallel efficiency of 96% for calculations on a solvated mesoporous silica nanoparticle system with ~67,000 basis functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of Cation-Induced Conformational Changes in Nanometer-Sized Uranyl Peroxide Clusters

Conformational changes of the pyrophosphate (Pp)-functionalized uranyl peroxide nanocluster [(UO 2 ) 24 (O 2 ) 24 (P 2 O 7 ) 12 ] 48– ({U 24 Pp 12 }), dissolved as a Li/Na salt, can be induced by the titration of alkali cations into solution. The most symmetric conformer of the molecule has idealized octahedral ( O h ) molecular symmetry. One-dimensional 31 P NMR experiments provide direct evidence that both K + and Rb + ions trigger an O h -to- D 4 h conformational change within {U 24 Pp 12 }. Variable-temperature 31 P NMR experiments conducted on partially titrated {U 24 Pp 12 } systems show an effect on the rates; increased activation enthalpy and entropy for the D 4 h -to- O h transition is observed in the presence of Rb + compared to K + . Two-dimensional, exchange spectroscopy 31 P NMR revealed that magnetization transfer links chemically unique Pp bridges that are present in the D 4 h conformation and that this magnetization transfer occurs via a conformational rearrangement mechanism as the bridges interconvert between two symmetries. The interconversion is triggered by the departure and reentry of K (or Rb) cations out of and into the cavity of the cluster. This rearrangement allows Pp bridges to interconvert without the need to break bonds. Cs ions exhibit unique interactions with {U 24 Pp 12 } clusters and cause only minor changes in the solution 31 P NMR signatures, suggesting that O h symmetry is conserved. Single-crystal X-ray diffraction measurements reveal that the mixed Li/Na/Cs salt adopts D 2 h molecular symmetry, implying that while solvated, this cluster is in equilibrium with a more symmetric form. Thus, these results highlight the unusually flexible nature of the actinide-based {U 24 Pp 12 } and its sensitivity to countercations in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion–ion and ion–solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K + ion. Diffusion-based Nernst–Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

25 ENERGY STORAGE↗

Gas-phase stability of large lanthanide:diglycolamide clusters evaluated using collision-induced dissociation

Efficient utilization of long-term deep geologic storage repositories is critical for the large-scale deployment of nuclear energy. As the minor actinides (americium and curium) are the largest contributor to the decay heat of used nuclear fuel after a few hundred years, their removal from used fuel prior to disposal can significantly increase the amount of waste that can be stored in a given volume. The development of ligands that can effectively separate the minor actinides from other components of used nuclear fuel are crucial to efficient utilization of geological storage repositories. N,N,N’,N’-tetraoctyl diglycolamide (TODGA) is a promising extractant for the separation of lanthanides and the minor actinides from other components of used nuclear fuel. While the use of TODGA in f-element separations has been investigated in process-based formulations, gas-phase metal ion cluster experiments enable the study of covalent interactions in reprocessing systems absent from solvent effects. Exploring fundamental differences in lanthanide-ligand covalent interactions can impact the development and implementation of actinide-lanthanide separation systems in nuclear fuels reprocessing. In this study, lanthanide-TODGA clusters were synthesized in the gas-phase and identified using mass spectrometry fragmentation experiments. Large europium and samarium clusters were identified that contained up to 8 and 10 bound TODGA ligands, respectively; this was surprising due to the size and multidentate binding that is normal for TODGA. In addition, while smaller clusters showed evidence of sequential ligand fragmentation, larger clusters displayed the loss of neutral TODGA with applied collision voltage. Interestingly, the voltage required for this removal decreased as more TODGA ligands were bound to the metal, suggesting that the metal coordination sphere was becoming more saturated and TODGA ligands were more weakly bound as the clusters got larger.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence for Low-Valent Electronic Configurations in Iron–Sulfur Clusters

Although biological iron-sulfur (Fe–S) clusters perform some of the most difficult redox reactions in Nature, they are thought to be composed exclusively of Fe 2+ and Fe 3+ ions, as well as mixed-valent pairs with average oxidation states of Fe 2.5+ . We herein show that Fe–S clusters formally composed of these valences can access a wider range of electronic configurations—in particular, those featuring low-valent Fe 1+ centers. We demonstrate that CO binding to a synthetic [Fe 4 S 4 ] o cluster supported by N-heterocyclic carbene ligands induces generation of Fe 1+ centers via intracluster electron transfer, wherein a neighboring pair of Fe 2+ sites reduces the CO-bound site to a low-valent Fe 1+ state. Similarly, CO binding to an [Fe 4 S 4 ] + cluster induces electron delocalization with a neighboring Fe site to form a mixed-valent Fe 1.5+ Fe 2.5+ pair in which the CO-bound site adopts partial low-valent character. Furthermore, these low-valent configurations engender remarkable C–O bond activation without having to traverse highly negative and physiologically inaccessible [Fe 4 S 4 ] o /[Fe 4 S 4 ] – redox couples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward the “platinum standard” of quantum chemistry on quantum computers: Perturbative quadruple corrections in unitary coupled cluster theory

We propose a non-iterative, post-hoc correction to the unitary coupled cluster theory with the single, double, and triple excitations (UCCSDT) Ansatz, which considers the leading-order effects of neglected quadruple excitations. We present two ways to derive this correction, henceforth referred to as [Q-6], which leads to an improvement in the correlation energy shown to be truncated to sixth-order in many-body perturbation theory. Furthermore, a comparison between the UCC-based [Q-6] correction proposed in this work and analogous, “platinum standard” quadruple corrections proposed in conventional coupled cluster theory recognizes that [Q-6] is distinct from prior corrections since it is constructed entirely from internally connected components. Although trotterized (t) and full operator variants of UCCSDT exhibit errors in scans of small molecule potential energy surfaces that routinely exceed 1.6 mH, we find that t/UCCSDT[Q-6] is, nevertheless, able to achieve chemical accuracy as measured by the mean unsigned error.

Correlation energy↗

Photogenerated Charge Separation in a CdSe Nanocluster Encapsulated in a Metal–Organic Framework for Improved Photocatalysis

Nanoclusters encapsulated in a metal–organic framework (MOF) are the subject of intense interest because of their high promise for photocatalysis. We show computationally that a CdSe nanocluster incorporated into NU-1000 (Cd 6 Se 6 @NU-1000) allows electron transfer from a visible-light photo-excited organic linker to the lowest unoccupied crystal orbital of the inorganic cluster, which can lead to charge separation. Substituents on the linker can be used to tune the band position and the band gap of the whole composite (Cd 6 Se 6 @NU-1000), and this makes it possible to perform the water splitting reaction under visible light. Here, the achievement of charge separation within the MOF leads to longer life of the excited states, and as a result, it can enhance the photocatalytic activity of the system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the Role of Neutral 4-Amino-1,2,4-triazole in the Formation of Hexanuclear f-Element Hydrolysis Products

Our recent observations of an unexpected Ce(III) hydrolysis product from the reaction of 4-amino-1,2,4-triazole (4-NH 2 -1,2,4-Triaz) with CeCl 3 ·7H 2 O, [Ce 6 (μ 3 -O) 4 (μ 3 -OH) 2 (μ 3 -Cl) 2 (Cl) 6 (μ 2 -4-NH 2 -1,2,4-Triaz) 12 ]·7H 2 O, the first high-nuclearity lanthanide complex where all Ln atoms are connected pairwise through 12 N-donor ligands or 12 neutral bridging ligands of any type, prompted us to explore the utility of this ligand in trapping additional f-element examples. Reactions of LnCl 3 ·6H 2 O (Ln = Nd, Eu, Ho) with a large excess of 4-NH 2 -1,2,4-Triaz (20 equiv) and with the addition of small amounts of water to help solubilize the metal salts led to the isolation of the unique hydrolysis products [Nd 6 (μ 3 -OH) 8 Cl 6 (μ 2 -4-NH 2 -1,2,4-Triaz) 12 ][Cl 4 ]·2H 2 O, [Eu 6 (μ 6 -Cl) 0.23 (μ 3 -O 0.77 ) 4 (μ 3 -O) 2.6 (μ 3 -Cl) 0.4 Cl 6 (μ 2 -4-NH 2 -1,2,4-Triaz) 12 ], and [Ho 6 (μ 6 -Cl) 0.21 (μ 3 -O 0.79 ) 4 (μ 3 -OH) 2 Cl 6 (μ 2 -4-NH 2 -1,2,4-Triaz) 12 ][Cl] 3.4 . Here, we also report a Ce(III) analogue prepared in glassware contaminated with Pb(OAc) 2 , namely, [Ce 6 (μ 3 -OH) 8 (BrPbBr 5 )(μ 2 -4-NH 2 -1,2,4-Triaz) 11.5 (OH 2 ) 6 ][Pb 0.84 Br 4.2 ][Br] 3.8 ·2(4-NH 2 -1,2,4-Triaz)·3.6H 2 O. The Nd(III) complex is the structurally most ordered with a clear [Nd 6 (μ 3 -OH) 8 ] cluster core, while the Eu(III) and Ho(III) compounds contain partial occupancy of a μ 6 position and thus result in an incomplete Ln 6 O 9 cluster core formation. The crystallographic results suggest that the 4-NH 2 -1,2,4-Triaz ligand brings Ln(III) ions together, followed by the formation of an Ln 6 O 8 or Ln 6 O 9 core with whatever remaining anions or ligands can be incorporated. Given the complexity of the hydrolysis products of nuclear waste, we expect to continue to find a myriad of closely related complex structures of these types for the f-elements.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Photoinduced Electron Transfer Informs on Pathway Coupling in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is an enzymatic mechanism that generates extremely high-energy electrons to drive unfavorable chemical reactions. It is utilized by the NADHdependent ferredoxin:NADP + -oxidoreductase (Nfn) enzyme in hyperthermophile Pyrococcus furiosus to bifurcate electrons from NADPH into the coupled low-potential (endergonic) and highpotential (exergonic) pathways. This process enables P. furiosus to live in harsh and uninviting environments. Despite its biological importance, the mechanisms used by Nfn to facilitate exceptional directional control over short-lived, high-energy electrons and to prevent undesired transfer, particularly along the low-potential pathway, are still not well understood. To elucidate how the protein environment contributes to electronic control in the lowpotential pathway, new techniques must be utilized to probe these unstable intermediates. In this study, we have adapted lowtemperature photoexcitation combined with electron paramagnetic resonance (EPR) to accumulate the short-lived intermediate and place it in the context of the other cofactors involved in the low-potential pathway of Nfn. We observed coincident growth of both the radical intermediate and its nearby [4Fe-4S] cluster over 4.5 h of illumination with NADPH at cryogenic temperatures. The photogenerated paramagnetic species were stable in LN 2 storage indefinitely and recombined when warmed to higher temperatures. The results provide insights into the electron transfer steps and cofactor interactions along the low potential pathway, facilitating a more robust mechanistic understanding of the high-energy events of electron bifurcation. Furthermore, through comparison of cryogenic and room temperature experiments, a potential gating step involving the movement of key residues important for the reversibility of electron flow along this pathway is suggested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Venturing Past Uranium: Synthesis of a Np(IV) Polyoxomolybdate–Alkoxide Sandwich Complex

The synthesis of a Np(IV) polyoxomolybdate– alkoxide sandwich complex, (TBA) 2 [Np{Mo 5 O 13 (OMe) 4 NO} 2 ] (TBA = tetrabutylammonium), is reported. This compound represents a rare example of a neptunium polyoxometalate cluster isolated outside of water, allowing for characterization of its electrochemical properties in nonaqueous solvents. Complexation of An(IV) cations fine-tunes the redox properties of the cluster, with the observed four reversible reductive events varying slightly both in potential and peak separation depending on the actinide present. The new Np(IV) complex also shows an irreversible event assigned to oxidation of Np(IV) to Np(V). New methodology for facile 17 O enrichment of (TBA) 2 [Mo 5 O 13 (OMe) 4 NO][Na(MeOH)] is presented, which provides a simple pathway to 17 O enriched analogues of the sandwich complexes discussed (Zr(IV), Hf(IV), Th(IV), U(IV), U(V), Np(IV)). 17 O NMR spectroscopy subsequently provides insights into both the nature of metal–oxygen bonding, as well as the influence of unpaired f-electrons on the local environment of the oxygen nuclei.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Main and Satellite Features in the Ni 2p XPS of NiO

The origin and assignment of the complex main and satellite XPS features of the cations in ionic compounds has been the subject of extensive theoretical studies using different methods. There is agreement that within a molecular orbital model one needs to take into account different types of configurations. Specifically, those where a core electron is removed but no other configuration changes are made and those where in addition to ionization there are also shake or charge transfer changes to the ionic configuration. However, there are strong disagreements about the assignment of XPS features to these configurations. The present work is directed toward resolving the origin of main and satellite features for the Ni 2p XPS of NiO based on ab initio molecular orbital wavefunctions for a cluster model of NiO. A major problem in earlier ab initio XPS studies of ionic compounds has been the use of a common set of orbitals that was not able to properly describe all the ionic configurations that contribute to the full XPS spectra. This is resolved in the present work by using orbitals that are optimized for averages of the occupations of the different configurations that contribute to the XPS. The approach of using State Averaged orbitals is validated through comparisons between different averages and through use of higher order excitations in the wavefunctions for the ionic states. It represents a major extension of our earlier work on the main and satellite features of the Fe 2p XPS of Fe 2 O 3 and proves the reliability and the generality of the assignments of the character and origin of the different features of the XPS obtained with orbitals optimized for State Averages. These molecular orbital methods permit the characterization of the ionic states in terms of the importance of shake excitations and of the coupling of ionization of 2p 1/2 and 2p 3/2 spin-orbit split sub shells. Further, the work lays the foundation for definitive assignments of the character of main and satellite XPS features and points to their origin in the electronic structure of the material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gas-Phase Stability of Large Lanthanide:Ligand Clusters Evaluated Using Collision-Induced Dissociation

Introduction In the reprocessing of f-elements present in used nuclear fuels, a variety of diglycolamides (DGA’s) are used as extractants for actinide partitioning. In particular, the Actinide-Lanthanide Separation (ALSEP) process typically utilizes either the N,N,N’,N’-tetraoctyl diglycolamide (TODGA) or N,N,N',N'-tetra-2-ethylhexyl diglycolamide (T2EHDGA) extractant ligands following the partitioning of uranium and plutonium from used nuclear fuel. To better understand fundamental interactions in these processes, covalent bonding of several f-elements with diglycolamides, primarily TODGA, is investigated in the gas phase using nanospray ionization and a quadrupole time-of-flight mass spectrometer. Further, analysis of the identity and relative strength of the cluster is enabled by MS2 isolation and collision induced dissociation. Methods Metal ion cluster analysis was completed using a Bruker (Billerica, MA, USA) mircOTOF-Q II quadrupole time-of-flight mass spectrometer with a CaptiveSpray nanospray ion source. Detection was accomplished using positive ionization mode. Metal: ligand solutions were assembled as 30 µM europium nitrate, samarium nitrate, cerium nitrate, or holmium nitrate and 3 µM DGA in acetonitrile or a 50:50 mixture of acetonitrile: isopropanol. Preliminary data The samarium cluster experiments yielded clusters with a samarium:TODGA ratio of up to 1:7 able to be isolated and evidence of greater ratios present in the mass spectrum. This is surprising, as metal clusters are not expected to have a coordination space able to accommodate this many TODGA ligands, due to its size and tridenticity. Collisional activation of [Sm(TODGA)3]3+ suggested loss of a TODGA radical cation, in addition to ligand fragmentation. In contrast, activation of clusters with higher Sm:TODGA ratios resulted in loss of entire ligands, with no evidence of fragmentation. A lower collision energy was required to remove ligands as the number of bound TODGAs increased, suggesting that in larger clusters, ligands are more delicately complexed to the metal. In addition, several clusters were observed with the composition [Sm(NO3)x(TODGA)n x]+3 x. With a single nitrate ion, clusters with up to six TODGAs were able to be isolated. In a similar pattern to the samarium clusters containing only TODGA, less collision energy was required to eliminate one or more TODGAs with increasing size. Clusters with composition [Sm(NO3)(TODGA)n-1]2+ appeared in lower abundance and were more collisionally stable than [Sm(TODGA)n]3+ clusters. With two nitrate ions, only clusters with a single TODGA were able to be isolated. Analogous europium experiments resulted in similar clusters. Ratios of up to 1:7 Eu:TODGA and clusters with one nitrate and up to five TODGAs were isolated. In clusters with two nitrate ions, only one or two TODGAs were observed to be bound. Similar to samarium, MS2 experiments with the Eu clusters suggested that larger clusters required less collision energy to eliminate TODGA. Europium clusters with the composition [Eu(NO3)(TODGA)n-1]2+ were observed in greater abundance and with greater stability than the equivalent cluster with the composition [Eu(TODGA)n]3+. Novel Aspect These are the first reported Ln:TODGA clusters, allowing us to begin to investigate intrinsic complexation of lanthanides with process-relevant ligands.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, in this work, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1–, and 2–, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Nanosized {Ni II 18 } Cluster with a ‘Flying Saucer’ Topology Exhibiting Slow Relaxation of Magnetisation Phenomena at Both 15 K and 1.3 K

A high-nuclearity {Ni 18 } complex (1) with a unique ‘flying saucer’ motif has been prepared from the organic chelate, α-methyl-2-pyridine-methanol (mpmH), in conjunction with bridging azido (N 3 - ) and peroxido (O 2 2- ) ligands. Magnetic susceptibility measurements revealed the presence of both ferro- and antiferromagnetic exchange interactions between the metal centres in 1, and the stabilization of spin states with appreciable S values at two different temperature regimes. The end-on bridging azido and alkoxido groups are in all likelihood the ferromagnetic mediators, while the η 3 :η 3 :μ 6 -bridging peroxides most likely promote the antiparallel alignment of the metals’ spin vectors, yielding an overall non-zero spin ground state for the centrosymmetric compound 1. Furthermore, the {Ni 18 } nanosized cluster behaves as a single-molecule magnet, exhibiting magnetic hysteresis at low temperatures and two relaxation processes at 15 K and 1.3 K, a very rare phenomenon in polynuclear magnetic 3d-metal clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-Principles Simulations of CuCl in High-Temperature Water Vapor

Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H 2 O) n , where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta2014, 129, 33–53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids2018, 2018, 4279124]. Here, we have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H 2 O–Cu–Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗