SEARCH · Engineering Papers

Results for “bond angles”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗

DOE OSTI2020

Searching for correlations between vibrational spectral features and structural parameters of silicate glass network

Infrared (IR) and Raman spectroscopic features of silicate glasses are usually interpreted based on the analogy with those of smaller molecules, molecular clusters, or crystalline counterparts; this study tests the accuracy and validity of these widely cited peak assignment schemes by comparing vibrational spectral features with bond parameters of the glass network created by molecular dynamics (MD) simulations. A series of sodium silicate glasses with compositions of [Na 2 O] x [Al 2 O 3 ] 2 [SiO 2 ] 98- x with x = 7, 12, 17, and 22 were synthesized and analyzed with IR and Raman. A silica glass substrate and a crystalline quartz were also analyzed for comparison. Glass structures with the same compositions were generated with MD simulations using three types of potentials: fixed partial charge pairwise (Teter), partial diffuse charge potential (MGFF), and bond order-based charge transfer potential (ReaxFF). The comparison of simulated and experimental IR spectra showed that, among these three potentials tested, ReaxFF reproduces the concentration dependence of spectral features closest to the experimentally observed trend. Hence, the bond length and angle distributions as well as Si– Q n species and ring size distributions of silica and sodium silicate glasses were obtained from ReaxFF-MD simulations and further compared with the peak assignment or deconvolution schemes—which have been widely used since 1970s and 1980s—(a) correlation between the IR peak position in the Si–O stretch region (1050-1120 cm -1 ) and the Si–O–Si bond angle; (b) deconvolution of the Raman bands in the Si–O stretch region with the Q n speciation; and (c) assignment of the Raman bands in the 420-600 cm -1 region to the bending modes of (SiO) n rings with different sizes (typically, n = 3-6). The comparisons showed that none of these widely used methods is congruent with the bond parameters or structures of silicate glass networks produced via ReaxFF-MD simulations. This result invokes that the adequacy of these spectral interpretation methods must be questioned. Alternative interpretations are proposed, which are to be tested independently in future studies.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P 4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C–C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp 3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P 4 to ‘break’ their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A Diazo Linker Ligand Promotes Flexibility and Induced Fit Binding in a Microporous Copper Coordination Network

Abstract Flexible organic linkers represent an intuitive and effective strategy to design flexible metal–organic materials. We report herein a systematic study concerning the effect of varying the central bond of mixed pyridyl‐benzoate linkers, L, upon the flexibility of three isostructuralkddtopology microporous coordination networks (CNs) of formula ML 2 :X‐kdd‐1‐Cu,1= L = (E)‐4‐(pyridin‐4‐yldiazenyl)benzoate;X‐kdd‐2‐Cu,2= L = (E)‐4‐(2‐(pyridin‐4‐yl)vinyl)benzoate; the previously reportedX‐kdd‐3‐Cu,3= L = 4‐(pyridin‐4‐ylethynyl)benzoate. As revealed by single crystal x‐ray diffraction (SCXRD) and gas sorption studies,X‐kdd‐1‐Cu, exhibited gate‐opening during CO 2 and hydrocarbon (C2 and C8) sorption experiments whereas the other two CNs did not. Insight into these phase transformations was gained from in situ variable‐pressure and variable temperature powder X‐ray diffraction (PXRD), SCXRD, and modeling. Rotation of ligand1around the diazo bond, torsion angle changes between phenyl and carboxylate moieties, and deformation of the Cu‐based rod building blocks enabled activatedX‐kdd‐1‐Cuto form new phases with C8 isomers and CH 2 Cl 2 , CH 2 Cl 2 inducing contraction of the activated phase. Computational studies suggest that1enables flexibility thanks to its lower barrier of deformation versus2or3. This study teaches that diazo moieties could offer a general strategy to enhance the flexibility of CNs.

Chemistry↗

DOE OSTI2025

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗

DOE OSTI2024

Understanding Strain and Failure of a Knot in Polyethylene Using Molecular Dynamics with Machine-Learned Potentials

A neural network potential (NNP) has been developed by fitting to ab initio electronic structure data on hydrocarbons and is used to study failure of linear and knotted polyethylene (PE) chains. A linear PE chain must be highly strained before breaking as the stress is equally distributed across the chain. In contrast, the stress in a PE chain with a 31 or overhand knot, accumulates at the knot’s entrance/exit. We find the strain energy is greatest when the bond length and angle are strained simultaneously, and that the knot weakens the chain by increasing the variance of the C–C–C angle, thereby allowing rupture at lower bond strains. Here, we extend our analysis to both 51 and 52 knots and find that both break at the entrance/exit of a loop. Notably, molecular scale PE knots exhibit many of the same characteristics as knots in a macroscopic rope, with stick–slip phenomena upon tightening and similar points of failure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Gaussian approximation potential for amorphous Si : H

