Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ALKANE”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Improving the efficiency of configurational-bias Monte Carlo: A density-guided method for generating bending angle trials for linear and branched molecules

A new method has been developed to generate bending angle trials to improve the acceptance rate and the speed of configurational-bias Monte Carlo. Whereas traditionally the trial geometries are generated from a uniform distribution, in this method we attempt to use the exact probability density function so that each geometry generated is likely to be accepted. In actual practice, due to the complexity of this probability density function, a numerical representation of this distribution function would be required. This numerical table can be generated a priori from the distribution function. This method has been tested on a united-atom model of alkanes including propane, 2-methylpropane, and 2,2-dimethylpropane, that are good representatives of both linear and branched molecules. It has been shown from these test cases that reasonable approximations can be made especially for the highly branched molecules to reduce drastically the dimensionality and correspondingly the amount of the tabulated data that is needed to be stored. Despite these approximations, the dependencies between the various geometrical variables can be still well considered, as evident from a nearly perfect acceptance rate achieved. For all cases, the bending angles were shown to be sampled correctly by this method with an acceptance rate of at least 96% for 2,2-dimethylpropane to more than 99% for propane. Since only one trial is required to be generated for each bending angle (instead of thousands of trials required by the conventional algorithm), this method can dramatically reduce the simulation time. The profiling results of our Monte Carlo simulation code show that trial generation, which used to be the most time consuming process, is no longer the time dominating component of the simulation.

2-2-DIMETHYLPROPANE↗
DOE OSTI2014

NMR characterization of hydrocarbon adsorption on calcite surfaces: A first principles study

The electronic and coordination environment of minerals surfaces, as calcite, are very difficult to characterize experimentally. This is mainly due to the fact that there are relatively few spectroscopic techniques able to detect Ca{sup 2+}. Since calcite is a major constituent of sedimentary rocks in oil reservoir, a more detailed characterization of the interaction between hydrocarbon molecules and mineral surfaces is highly desirable. Here we perform a first principles study on the adsorption of hydrocarbon molecules on calcite surface (CaCO{sub 3} (101{sup ¯}4)). The simulations were based on Density Functional Theory with Solid State Nuclear Magnetic Resonance (SS-NMR) calculations. The Gauge-Including Projector Augmented Wave method was used to compute mainly SS-NMR parameters for {sup 43}Ca, {sup 13}C, and {sup 17}O in calcite surface. It was possible to assign the peaks in the theoretical NMR spectra for all structures studied. Besides showing different chemical shifts for atoms located on different environments (bulk and surface) for calcite, the results also display changes on the chemical shift, mainly for Ca sites, when the hydrocarbon molecules are present. Even though the interaction of the benzene molecule with the calcite surface is weak, there is a clearly distinguishable displacement of the signal of the Ca sites over which the hydrocarbon molecule is located. A similar effect is also observed for hexane adsorption. Through NMR spectroscopy, we show that aromatic and alkane hydrocarbon molecules adsorbed on carbonate surfaces can be differentiated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A model-free temperature-dependent conformational study of n-pentane in nematic liquid crystals

The proton NMR spectra of n-pentane orientationally ordered in two nematic liquid-crystal solvents are studied over a wide temperature range and analysed using covariance matrix adaptation evolutionary strategy. Since alkanes possess small electrostatic moments, their anisotropic intermolecular interactions are dominated by short-range size-and-shape effects. As we assumed for n-butane, the anisotropic energy parameters of each n-pentane conformer are taken to be proportional to those of ethane and propane, independent of temperature. The observed temperature dependence of the n-pentane dipolar couplings allows a model-free separation between conformer degrees of order and conformer probabilities, which cannot be achieved at a single temperature. In this way for n-pentane 13 anisotropic energy parameters (two for trans trans, tt, five for trans gauche, tg, and three for each of gauche{sub +} gauche{sub +}, pp, and gauche{sub +} gauche{sub −}, pm), the isotropic trans-gauche energy difference E{sub tg} and its temperature coefficient E{sub tg}{sup ′} are obtained. The value obtained for the extra energy associated with the proximity of the two methyl groups in the gauche{sub +} gauche{sub −} conformers (the pentane effect) is sensitive to minute details of other assumptions and is thus fixed in the calculations. Conformer populations are affected by the environment. In particular, anisotropic interactions increase the trans probability in the ordered phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Molecular dynamics study on condensation/evaporation coefficients of chain molecules at liquid–vapor interface

