Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERMOLECULAR FORCES”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Rotational study of the NH{sub 3}–CO complex: Millimeter-wave measurements and ab initio calculations

The rotational spectrum of the van der Waals complex NH{sub 3}–CO has been measured with the intracavity OROTRON jet spectrometer in the frequency range of 112–139 GHz. Newly observed and assigned transitions belong to the K = 0–0, K = 1–1, K = 1–0, and K = 2–1 subbands correlating with the rotationless (j{sub k}){sub NH3} = 0{sub 0} ground state of free ortho-NH{sub 3} and the K = 0–1 and K = 2–1 subbands correlating with the (j{sub k}){sub NH3} = 1{sub 1} ground state of free para-NH{sub 3}. The (approximate) quantum number K is the projection of the total angular momentum J on the intermolecular axis. Some of these transitions are continuations to higher J values of transition series observed previously [C. Xia et al., Mol. Phys. 99, 643 (2001)], the other transitions constitute newly detected subbands. The new data were analyzed together with the known millimeter-wave and microwave transitions in order to determine the molecular parameters of the ortho-NH{sub 3}–CO and para-NH{sub 3}–CO complexes. Accompanying ab initio calculations of the intermolecular potential energy surface (PES) of NH{sub 3}–CO has been carried out at the explicitly correlated coupled cluster level of theory with single, double, and perturbative triple excitations and an augmented correlation-consistent triple zeta basis set. The global minimum of the five-dimensional PES corresponds to an approximately T-shaped structure with the N atom closest to the CO subunit and binding energy D{sub e} = 359.21 cm{sup −1}. The bound rovibrational levels of the NH{sub 3}–CO complex were calculated for total angular momentum J = 0–6 on this intermolecular potential surface and compared with the experimental results. The calculated dissociation energies D{sub 0} are 210.43 and 218.66 cm{sup −1} for ortho-NH{sub 3}–CO and para-NH{sub 3}–CO, respectively.

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

Rotational study of the CH{sub 4}–CO complex: Millimeter-wave measurements and ab initio calculations

The rotational spectrum of the van der Waals complex CH{sub 4}–CO has been measured with the intracavity OROTRON jet spectrometer in the frequency range of 110–145 GHz. Newly observed and assigned transitions belong to the K = 2–1 subband correlating with the rotationless j{sub CH4} = 0 ground state and the K = 2–1 and K = 0–1 subbands correlating with the j{sub CH4} = 2 excited state of free methane. The (approximate) quantum number K is the projection of the total angular momentum J on the intermolecular axis. The new data were analyzed together with the known millimeter-wave and microwave transitions in order to determine the molecular parameters of the CH{sub 4}–CO complex. Accompanying ab initio calculations of the intermolecular potential energy surface (PES) of CH{sub 4}–CO have been carried out at the explicitly correlated coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)-F12a] and an augmented correlation-consistent triple zeta (aVTZ) basis set. The global minimum of the five-dimensional PES corresponds to an approximately T-shaped structure with the CH{sub 4} face closest to the CO subunit and binding energy D{sub e} = 177.82 cm{sup −1}. The bound rovibrational levels of the CH{sub 4}–CO complex were calculated for total angular momentum J = 0–6 on this intermolecular potential surface and compared with the experimental results. The calculated dissociation energies D{sub 0} are 91.32, 94.46, and 104.21 cm{sup −1} for A (j{sub CH4} = 0), F (j{sub CH4} = 1), and E (j{sub CH4} = 2) nuclear spin modifications of CH{sub 4}–CO, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2003

Gibbs-Thomson Law for Singular Step Segments: Thermodynamics Versus Kinetics