Hydrogenation of amorphous silicon (a–Si : H) is critical for reducing defect densities, passivating midgap states and surfaces, and improving photoconductivity in silicon-based electro-optical devices. Modeling the atomic-scale structure of this material is critical to understanding these processes, which in turn is needed to describe c–Si/a–Si : H heterojunctions that are at the heart of modern solar cells with world-record efficiency. Density functional theory (DFT) studies achieve the required high accuracy but are limited to moderate system sizes of 100 atoms or so by their high computational cost. Simulations of amorphous materials have been hindered by this high cost because large structural models are required to capture the medium-range order that is characteristic of such materials. Empirical potential models are much faster, but their accuracy is not sufficient to correctly describe the frustrated local structure. Data-driven, machine-learned interatomic potentials have broken this impasse and have been highly successful in describing a variety of amorphous materials in their elemental phase. Here, we extend the Gaussian approximation potential (GAP) for silicon by incorporating the interaction with hydrogen, thereby significantly improving the degree of realism with which amorphous silicon can be modeled. We show that our Si : H GAP enables the simulation of hydrogenated silicon with an accuracy very close to DFT but with computational expense and run times reduced by several orders of magnitude for large structures. Here, we demonstrate the capabilities of the Si : H GAP by creating models of hydrogenated liquid and amorphous silicon and showing that their energies, forces, and stresses are in excellent agreement with DFT results, and their structure as captured by bond and angle distributions are in agreement with both DFT and experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Ig-VAE: Generative modeling of protein structure by direct 3D coordinate generation

While deep learning models have seen increasing applications in protein science, few have been implemented for protein backbone generation—an important task in structure-based problems such as active site and interface design. We present a new approach to building class-specific backbones, using a variational auto-encoder to directly generate the 3D coordinates of immunoglobulins. Our model is torsion- and distance-aware, learns a high-resolution embedding of the dataset, and generates novel, high-quality structures compatible with existing design tools. We show that the Ig-VAE can be used with Rosetta to create a computational model of a SARS-CoV2-RBD binder via latent space sampling. We further demonstrate that the model’s generative prior is a powerful tool for guiding computational protein design, motivating a new paradigm under which backbone design is solved as constrained optimization problem in the latent space of a generative model.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗

DOE OSTI2024

Optimizing Aluminum Bonding: Exploring Surface Roughness and Contact Angle Effects through Plasma and Acid Etching Contrasts

The use of adhesives to replace traditional techniques like welding and riveting has been on the rise in various industries as increasingly complex geometries and situations for bonding is required. It’s well known in the scientific community that an adhesive bond is only as efficient as the surface preparation used beforehand. Plasma treatments have been shown to gently clean the surface by only removing the top-most layer from the adherend and “charging” the surface instead of berating it. Whereas chemical treatments “etch” the surface in hopes of increasing surface area by increasing surface roughness. This method can be difficult to control quantitatively and can generate waste leading to serious environmental concerns. This study aims to address the encountered differences between plasma treatments and acid etching on aluminum substrates for adhesive bonding. Treating the surface of adherends plays a crucial role in promoting a good bond between adhesive and adherend. In this study, the effects of plasma and acid treatments are categorized using surface roughness measurements as well as contact angle which have been shown to correlate to mechanical strength. When evaluating different surface treatments, it is notable to mention that not only surface roughness and contact angle play a role in determining which method is most effective but also the consequences of the treatment itself. It was found that although chemical etching was the most effective in increasing roughness and decreasing contact angle, it generates a large amount of waste and is not sustainable compared to plasma work which uses a renewable process. The results from this work provide a deeper understanding of the relationship between surface treatments and surface roughness/contact angle and how this relates to a stronger adhesive bond.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Optimizing Aluminum Bonding: Exploring Surface Roughness and Contact Angle Effects through Plasma and Acid Etching Contrasts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Nonlinear deformation and elasticity of BCC refractory metals and alloys