The structure and thermodynamic properties of the liquid–vapor interface are of fundamental interest for numerous technological implications. For simple molecules, e.g., argon and water, the molecular condensation/evaporation behavior depends strongly on their translational motion and the system temperature. Existing molecular dynamics (MD) results are consistent with the theoretical predictions based on the assumption that the liquid and vapor states in the vicinity of the liquid–vapor interface are isotropic. Additionally, similar molecular condensation/evaporation characteristics have been found for long-chain molecules, e.g., dodecane. It is unclear, however, whether the isotropic assumption is valid and whether the molecular orientation or the chain length of the molecules affects the condensation/evaporation behavior at the liquid–vapor interface. In this study, MD simulations were performed to study the molecular condensation/evaporation behavior of the straight-chain alkanes, i.e., butane, octane, and dodecane, at the liquid–vapor interface, and the effects of the molecular orientation and chain length were investigated in equilibrium systems. The results showed that the condensation/evaporation behavior of chain molecules primarily depends on the molecular translational energy and the surface temperature and is independent of the molecular chain length. Furthermore, the orientation at the liquid–vapor interface was disordered when the surface temperature was sufficiently higher than the triple point and had no significant effect on the molecular condensation/evaporation behavior. The validity of the isotropic assumption was confirmed, and we conclude that the condensation/evaporation coefficients can be predicted by the liquid-to-vapor translational length ratio, even for chain molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Thermodynamic and kinetic anisotropies in octane thin films

Confinement breaks the translational symmetry of materials, making all thermodynamic and kinetic quantities functions of position. Such symmetry breaking can be used to obtain configurations that are not otherwise accessible in the bulk. Here, we use computer simulations to explore the effect of substrate-liquid interactions on thermodynamic and kinetic anisotropies induced by a solid substrate. We consider n-octane nano-films that are in contact with substrates with varying degrees of attraction, parameterized by an interaction parameter ϵ{sub S}. Complete freezing of octane nano-films is observed at low temperatures, irrespective of ϵ{sub S}, while at intermediate temperatures, a frozen monolayer emerges at solid-liquid and vapor-liquid interfaces. By carefully inspecting the profiles of translational and orientational relaxation times, we confirm that the translational and orientational degrees of freedom are decoupled at these frozen monolayers. At sufficiently high temperatures, however, free interfaces and solid-liquid interfaces close to loose (low-ϵ{sub S}) substrates undergo “pre-freezing,” characterized by mild peaks in several thermodynamic quantities. Two distinct dynamic regimes are observed at solid-liquid interfaces. The dynamics is accelerated in the vicinity of loose substrates, while sticky (high-ϵ{sub S}) substrates decelerate dynamics, sometimes by as much as two orders of magnitude. These two distinct dynamical regimes have been previously reported by Haji-Akbari and Debenedetti [J. Chem. Phys. 141, 024506 (2014)] for a model atomic glass-forming liquid. We also confirm the existence of two correlations—proposed in the above-mentioned work—in solid-liquid subsurface regions of octane thin films, i.e., a correlation between atomic density and normal stress, and between atomic translational relaxation time and lateral stress. Finally, we inspect the ability of different regions of an octane film to explore the potential energy landscape by performing inherent structure calculations, and observe no noticeable difference between the free surface and the bulk in efficiently exploring the potential energy landscape. This is unlike the films of model atomic glass formers that tend to sample their respective landscape more efficiently at free surfaces. We discuss the implications of this finding to the ability of octane—and other n-alkanes—to form ultrastable glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Anharmonic modeling of the conformation-specific IR spectra of ethyl, n-propyl, and n-butylbenzene

In this paper, conformation-specific UV-IR double resonance spectra are presented for ethyl, n-propyl, and n-butylbenzene. With the aid of a local mode Hamiltonian that includes the effects of stretch-scissor Fermi resonance, the spectra can be accurately modeled for specific conformers. These molecules allow for further development of a first principles method for calculating alkyl stretch spectra. Across all chain lengths, certain dihedral patterns impart particular spectral motifs at the quadratic level. However, the anharmonic contributions are consistent from molecule to molecule and conformer to conformer. This transferability of anharmonicities allows for the Hamiltonian to be constructed from only a harmonic frequency calculation, reducing the cost of the model. Finally, the phenyl ring alters the frequencies of the CH 2 stretches by about 15 cm -1 compared to their n-alkane counterparts in trans configurations. Conformational changes in the chain can lead to shifts in frequency of up to 30 cm -1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

On the determination of the crystal-vapor surface free energy, and why a Gaussian expression can be accurate for a system far from Gaussian

The crystal-vapor surface free energy γ is an important physical parameter governing physical processes, such as wetting and adhesion. We explore exact and approximate routes to calculate γ based on cleaving an intact crystal into non-interacting sub-systems with crystal-vapor interfaces. We do this by turning off the interactions, ΔV, between the sub-systems. Using the soft-core scheme for turning off ΔV, we find that the free energy varies smoothly with the coupling parameter λ, and a single thermodynamic integration yields the exact γ. We generate another exact method, and a cumulant expansion for γ by expressing the surface free energy in terms of an average of e{sup −βΔV} in the intact crystal. The second cumulant, or Gaussian approximation for γ is surprisingly accurate in most situations, even though we find that the underlying probability distribution for ΔV is clearly not Gaussian. We account for this fact by developing a non-Gaussian theory for γ and find that the difference between the non-Gaussian and Gaussian expressions for γ consist of terms that are negligible in many situations. Exact and approximate methods are applied to the (111) surface of a Lennard-Jones crystal and are also tested for more complex molecular solids, the surface of octane and nonadecane. Alkane surfaces were chosen for study because their crystal-vapor surface free energy has been of particular interest for understanding surface freezing in these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Reduced scaling Hilbert space variational Monte Carlo

