From nematic shells to nematic droplets: energetics and defect transitions.

In this work, we investigate the possibility of inducing valence transitions, i.e. transitions between different defect configurations, by transforming a nematic shell into a nematic droplet. Our shells are liquid crystal droplets containing a smaller aqueous droplet inside, which are suspended in an aqueous phase. When osmotically de-swelling the inner droplet, the shell progressively increases its thickness until it eventually becomes a single droplet. During the process, the shell energy landscape evolves, triggering a response in the system. We observe two different scenarios. Either the inner droplet progressively shrinks and disappears, inducing a defect reorganization, or it is expelled from the shell at a critical radius of the inner droplet, abruptly changing the geometry of the system. We use numerical simulations and modeling to investigate the origin of these behaviors. We find that the selected route depends on the defect structure and the energetics of the system as it evolves. The critical inner radius and time for expulsion depend on the osmotic pressure of the outer phase, suggesting that the flow through the shell plays a role in the process.

Amphiphile-Induced Phase Transition of Liquid Crystals at Aqueous Interfaces.

Monolayer assemblies of amphiphiles at planar interfaces between thermotropic liquid crystals (LCs) and an aqueous phase can give rise to configurational transitions of the underlying LCs. A common assumption has been that a reconfiguration of the LC phase is caused by an interdigitation of the hydrophobic tails of amphiphiles with the molecules of the LC at the interface. A different mechanism is discovered here, whereby reorientation of the LC systems is shown to occur through lowering of the orientation-dependent surface energy of the LC due to formation of a thin isotropic layer at the aqueous interface. Using a combination of atomistic molecular dynamics simulations and experiments, we demonstrate that a monolayer of specific amphiphiles at an aqueous interface can cause a local nematic-to-isotropic phase transition of the LC by disturbing the antiparallel configuration of the LC molecules. These results provide new insights into the interfacial, molecular-level organization of LCs that can be exploited for rational design of biological sensors and responsive systems.

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard-Jones potential.

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard-Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.

SSAGES: Software Suite for Advanced General Ensemble Simulations.

Molecular simulation has emerged as an essential tool for modern-day research, but obtaining proper results and making reliable conclusions from simulations requires adequate sampling of the system under consideration. To this end, a variety of methods exist in the literature that can enhance sampling considerably, and increasingly sophisticated, effective algorithms continue to be developed at a rapid pace. Implementation of these techniques, however, can be challenging for experts and non-experts alike. There is a clear need for software that provides rapid, reliable, and easy access to a wide range of advanced sampling methods and that facilitates implementation of new techniques as they emerge. Here we present SSAGES, a publicly available Software Suite for Advanced General Ensemble Simulations designed to interface with multiple widely used molecular dynamics simulations packages. SSAGES allows facile application of a variety of enhanced sampling techniques-including adaptive biasing force, string methods, and forward flux sampling-that extract meaningful free energy and transition path data from all-atom and coarse-grained simulations. A noteworthy feature of SSAGES is a user-friendly framework that facilitates further development and implementation of new methods and collective variables. In this work, the use of SSAGES is illustrated in the context of simple representative applications involving distinct methods and different collective variables that are available in the current release of the suite. The code may be found at: https://github.com/MICCoM/SSAGES-public.

Nature of peptide wrapping onto metal nanoparticle catalysts and driving forces for size control.

Colloidal metal nanocrystals find many applications in catalysis, energy conversion devices, and therapeutics. However, the nature of ligand interactions and implications on shape control have remained uncertain at the atomic scale. Large differences in peptide adsorption strength and facet specificity were found on flat palladium surfaces versus surfaces of nanoparticles of 2 to 3 nm size using accurate atomistic simulations with the Interface force field. Folding of longer peptides across many facets explains the formation of near-spherical particles with local surface disorder, in contrast to the possibility of nanostructures of higher symmetry with shorter ligands. The average particle size in TEM correlates inversely with the surface coverage with a given ligand and with the strength of ligand adsorption. The role of specific amino acids and sequence mutations on the nanoparticle size and facet composition is discussed, as well as the origin of local surface disorder that leads to large differences in catalytic reactivity.

Segregation of liquid crystal mixtures in topological defects.