Application of isotropic pressure or uniaxial strain alters the elastic properties of materials; sufficiently large strains can drive structural transformations. Linear elasticity describes stability against infinitesimal strains, while nonlinear elasticity describes the response to finite deformations. Here, it was previously shown that uniaxial strain along [100] drives refractory metals and alloys towards mechanical instabilities. These include an extensional instability, and a symmetry-breaking orthorhombic distortion caused by a Jahn-Teller-Peierls instability that splays the cubic lattice vectors. Here we analyze these transitions in depth. Eigenvalues and eigenvectors of the Wallace tensor identify and classify linear instabilities in the presence of strain. We show that both instabilities are discontinuous, leading to discrete jumps in the lattice parameters. We provide physical intuition for the instabilities by analyzing the changes in first-principles energy, stress, bond lengths, and angles upon application of strain. Electronic band structure calculations show differential occupation of bonding and antibonding orbitals, driven by the changing bond lengths and leading to the structural transformations. Strain thresholds for these instabilities depend on the valence electron count.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Uranyl Naphthylsalophen and Pyrasal Complexes: Oxo Ligands Acting as Hydrogen Bond Acceptors in the Solid State

Uranium is most stable when it is exposed to oxygen or water in its +6 oxidation state as the uranyl (UO 2 2+ ) ion. This ion is subsequently particularly stable and very resistant to functionalization due to the inverse trans effect. Uranyl oxo ligands are typically not considered good hydrogen bond acceptors due to their weak Lewis basicity; however, the ligands bound in the equatorial plane greatly affect the strength of the oxo ligands’ hydrogen bonding. Here, in this work, new naphthylsalophen and pyrasal complexes of uranium were synthesized and crystallized for characterization in the solid state. The bond lengths and angles of the uranyl ion and the ligand conformation are compared. In the solid state, one of the pyrasal complexes showed a hydrogen bond directly from a water molecule to the uranyl oxo ligand, which resulted in an asymmetric lengthening of the U–O yl bonds from 1.789 to 1.862 Å and 1.784 to 1.844 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Synchrotron X-ray diffraction of pyrolusite (MnO 2 ) and rutile (TiO 2 ) during heating to ~1000 °C

To compare thermal expansion behaviors in isomorphic structures with different redox behaviors, we separately heated pyrolusite (Mn 4+ O 2 ) and rutile (Ti 4+ O 2 ) powders from 25 to ~1000 °C at a rate of 0.05 °C/s at ambient pressure and analyzed them using synchrotron X-ray diffraction and Rietveld refinement. The pyrolusite exhibited two reduction-induced phase decompositions. Between 531 and 583 °C, pyrolusite decomposed to bixbyite (Mn 3+ 2 O 3 ), and between 972 and 998 °C, bixbyite decomposed to hausmannite (Mn 2+ Mn 3+ 2 O 4 ). The rutile experienced no phase decompositions or transitions over this same temperature interval. The refined bond lengths and angles for pyrolusite showed that the Mn 4+ coordination octahedron became more distorted near the phase decomposition to bixbyite due to the approach of two coordinating oxygen atoms. The 3d electrons of each Mn 4+ ion in pyrolusite make π bonds with the 2p electrons of the surrounding O 2- ions, whereas the Ti 4+ of rutile has no 3d electrons. Thus, π bonding between octahedral Mn 4+ and the surrounding O 2- anions increases the strength of the Mn–O bonds of pyrolusite relative to the Ti–O bonds of rutile. However, we observed a small decrease in O–O distances in pyrolusite before the decomposition, suggesting that at high temperature, increased π bonding between adjacent O anions precedes the release of O 2 and the reduction of Mn 4+ to Mn 3+ . Analogous behavior was observed for bixbyite before its reductive phase decomposition to hausmannite. In contrast, no anomalous changes in the O–O distances occurred for rutile. Our X-ray diffraction analyses provided accurate thermal expansion coefficients for these materials over a broader temperature range than reported in previous studies. The Mn–O bond lengths in pyrolusite were shorter and stronger than the Ti–O bonds of rutile; for example, at 60 °C, bond distances in pyrolusite refined to 1.8854(6) Å whereas the bond lengths of rutile were 1.9575(4) Å. Accordingly, the axial and volumetric thermal expansion coefficients refined for rutile were greater than those for pyrolusite. Further, the axial coefficients of thermal expansion (CTE) for pyrolusite were 96% (a-axis) and 69% (c-axis) of the corresponding values for rutile, and the volumetric CTE for pyrolusite was 85% of that for rutile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Non‐Linear Optical Properties of the ( RE ) 3 CuGeS 7 Family of Compounds

Abstract Non‐linear optical materials must possess a balanced combination of laser‐induced damage threshold (LDT) and second‐harmonic generation (SHG) and be phase matchable. In our previous work, chiral and polar La 3 CuGeS 7 was identified as a promising non‐linear optical material. Herein, we report the optimization of non‐linear optical properties through replacement of La with smaller lanthanides. It is determined that Gd 3 CuGeS 7 exhibits the best combination of SHG (1.6× AgGaS 2 at 88–105 μm particle size) and LDT (3× AgGaS 2 , 89 MW/cm 2 ) and is phase matchable. Based on changes in metal‐sulfur bond lengths and angles, we further propose structural optimization through solid‐solution formation and doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