Classical Burton-Cabrera-Frank theory presumes that thermal fluctuations are so fast that at any time density of kinks on a step is comparable with the reciprocal intermolecular distance, so that the step rate is about isotropic within the crystal plane. Such azimuthal isotropy is, however, often not the case: Kink density may be much lower. In particular, it was recently found on the (010) face of orthorhombic lysozyme that interkink distance may exceed 500-600 intermolecular distances. Under such conditions, Gibbs-Thomson law (GTL) may not be applicable: On a straight step segment between two corners, communication between the comers occurs exclusively by kink exchange. Annihilation between kinks of opposite sign generated at the comers results in the grain in step energy entering GTL. If the step segment length l much greater than D/v, where D and v are the kink diffusivity and propagation rate, respectively, the opposite kinks have practically no chance to annihilate and GTL is not applicable. The opposite condition of the GTL applicability, l much less than D/v, is equivalent to the requirement that relative supersaturation Delta(sub mu)/kT much less than alpha/l, where alpha is molecular size. Thus, GTL may be applied to a segment of 10(exp 3)alpha approx. 3 x 10(exp -5)cm approx 0.3 micron only if supersaturation is less than 0.1%, while practically used driving forces for crystallization are much larger. Relationships alternative to the GTL for different, but low, kink density have been discussed. They confirm experimental evidences that the Burton-Cabrera-Frank theory of spiral growth is growth rates twice as low as compared to the observed figures. Also, application of GTL results in unrealistic step energy while suggested kinetic law give reasonable figures.

Chernov, A. A.↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

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

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

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

Combining density functional and incremental post-Hartree-Fock approaches for van der Waals dominated adsorbate-surface interactions: Ag{sub 2}/graphene

A combined density functional (DFT) and incremental post-Hartree-Fock (post-HF) approach, proven earlier to calculate He-surface potential energy surfaces [de Lara-Castells et al., J. Chem. Phys. 141, 151102 (2014)], is applied to describe the van der Waals dominated Ag{sub 2}/graphene interaction. It extends the dispersionless density functional theory developed by Pernal et al. [Phys. Rev. Lett. 103, 263201 (2009)] by including periodic boundary conditions while the dispersion is parametrized via the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)]. Starting with the elementary cluster unit of the target surface (benzene), continuing through the realistic cluster model (coronene), and ending with the periodic model of the extended system, modern ab initio methodologies for intermolecular interactions as well as state-of-the-art van der Waals-corrected density functional-based approaches are put together both to assess the accuracy of the composite scheme and to better characterize the Ag{sub 2}/graphene interaction. The present work illustrates how the combination of DFT and post-HF perspectives may be efficient to design simple and reliable ab initio-based schemes in extended systems for surface science applications.

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

Rotational Spectroscopy of the NH{sub 3}–H{sub 2} Molecular Complex

We report the first high resolution spectroscopic study of the NH{sub 3}–H{sub 2} van der Waals molecular complex. Three different experimental techniques, a molecular beam Fourier transform microwave spectrometer, a millimeter-wave intracavity jet OROTRON spectrometer, and a submillimeter-wave jet spectrometer with multipass cell, were used to detect pure rotational transitions of NH{sub 3}–H{sub 2} in the wide frequency range from 39 to 230 GHz. Two nuclear spin species, ( o )-NH{sub 3}–( o )-H{sub 2} and ( p )-NH{sub 3}–( o )-H{sub 2}, have been assigned as carriers of the observed lines on the basis of accompanying rovibrational calculations performed using the ab initio intermolecular potential energy surface (PES) of Maret et al. The experimental spectra were compared with the theoretical bound state results, thus providing a critical test of the quality of the NH{sub 3}–H{sub 2} PES, which is a key issue for reliable computations of the collisional excitation and de-excitation of ammonia in the dense interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Tailoring the Weight of Surface and Intralayer Edge States to Control LUMO Energies

Abstract The energies of the frontier molecular orbitals determine the optoelectronic properties in organic films, which are crucial for their application, and strongly depend on the morphology and supramolecular structure. The impact of the latter two properties on the electronic energy levels relies primarily on nearest‐neighbor interactions, which are difficult to study due to their nanoscale nature and heterogeneity. Here, an automated method is presented for fabricating thin films with a tailored ratio of surface to bulk sites and a controlled extension of domain edges, both of which are used to control nearest‐neighbor interactions. This method uses a Langmuir–Schaefer‐type rolling transfer of Langmuir layers (rtLL) to minimize flow during the deposition of rigid Langmuir layers composed of π‐conjugated molecules. Using UV–vis absorption spectroscopy, atomic force microscopy, and transmission electron microscopy, it is shown that the rtLL method advances the deposition of multi‐Langmuir layers and enables the production of films with defined morphology. The variation in nearest‐neighbor interactions is thus achieved and the resulting systematically tuned lowest unoccupied molecular orbital (LUMO) energies (determined via square‐wave voltammetry) enable the establishment of a model that functionally relates the LUMO energies to a morphological descriptor, allowing for the prediction of the range of accessible LUMO energies.