The structure and physical properties of liquid crystal (LC) mixtures are a function of composition, and small changes can have pronounced effects on observables, such as phase-transition temperatures. Traditionally, LC mixtures have been assumed to be compositionally homogenous. The results of chemically detailed simulations presented here show that this is not the case; pronounced deviations of the local order from that observed in the bulk at defects and interfaces lead to significant compositional segregation effects. More specifically, two disclination lines are stabilized in this work by introducing into a nematic liquid crystal mixture a cylindrical body that exhibits perpendicular anchoring. It is found that the local composition deviates considerably from that of the bulk at the interface with the cylinder and in the defects, thereby suggesting new assembly and synthetic strategies that may capitalize on the unusual molecular environment provided by liquid crystal mixtures.

Molecular Structure of Canonical Liquid Crystal Interfaces.

Numerous applications of liquid crystals rely on control of molecular orientation at an interface. However, little is known about the precise molecular structure of such interfaces. In this work, synchrotron X-ray reflectivity measurements, accompanied by large-scale atomistic molecular dynamics simulations, are used for the first time to reconstruct the air-liquid crystal interface of a nematic material, namely, 4-pentyl-4'-cyanobiphenyl (5CB). The results are compared to those for 4-octyl-4'-cyanobiphenyl (8CB) which, in addition to adopting isotropic and nematic states, can also form a smectic phase. Our findings indicate that the air interface imprints a highly ordered structure into the material; such a local structure then propagates well into the bulk of the liquid crystal, particularly for nematic and smectic phases.

Water Flux Induced Reorientation of Liquid Crystals.

It is well understood that the adsorption of solutes at the interface between a bulk liquid crystal phase and an aqueous phase can lead to orientational or anchoring transitions. A different principle is introduced here, whereby a transient reorientation of a thermotropic liquid crystal is triggered by a spontaneous flux of water across the interface. A critical water flux can be generated by the addition of an electrolyte to the bulk aqueous phase, leading to a change in the solvent activity; water is then transported through the liquid crystal phase and across the interface. The magnitude of the spontaneous water flux can be controlled by the concentration and type of solutes, as well as the rate of salt addition. These results present new, previously unappreciated fundamental principles that could potentially be used for the design of materials involving transient gating mechanisms, including biological sensors, drug delivery systems, separation media, and molecular machines.

Understanding Atomic-Scale Behavior of Liquid Crystals at Aqueous Interfaces.

The ordered environment presented by liquid crystals at interfaces enables a range of novel functionalities that is only now beginning to be exploited in applications ranging from light focusing devices to biosensors. One key feature of liquid crystals is that molecular events occurring at an interface propagate over large distances through the bulk. In spite of their importance, our fundamental understanding of liquid crystal-water and liquid crystal-air interfaces remains limited. In this work, we present results from large-scale atomistic molecular dynamics simulations on the organization of the nematic and isotropic phases of the nitrile-containing mesogenic molecule 4-cyano-4'-pentylbiphenyl (5CB) in the vicinity of vacuum and aqueous interfaces. Hybrid boundary conditions are imposed by confining 5CB films between vacuum and an aqueous medium to examine how those two types of interfaces influence the specific structural arrangement and ordering of 5CB. Consistent with experiments, our results indicate that 5CB exhibits homeotropic anchoring at the vacuum interface, and planar alignment at aqueous interfaces. Two-dimensional molecular dynamics potential of mean force calculations and average polarization densities show that the polar nitrile group of 5CB remains hydrated near the aqueous interface, where it modulates the orientation of water molecules. Estimates of the anchoring strength reveal an oscillatory decay and a semilinear decay with distance from the interface in vacuum and water, respectively.

Simulations of inorganic-bioorganic interfaces to discover new materials: insights, comparisons to experiment, challenges, and opportunities.