We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through various tensor decomposition methods, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons when using a minimal basis. Further study will be needed to assess where the crossover occurs in more compact molecular geometries and larger basis sets and to explore how in that context the crossover can be accelerated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Vibrational coherence transfer illuminates dark modes in models of the FeFe hydrogenase active site

Within the conceptual framework of Redfield theory, the optical response function arises from the dynamical evolution of the system’s density operator, where nonunitary relaxation is encoded in the Redfield relaxation superoperator. In the conventional approach, the so-called secular approximation neglects terms that induce transitions between distinct coherences and among coherences and populations. The rationale is that these nonsecular terms are small in comparison to the far more dominant population relaxation and coherence dephasing contributions. Since two-dimensional infrared (2D-IR) spectroscopy has significant contributions arising from population relaxation and transfer pathways, it can be challenging to isolate signatures of the nonsecular relaxation. We report om this paper that in three diiron dithiolate hexacarbonyl complexes that serve as small-molecule models of the [FeFe] hydrogenase H-cluster subsite, a fortuitous vibrational energy structure enables direct and clear signatures of vibrational coherence transfer in alkane solution. This finding holds promise towards developing a molecularly detailed understanding of the mechanism of vibrational coherence transfer processes, thanks to the ease of synthesizing derivatives based on the chemical modularity of these well studied diiron compounds. In addition to the fundamental need to characterize coherence transfer in molecular spectroscopy, we find in this set of molecules a practical utility for the nonsecular dynamics: the ability to determine the frequency of an IR-inactive mode. A coherence generated during the waiting time of the 2D-IR measurement transfers to a coherence involving the single dark CO stretching mode, which modulates some peak amplitudes in the 2D spectrum, revealing its transient excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C160H322 n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580–6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Kumar, Ashutosh (ORCID:0000000175896030)↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Characterization of a new ultra-high pressure shock tube facility for combustion and propulsion studies

A new shock tube facility has been designed, constructed, and characterized at the University of Central Florida. This facility is capable of withstanding pressures of up to 1000 atm, allowing for combustion diagnostics of extreme conditions, such as in rocket combustion chambers or in novel power conversion cycles. For studies with toxic gas impurities, the high initial pressures required the development of a gas delivery system to ensure the longevity of the facility and the safety of the personnel. Data acquisition and experimental propagation were implemented with remote access to ensure safety, paired with a LabVIEW- and Python-based user interface. Thus far, test pressures of 270 atm, blast pressures of 730 atm, and temperatures approaching 10 000 K have been achieved. The extreme limitations of this facility allow for emission spectroscopy to be performed during the oxidation of fuel mixtures, e.g., alkanes diluted in argon and carbon dioxide. Ignition delay times were determined and compared to simulations using chemical kinetic mechanisms. The design, experimental procedures, processes of analysis, and uncertainty determination are outlined, and typical pressure profiles are compared with a new gas dynamics solver and empirical correlations developed across multiple shock tube facilities. Preliminary reactive mixture analyses are included with further investigation of the mixtures outlined.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Interpretation of autoencoder-learned collective variables using Morse–Smale complex and sublevelset persistent homology: An application on molecular trajectories

Dimensionality reduction often serves as the first step toward a minimalist understanding of physical systems as well as the accelerated simulations of them. In particular, neural network-based nonlinear dimensionality reduction methods, such as autoencoders, have shown promising outcomes in uncovering collective variables (CVs). However, the physical meaning of these CVs remains largely elusive. In this work, we constructed a framework that (1) determines the optimal number of CVs needed to capture the essential molecular motions using an ensemble of hierarchical autoencoders and (2) provides topology-based interpretations to the autoencoder-learned CVs with Morse–Smale complex and sublevelset persistent homology. Furthermore, this approach was exemplified using a series of n-alkanes and can be regarded as a general, explainable nonlinear dimensionality reduction method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Linear-scaling quadruple excitations in local pair natural orbital coupled-cluster theory

Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.

Auxiliary functions↗
DOE OSTI2016

Self-assembly of crystalline nanotubes from monodisperse amphiphilic diblock copolypeptoid tiles

Significance A fundamental challenge in materials science is to create synthetic nanoarchitectures that rival the structural complexity found in nature. A promising bioinspired approach is to synthesize sequence-defined polymer chains that fold into precise protein-like structures. Here, we discovered a family of sequence-defined, amphiphilic diblock copolypeptoids that spontaneously form crystalline, hollow nanotubes of uniform diameter, without a central hydrophobic core in aqueous solution. These low–molecular-weight polymers have two chemically distinct domains that are congruent in size and shape, which allows them to behave like molecular tile units that can form brick-like supramolecular lattices. This demonstrates the remarkable structure-directing influence of n -alkane and ethyleneoxy side chains, and that flexible polymers can serve as tile units that can assemble into precision supramolecular architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