36 MATERIALS SCIENCE↗
DOE OSTI2018

In vivo cross-linking supports a head-to-tail mechanism for regulation of the plant plasma membrane P-type H + -ATPase

In higher plants, a P-type proton pumping ATPase generates the proton-motive force essential for the function of all other transporters and for proper growth and development. X-ray crystallographic studies of the plant plasma membrane proton pump have provided information on amino acids involved in ATP catalysis but provided no information on the structure of the C-terminal regulatory domain. Despite progress in elucidating enzymes involved in the signaling pathways that activate or inhibit this pump, the site of interaction of the C-terminal regulatory domain with the catalytic domains remains a mystery. Genetic studies have pointed to amino acids in various parts of the protein that may be involved but direct chemical evidence for which ones are specifically interacting with the C-terminus is lacking. In this study we used in vivo cross-linking experiments with a photo-reactive unnatural amino acid, p-benzoyl-phenylalanine (BPA) and tandem mass spectrometry, to obtain direct evidence that the carboxyl terminal regulatory domain directly interacts with amino acids located within the N-terminal actuator domain. Our observations are consistent with a mechanism in which intermolecular, rather than intramolecular, interactions are involved. Our model invokes a ‘head-to-tail’ organization of ATPase monomers in which the carboxyl terminal domain of one ATPase molecule interacts with the actuator domain of another ATPase molecule. This model serves to explain why cross-linked peptides are found only in dimers and trimers, and is consistent with prior studies suggesting that within the membrane, the protein can be organized as homopolymers, including dimers, trimers and hexamers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Character of intermolecular interaction in pyridine-argon complex: Ab initio potential energy surface, internal dynamics, and interrelations between SAPT energy components

The pyridine-Ar (PAr) van der Waals (vdW) complex is studied using a high level ab initio method. Its structure, binding energy, and intermolecular vibrational states are determined from the analytical potential energy surface constructed from interaction energy (IE) values computed at the coupled cluster level of theory with single, double, and perturbatively included triple excitations with the augmented correlation consistent polarized valence double-ζ (aug-cc-pVDZ) basis set complemented by midbond functions. The structure of the complex at its global minimum with Ar at a distance of 3.509 Å from the pyridine plane and shifted by 0.218 Å from the center of mass towards nitrogen agrees well with the corresponding equilibrium structure derived previously from the rotational spectrum of PAr. The PAr binding energy D{sub e} of 392 cm{sup −1} is close to that of 387 cm{sup −1} calculated earlier at the same ab initio level for the prototypical benzene-Ar (BAr) complex. However, under an extension of the basis set, D{sub e} for PAr becomes slightly lower than D{sub e} for BAr. The ab initio vdW vibrational energy levels allow us to estimate the reliability of the methods for the determination of the vdW fundamentals from the rotational spectra. To disclose the character of the intermolecular interaction in PAr, the symmetry-adapted perturbation theory (SAPT) is employed for the analysis of different physical contributions to IE. It is found that SAPT components of IE can be approximately expressed in the binding region by only two of them: the exchange repulsion and dispersion energy. The total induction effect is negligible. The interrelations between various SAPT components found for PAr are fulfilled for a few other complexes involving aromatic molecules and Ar or Ne, which indicates that they are valid for all rare gas (Rg) atoms and aromatics.

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

Dynamical screening of van der Waals interactions in nanostructured solids: Sublimation of fullerenes

Sublimation energy is one of the most important properties of molecular crystals, but it is difficult to study, because the attractive long-range van der Waals (vdW) interaction plays an important role. Here, we apply efficient semilocal density functional theory (DFT), corrected with the dynamically screened vdW interaction (DFT + vdW), the Rutgers-Chalmers nonlocal vdW-DF, and the pairwise-based dispersion-corrected DFT-D2 developed by Grimme and co-workers, to study the sublimation of fullerenes. We find that the short-range part, which accounts for the interaction due to the orbital overlap between fullerenes, is negligibly small. Our calculation shows that there exists a strong screening effect on the vdW interaction arising from the valence electrons of fullerenes. On the other hand, higher-order contributions can be as important as the leading-order term. The reasons are that (i) the surface of fullerene molecules is metallic and thus highly polarizable, (ii) the band gap of fullerene solids is small (less than 2 eV), and (iii) fullerene molecules in the solid phase are so densely packed, yielding the high valence electron density and small equilibrium intermolecular distances (the first nearest neighbor distance is only about 10 Å for C{sub 60}). However, these two effects make opposite contributions, leading to significant error cancellation between these two contributions. We demonstrate that, by considering higher-order contributions and the dynamical screening, the DFT + vdW method can yield sublimation energies of fullerenes in good agreement with reference values, followed by vdW-DF and DFT-D2. The insights from this study are important for a better understanding of the long-range nature of vdW interactions in nanostructured solids.

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