One Electron Multiple Proton Transfer in Model Organic Donor–Acceptor Systems: Implications for High-Frequency EPR

EPR spectroscopy is an important spectroscopic method for identification and characterization of radical species involved in many biological reactions. The tyrosyl radical is one of the most studied amino acid radical intermediates in biology. Often in conjunction with histidine residues, it is involved in many fundamental biological electron and proton transfer processes, such as in the water oxidation in photosystem II. As biological processes are typically extremely complicated and hard to control, molecular bio-mimetic model complexes are often used to clarify the mechanisms of the biological reactions. In this work, we present theoretical calculations to investigate the sensitivity of magnetic resonance parameters to proton-coupled electron transfer events, as well as conformational substates of the molecular constructs which mimic the tyrosine-histidine (Tyr-His) pairs found in a large variety of proteins. Upon oxidation of the phenol, the Tyr analog, these complexes can perform not only one-electron one-proton transfer (EPT), but also one-electron two-proton transfers (E2PT). It is shown that in aprotic environment the g X -components of the electronicg-tensor are extremely sensitive to the first proton transfer from the phenoxyl oxygen to the imidazole nitrogen (EPT product), leading to a significant increase of the g X -value of up to 0.003, but are not sensitive to the second proton transfer (E2PT). In the latter case, the change of the g X -value is much smaller (ca. 0.0001), which is too small to be distinguished even by high-frequency EPR. The 14 N hyperfine values are also too similar to allow differentiation between the different protonation states in EPT and E2PT. The magnetic resonance parameters were also calculated as a function of the rotation angles around single bonds. It was demonstrated that rotation of the phenoxyl group results in large positive changes (> 0.001) in the g X -values. Analysis of the data reveals that the main source of these changes is related to the strength of the H-bond between phenoxyl oxygen and the proton(s) on N 1 and N 2 positions of the imidazole.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Physical structure of constitutional isomers influences antiproliferation activity of thiosemicarbazone-alkylthiocarbamate copper complexes

A series of hybrid thiosemicarbazone-alkylthiocarbamate copper complexes with similar electronic environments but distinct physical structures have been prepared, characterized, and evaluated for antiproliferation activity. The complexes include the constitutional isomers (1-phenylpropane-1-imine-(O-ethylthiocarbamato)-2-one-(N-methylthiosemicarbazonato))copper(II) (CuL 1 ) and (1-phenylpropane-1-one-(N-methylthiosemicarbazonato)-2-imine-(O-ethylthiocarbamato))copper(II) (CuL 2 ) along with (1-propane-1-imine-(O-ethylthiocarbamato)-2-one-(N-methylthiosemicarbazonato))copper(II) (CuL 3 ). Complexes CuL 1 and CuL 2 differ in the positions of the pendent thiosemicarbazone (TSC) and alkylthiocarbamate (ATC) moieties on the 1-phenylpropane backbone. Complex CuL 3 employs a propane backbone with the TSC in the 2-position as in CuL 1 . The isomer pair CuL 1 and CuL 2 have equivalent electronic environments with indistinguishable Cu II/I potentials (E 1/2 = -0.86 V vs. ferrocenium/ferrocene) and electron paramagnetic resonance (EPR) spectra (g ∥ = 2.26, g ⊥ = 2.08). The electronic structure of CuL 3 has a similar E 1/2 of -0.84 V and identical EPR parameters to CuL 1, 2 . Single crystal X-ray diffraction studies confirm a consistent donor environment with no substantial variation in the Cusingle bondN or Cusingle bondS bond distances and angles between the complexes. The antiproliferation activities of the CuL 1–3 were evaluated against the lung adenocarcinoma cell line (A549) and nonmalignant lung fibroblast cell line (IMR-90) using the MTT assay. CuL 1 had the highest A549 activity ( A549 EC 50 = 0.065 μM) and selectivity ( IMR-90 EC 50 / A549 EC 50 = 20). The constitutional isomer CuL 2 displayed decreased A549 activity (0.18 μM) and selectivity (10.6). The complex CuL 3 displayed activity (0.009 μM) similar to CuL 1 but with a lack of selectivity (1.0). Cellular copper loading determined by ICP-MS was consistent with the activity and selectivity trends. The complexes CuL 1–3 did not induce reactive oxygen species (ROS) generation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