Natural and man-made materials often rely on functional interfaces between inorganic and organic compounds. Examples include skeletal tissues and biominerals, drug delivery systems, catalysts, sensors, separation media, energy conversion devices, and polymer nanocomposites. Current laboratory techniques are limited to monitor and manipulate assembly on the 1 to 100 nm scale, time-consuming, and costly. Computational methods have become increasingly reliable to understand materials assembly and performance. This review explores the merit of simulations in comparison to experiment at the 1 to 100 nm scale, including connections to smaller length scales of quantum mechanics and larger length scales of coarse-grain models. First, current simulation methods, advances in the understanding of chemical bonding, in the development of force fields, and in the development of chemically realistic models are described. Then, the recognition mechanisms of biomolecules on nanostructured metals, semimetals, oxides, phosphates, carbonates, sulfides, and other inorganic materials are explained, including extensive comparisons between modeling and laboratory measurements. Depending on the substrate, the role of soft epitaxial binding mechanisms, ion pairing, hydrogen bonds, hydrophobic interactions, and conformation effects is described. Applications of the knowledge from simulation to predict binding of ligands and drug molecules to the inorganic surfaces, crystal growth and shape development, catalyst performance, as well as electrical properties at interfaces are examined. The quality of estimates from molecular dynamics and Monte Carlo simulations is validated in comparison to measurements and design rules described where available. The review further describes applications of simulation methods to polymer composite materials, surface modification of nanofillers, and interfacial interactions in building materials. The complexity of functional multiphase materials creates opportunities to further develop accurate force fields, including reactive force fields, and chemically realistic surface models, to enable materials discovery at a million times lower computational cost compared to quantum mechanical methods. The impact of modeling and simulation could further be increased by the advancement of a uniform simulation platform for organic and inorganic compounds across the periodic table and new simulation methods to evaluate system performance in silico.

Three-dimensional coordinates of individual atoms in materials revealed by electron tomography.

Crystallography, the primary method for determining the 3D atomic positions in crystals, has been fundamental to the development of many fields of science. However, the atomic positions obtained from crystallography represent a global average of many unit cells in a crystal. Here, we report, for the first time, the determination of the 3D coordinates of thousands of individual atoms and a point defect in a material by electron tomography with a precision of ∼19 pm, where the crystallinity of the material is not assumed. From the coordinates of these individual atoms, we measure the atomic displacement field and the full strain tensor with a 3D resolution of ∼1 nm(3) and a precision of ∼10(-3), which are further verified by density functional theory calculations and molecular dynamics simulations. The ability to precisely localize the 3D coordinates of individual atoms in materials without assuming crystallinity is expected to find important applications in materials science, nanoscience, physics, chemistry and biology.

Elucidation of peptide-directed palladium surface structure for biologically tunable nanocatalysts.

Peptide-enabled synthesis of inorganic nanostructures represents an avenue to access catalytic materials with tunable and optimized properties. This is achieved via peptide complexity and programmability that is missing in traditional ligands for catalytic nanomaterials. Unfortunately, there is limited information available to correlate peptide sequence to particle structure and catalytic activity to date. As such, the application of peptide-enabled nanocatalysts remains limited to trial and error approaches. In this paper, a hybrid experimental and computational approach is introduced to systematically elucidate biomolecule-dependent structure/function relationships for peptide-capped Pd nanocatalysts. Synchrotron X-ray techniques were used to uncover substantial particle surface structural disorder, which was dependent upon the amino acid sequence of the peptide capping ligand. Nanocatalyst configurations were then determined directly from experimental data using reverse Monte Carlo methods and further refined using molecular dynamics simulation, obtaining thermodynamically stable peptide-Pd nanoparticle configurations. Sequence-dependent catalytic property differences for C-C coupling and olefin hydrogenation were then elucidated by identification of the catalytic active sites at the atomic level and quantitative prediction of relative reaction rates. This hybrid methodology provides a clear route to determine peptide-dependent structure/function relationships, enabling the generation of guidelines for catalyst design through rational tailoring of peptide sequences.

Atomic-scale identification of Pd leaching in nanoparticle catalyzed C-C coupling: effects of particle surface disorder.

C-C coupling reactions are of great importance in the synthesis of numerous organic compounds, where Pd nanoparticle catalyzed systems represent new materials to efficiently drive these reactions. Despite their pervasive utility, the catalytic mechanism of these particle-based reactions remains highly contested. Herein we present evidence of an atom leaching mechanism for Stille coupling under aqueous conditions using peptide-capped Pd nanoparticles. EXAFS analysis revealed Pd coordination changes in the nanoparticle consistent with Pd atom abstraction, where sizing analysis by SAXS confirmed particle size changes associated with a leaching process. It is likely that recently discovered highly disordered surface Pd atoms are the favored catalytic active sites and are leached during oxidative addition, resulting in smaller particles. Probing the mechanism of nanoparticle-driven C-C coupling reactions through structural analyses provides fundamental information concerning these active sites and their reactivity at the atomic-scale, which can be used to improve catalytic performance to meet important sustainability goals.