Theoretical studies for the N{sub 2}–N{sub 2}O van der Waals complex: The potential energy surface, intermolecular vibrations, and rotational transition frequencies

Theoretical studies of the potential energy surface (PES) and bound states are performed for the N{sub 2}–N{sub 2}O van der Waals (vdW) complex. A four-dimensional intermolecular PES is constructed at the level of single and double excitation coupled-cluster method with a non-iterative perturbation treatment of triple excitations [CCSD(T)] with aug-cc-pVTZ basis set supplemented with bond functions. Two equivalent T-shaped global minima are located, in which the O atom of N{sub 2}O monomer is near the N{sub 2} monomer. The intermolecular fundamental vibrational states are assigned by inspecting the orientation of the nodal surface of the wavefunctions. The calculated frequency for intermolecular disrotation mode is 23.086 cm{sup −1}, which is in good agreement with the available experimental data of 22.334 cm{sup −1}. A negligible tunneling splitting with the value of 4.2 MHz is determined for the ground vibrational state and the tunneling splitting increases as the increment of the vibrational frequencies. Rotational levels and transition frequencies are calculated for both isotopomers {sup 14}N{sub 2}–N{sub 2}O and {sup 15}N{sub 2}–N{sub 2}O. The accuracy of the PES is validated by the good agreement between theoretical and experimental results for the transition frequencies and spectroscopic parameters.

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

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

In this work, synchrotron X-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol and D-sorbitol) from above the glass transition temperatures $T_g$ to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently reflecting the different natures of bonding. Upon cooling from 400 K, the O$\cdots$O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrest at $T_g$; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above $T_g$, glycerol loses the fewest HBs per OH, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C$\cdots$C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above $T_g$ causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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

Cytochrome Nanowire XRD

Microbial metabolism plays an important role in the global cycling of carbon, nutrients, and metals. Recent discovery of microorganisms sharing energy using direct electrical connections might represent a prominent strategy by which anaerobic microorganisms interact with the diversity of environments. The purpose of the proposed research is to apply these newly developed principles of electron exchange in microbial communities to the study of biogeochemical cycling in climate and subsurface systems. Previous studies revealed that soil bacteria Geobacter sulfurreducens produce conductive pili nanofilaments that facilitate long-range electron transport to extracellular electron acceptors and other cells. This discovery has already transformed our understanding of the function of microbial communities in diverse environments, making basic contributions to microbial ecology and helping to improve practical applications such as increasing current output of microbial fuel cells and enhancing the conversion of organic waste to methane. Microbial electron exchange via pili may also contribute to methane production in terrestrial environments that are major sources of atmospheric methane. However, the mechanism of electron transport through pili is not well understood. Initial studies have suggested that pili exhibit metallic-like conductivity similar to conducting polymer polyaniline. However, these studies were limited to electrical measurements. X-ray diffraction studies have suggested that aromatic amino acids tyrosine and phenylalanine are closely packed in pili. But these studies used homology modeling to predict the pilus assembly and information about the number of monomers in each pilus assembly unit and their organization are still lacking. The central hypothesis of this proposal is aromatics in pili enable intermolecular electron delocalization due to pi-stacking that can give rise to metallic-like conductivity. By analyzing pili structure using infrared scattering-scanning near field optical microscope (IR s-SNOM) and combining with solution NMR, Focused Ion Beam (FIB and computational capabilities, this work will provide comprehensive understanding of the conduction mechanism from a structure-function perspective. The specific aims are: 1) With infrared scattering scanning near field optical microscope (IR s-SNOM), image the organization and density of pilin monomers and pH-induced conformational changes in a pilus filament. 2) With solution NMR, resolve pH-induced structural changes in the pilin monomer. 3) Using Focused Ion Beam (FIB), draw electrodes on single pili to measure their conductivity. 4) With computational modeling using the NWChem, simulate the effect of environment on the electronic properties of pili as a function of pH. The proposed work will gain the nanoscale insight into the structure of the pili, as the basis for understanding interspecies electrical communications. The improved understanding will be harnessed to accelerate bioenergy production and bioremediation using pili as well as for predictive modeling of carbon cycling in diverse environments to reach the central goals of DOE-BER to develop sustainable energy sources, regulate the contaminants in the subsurface and understand the effects of greenhouse gas emissions on biosphere. The proposed studies will provide design principles for engineering of electrical interactions among these communities by manipulating the amino acid composition of pili. Moreover, these studies will identify the forces that maintain the function of these communities in the face of environmental perturbation. Therefore, these studies will help to address the need "to incorporate process understanding of biogeochemical cycling into a scalable hierarchy of predictive capabilities" specified in EMSL's Science Themes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Structure, Hydrodynamics, and Phase Transition of Freely Suspended Liquid Crystals