A rational biomimetic approach to structure defect generation in colloidal nanocrystals.

Controlling the morphology of nanocrystals (NCs) is of paramount importance for both fundamental studies and practical applications. The morphology of NCs is determined by the seed structure and the following facet growth. While means for directing facet formation in NC growth have been extensively studied, rational strategies for the production of NCs bearing structure defects in seeds have been much less explored. Here, we report mechanistic investigations of high density twin formation induced by specific peptides in platinum (Pt) NC growth, on the basis of which we derive principles that can serve as guidelines for the rational design of molecular surfactants to introduce high yield twinning in noble metal NC syntheses. Two synergistic factors are identified in producing twinned Pt NCs with the peptide: (1) the altered reduction kinetics and crystal growth pathway as a result of the complex formation between the histidine residue on the peptide and Pt ions, and (2) the preferential stabilization of {111} planes upon the formation of twinned seeds. We further apply the discovered principles to the design of small organic molecules bearing similar binding motifs as ligands/surfactants to create single and multiple twinned Pd and Rh NCs. Our studies demonstrate the rich information derived from biomimetic synthesis and the broad applicability of biomimetic principles to NC synthesis for diverse property tailoring.

Exploiting localized surface binding effects to enhance the catalytic reactivity of peptide-capped nanoparticles.

Peptide-based methods represent new approaches to selectively produce nanostructures with potentially important functionality. Unfortunately, biocombinatorial methods can only select peptides with target affinity and not for the properties of the final material. In this work, we present evidence to demonstrate that materials-directing peptides can be controllably modified to substantially enhance particle functionality without significantly altering nanostructural morphology. To this end, modification of selected residues to vary the site-specific binding strength and biological recognition can be employed to increase the catalytic efficiency of peptide-capped Pd nanoparticles. These results represent a step toward the de novo design of materials-directing peptides that control nanoparticle structure/function relationships.

Stability, surface features, and atom leaching of palladium nanoparticles: toward prediction of catalytic functionality.

Surfactant-stabilized metal nanoparticles have shown promise as catalysts although specific surface features and their influence on catalytic performance have not been well understood. We quantify the thermodynamic stability, the facet composition of the surface, and distinct atom types that affect rates of atom leaching for a series of twenty near-spherical Pd nanoparticles of 1.8 to 3.1 nm size using computational models. Cohesive energies indicate higher stability of certain particles that feature an approximate 60/20/20 ratio of {111}, {100}, and {110} facets while less stable particles exhibit widely variable facet composition. Unique patterns of atom types on the surface cause apparent differences in binding energies and changes in reactivity. Estimates of the relative rate of atom leaching as a function of particle size were obtained by the summation of Boltzmann-weighted binding energies over all surface atoms. Computed leaching rates are in good qualitative correlation with the measured catalytic activity of peptide-stabilized Pd nanoparticles of the same shape and size in Stille coupling reactions. The agreement supports rate-controlling contributions by atom leaching in the presence of reactive substrates. The computational approach provides a pathway to estimate the catalytic activity of metal nanostructures of engineered shape and size, and possible further refinements are described.

Tailoring molecular specificity toward a crystal facet: a lesson from biorecognition toward Pt{111}.

Surfactants with preferential adsorption to certain crystal facets have been widely employed to manipulate morphologies of colloidal nanocrystals, while mechanisms regarding the origin of facet selectivity remain an enigma. Similar questions exist in biomimetic syntheses concerning biomolecular recognition to materials and crystal surfaces. Here we present mechanistic studies on the molecular origin of the recognition toward platinum {111} facet. By manipulating the conformations and chemical compositions of a platinum {111} facet specific peptide, phenylalanine is identified as the dominant motif to differentiate {111} from other facets. The discovered recognition motif is extended to convert nonspecific peptides into {111} specific peptides. Further extension of this mechanism allows the rational design of small organic molecules that demonstrate preferential adsorption to the {111} facets of both platinum and rhodium nanocrystals. This work represents an advance in understanding the organic-inorganic interfacial interactions in colloidal systems and paves the way to rational and predictable nanostructure modulations for many applications.