Smectic liquid crystals are phases of rod shaped molecules organized into one dimensionally (1D) periodic arrays of layers, each layer being between one and two molecular lengths thick. In the least ordered smectic phases, the smectics A and C, each layer is a two dimensional (2D) liquid. Additionally there are a variety of more ordered smectic phases having hexatic short range translational order or 2D crystalline quasi long range translational order within the layers. The inherent fluid-layer structure and low vapor pressure of smectic liquid crystals enable the long term stabilization of freely suspended, single component, layered fluid films as thin as 30A, a single molecular layer. The layering forces the films to be an integral number of smectic layers thick, quantizing their thickness in layer units and forcing a film of a particular number of layers to be physically homogeneous with respect to its layer structure over its entire area. Optical reflectivity enables the precise determination of the number of layers. These ultrathin freely suspended liquid crystal films are structures of fundamental interest in condensed matter and fluid physics. They are the thinnest known stable condensed phase fluid structures and have the largest surface-to-volume ratio of any stable fluid preparation, making them ideal for the study of the effects of reduced dimensionality on phase behavior and on fluctuation and interface phenomena. Their low vapor pressure and quantized thickness enable the effective use of microgravity to extend the study of basic capillary phenomena to ultrathin fluid films. Freely suspended films have been a wellspring of new liquid crystal physics. They have been used to provide unique experimental conditions for the study of condensed phase transitions in two dimensions. They are the only system in which the hexatic has been unambiguously identified as a phase of matter, and the only physical system in which fluctuations of a 2D XY system and Kosterlitz Thouless phase transition has been observed and 2D XY quasi long range order verified. Smectic films have enabled the precise determination of smectic layer electron density and positional fluctuation profile and have been used to show that the interlayer interactions in anti-ferroelectric tilted smectics do not extend significantly beyond nearest neighbors. The interactions which are operative in liquid crystals are generally weak in comparison to those in crystalline phases, leading to the facile manipulation of the order in liquid crystals by external agents such as applied fields and surfaces. Effects arising from weak ordering are significantly enhanced in ultrathin free films and filaments wherein the intermolecular coupling is effectively reduced by loss of neighbors. Over the past four years this research, which we now detail, has produced a host of exciting new discoveries and unexpected results, maintaining the position of the study of freely suspended liquid crystal structures as one of most exciting and fruitful areas of complex fluid physics. In addition, several potentially interesting microgravity free film experiments have been identified.

Clark, Noel A.↗
NASA NTRS2004

Structure and Dynamics of Freely Suspended Liquid Crystals

Smectic liquid crystals are phases of rod shaped molecules organized into one dimensionally (1 D) periodic arrays of layers, each layer being between one and two molecular lengths thick. In the least ordered smectic phases, the smectics A and C, each layer is a two dimensional (2D) liquid. Additionally there are a variety of more ordered smectic phases having hexatic short range translational order or 2D crystalline or quasi long range translational order within the layers. The inherent fluid-layer structure and low vapor pressure of smectic liquid crystals enables the long term stabilization of freely suspended, single component, layered fluid films as thin as 30A, a single molecular layer. The layering forces the films to be an integral number of smectic layers thick, quantizing their thickness in layer units and forcing a film of a particular number of layers to be physically homogeneous with respect to its layer structure over its entire area. Optical reflectivity enables the precise determination of the number of layers. These ultrathin freely suspended liquid crystal films are structures of fundamental interest in condensed matter and fluid physics. They are the thinnest known stable fluid structures and have the largest surface-to-volume ratio of any stable fluid preparation, making them ideal for the study of the effects of reduced dimensionality on phase behavior and on fluctuation and interface phenomena. Their low vapor pressure and quantized thickness enable the effective use of microgravity to extend the study of basic capillary phenomena to ultrathin fluid films. Freely suspended films have been a wellspring of new LC physics. They have been used to provide unique experimental conditions for the study of condensed phase transitions in two dimensions. They are the only system in which the hexatic has been unambiguously identified as a phase of matter, and the only physical system in which fluctuations of a 2D XY system and Kosterlitz Thouless phase transition has been observed and 2D XY quasi long range order verified. Smectic films have enabled the precise determination of smectic layer electron density and positional fluctuation profiles and have been used to show that the interlayer interactions in antiferroelectric tilted smectics do not extend significantly beyond nearest neighbors. Freely suspended films played a pivotal role in the recent discovery of macroscopic chiral-polar ordering in fluids of achiral molecules. The interactions which are operative in liquid crystals are generally weak in comparison to those in crystalline phases, leading to the facile manipulation of the order in liquid crystals by external agents such as applied fields and surfaces. Effects arising from weak ordering are significantly enhanced in ultrathin free films and filaments, in which the intermolecular coupling is effectively further reduced by loss of neighbors. Over the past four years this research, which we now detail, has produced a host of exciting new discoveries and unexpected results, maintaining the study of freely suspended liquid crystal structures as one of most exciting and fruitful areas of complex fluid physics. In addition, a class of experiments on the behavior of 1D interfaces in 2D films have been pursued with results that point to potentially quite interesting effects in microgravity.

Clark, Noel A.↗
NASA NTRS1998

Modeling the Growth Rates of Tetragonal Lysozyme Crystal Faces

The measured macroscopic growth rates of the (110) and (101) faces of tetragonal lysozyme show an unexpectedly complex dependence on the supersaturation. The growth rates decay asymptotically to zero when the supersaturation is lowered to zero and increase rapidly when the supersaturation is increased. When supersaturations are increased still further the growth rates attain a maximum before starting to decrease. However, growth of these crystals is known to proceed by the classical dislocation and 2D nucleation growth mechanisms. This anomaly can be explained if growth is assumed to occur not by monomer units but by lysozyme aggregates. Analysis of the molecular packing of these crystals revealed that they were constructed of strongly bonded 4(sub 3) helices, while weaker bonds were responsible for binding the helices to each other. It follows that during crystal growth the stronger bonds are formed before the weaker ones. Thus, the growth of these crystals could be viewed as a two step process: aggregate growth units corresponding to the 4(sub 3) helix are first formed in the bulk solution by stronger intermolecular bonds and then attached to the crystal face by weaker bonds on dislocation hillocks or 2D islands. This will lead to a distribution of aggregates in the solution with monomers and lower order aggregates being predominant at low supersaturations and higher order aggregates being predominant at high supersaturations. If the crystal grows mostly by higher order aggregates, such as tetramers and octamers, it would explain the anomalous dependence of the growth rates on the supersaturation. Besides the analysis of molecular packing, a comprehensive analysis of the measured (110) and (101) growth rates was also undertaken in this study. The distribution of aggregates in lysozyme nutrient solutions at various solution conditions were determined from reversible aggregation reactions at equilibrium. The supersaturation was defined for each aggregate species with respect to its concentration at saturation in order to apply growth rate models to this process. The measured growth rates were then compared with the predicted ones from several dislocation and 2D nucleation growth models, employing tetramer and octamer growth units in polydisperse solutions and monomer units in monodisperse solutions. For the (110) face, the calculations consistently showed that the measured growth rates followed the expected model relations with octamer growth units. For the (101) face, it is not possible to obtain a clear agreement between the predicted and measured growth rates for a single growth unit as done for the (110) face. However, the calculations do indicate that the average size of the growth unit is between a tetramer and an octamer. This suggests that tetramers, octamers and other intermediate size growth units all participate in the growth process for this face. These calculations show that it is possible to model the macroscopic protein crystal growth rates if the molecular level processes can be account for, particularly protein aggregation processes in the bulk solution. Our recent investigations of tetragonal lysozyme crystals employing high resolution atomic force microscopy scans have further confirmed the growth of these crystals by aggregate growth units corresponding to 4(sub 3) helices.

Li, Meirong↗