skip to main content

Title: Development of a PointNet for Detecting Morphologies of Self-Assembled Block Oligomers in Atomistic Simulations
ABSTRACT: Molecular simulations with atomistic or coarse- 6 grained force fields are a powerful approach for understanding and 7 predicting the self-assembly phase behavior of complex molecules. 8 Amphiphiles, block oligomers, and block polymers can form 9 mesophases with different ordered morphologies describing the 10 spatial distribution of the blocks, but entirely amorphous nature for 11 local packing and chain conformation. Screening block oligomer 12 chemistry and architecture through molecular simulations to find 13 promising candidates for functional materials is aided by effective 14 and straightforward morphology identification techniques. Captur- 15 ing 3-dimensional periodic structures, such as ordered network 16 morphologies, is hampered by the requirement that the number of 17 molecules in the simulated system and the shape of the periodic simulation box need to be commensurate with those of the resulting 18 network phase. Common strategies for structure identification include structure factors and order parameters, but these fail to 19 identify imperfect structures in simulations with incorrect system sizes. Building upon pioneering work by DeFever et al. [Chem. Sci. 20 2019, 10, 7503−7515] who implemented a PointNet (i.e., a neural network designed for computer vision applications using point 21 clouds) to detect local structure in simulations of single-bead particles and water molecules, we present a PointNet for detection of 22 nonlocal ordered morphologies of complex more » block oligomers. Our PointNet was trained using atomic coordinates from molecular 23 dynamics simulation trajectories and synthetic point clouds for ordered network morphologies that were absent from previous 24 simulations. In contrast to prior work on simple molecules, we observe that large point clouds with 1000 or more points are needed 25 for the more complex block oligomers. The trained PointNet model achieves an accuracy as high as 0.99 for globally ordered 26 morphologies formed by linear diblock, linear triblock, and 3-arm and 4-arm star-block oligomers, and it also allows for the discovery 27 of emerging ordered patterns from nonequilibrium systems. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
The journal of physical chemistry
Sponsoring Org:
National Science Foundation
More Like this
  1. Identifying local structure in molecular simulations is of utmost importance. The most common existing approach to identify local structure is to calculate some geometrical quantity referred to as an order parameter. In simple cases order parameters are physically intuitive and trivial to develop ( e.g. , ion-pair distance), however in most cases, order parameter development becomes a much more difficult endeavor ( e.g. , crystal structure identification). Using ideas from computer vision, we adapt a specific type of neural network called a PointNet to identify local structural environments in molecular simulations. A primary challenge in applying machine learning techniques to simulation is selecting the appropriate input features. This challenge is system-specific and requires significant human input and intuition. In contrast, our approach is a generic framework that requires no system-specific feature engineering and operates on the raw output of the simulations, i.e. , atomic positions. We demonstrate the method on crystal structure identification in Lennard-Jones (four different phases), water (eight different phases), and mesophase (six different phases) systems. The method achieves as high as 99.5% accuracy in crystal structure identification. The method is applicable to heterogeneous nucleation and it can even predict the crystal phases of atoms near external interfaces.more »We demonstrate the versatility of our approach by using our method to identify surface hydrophobicity based solely upon positions and orientations of surrounding water molecules. Our results suggest the approach will be broadly applicable to many types of local structure in simulations.« less
  2. Self-assembling dendrimers have facilitated the discovery of periodic and quasiperiodic arrays of supramolecular architectures and the diverse functions derived from them. Examples are liquid quasicrystals and their approximants plus helical columns and spheres, including some that disregard chirality. The same periodic and quasiperiodic arrays were subsequently found in block copolymers, surfactants, lipids, glycolipids, and other complex molecules. Here we report the discovery of lamellar and hexagonal periodic arrays on the surface of vesicles generated from sequence-defined bicomponent monodisperse oligomers containing lipid and glycolipid mimics. These vesicles, known as glycodendrimersomes, act as cell-membrane mimics with hierarchical morphologies resembling bicomponent rafts. These nanosegregated morphologies diminish sugar–sugar interactions enabling stronger binding to sugar-binding proteins than densely packed arrangements of sugars. Importantly, this provides a mechanism to encode the reactivity of sugars via their interaction with sugar-binding proteins. The observed sugar phase-separated hierarchical arrays with lamellar and hexagonal morphologies that encode biological recognition are among the most complex architectures yet discovered in soft matter. The enhanced reactivity of the sugar displays likely has applications in material science and nanomedicine, with potential to evolve into related technologies.
  3. Relationships between average network productivity and accumulation or density aggregated 2 across spatially compact regions of urban networks—so called network Macroscopic Fundamental 3 Diagrams (MFDs)—have recently been shown to exist. Various analytical methods have been put 4 forward to estimate a network’s MFD as a function of network properties, such as average block 5 lengths, signal timings, and traffic flow characteristics on links. However, real street networks are 6 not homogeneous—they generally have a hierarchical structure where some streets (e.g., arterials) 7 promote higher mobility than others (e.g., local roads). This paper provides an analytical method 8 to estimate the MFDs of hierarchical street networks by considering features that are specific to 9 hierarchical network structures. Since the performance of hierarchical networks is driven by how 10 vehicles are routed across the different street types, two routing conditions— user equilibrium and 11 system optimal routing—are considered in the analytical model. The proposed method is first 12 implemented to describe the MFD of a hierarchical one-way limited access linear corridor and 13 then extended to a more realistic hierarchical two-dimensional grid network. For both cases, it is 14 shown that the MFD of a hierarchical network may no longer be unimodalmore »or concave as 15 traditionally assumed in most MFD-based modeling frameworks. These findings are verified using 16 simulations of hierarchical corridors. Finally, the proposed methodology is applied to demonstrate 17 how it can be used to make decisions related to the design of hierarchical street network structures.« less
  4. Context. Formaldehyde is a potential biogenic precursor involved in prebiotic chemical evolution. The cold conditions of the interstellar medium (ISM) allow H 2 CO to be reactive, playing a significant role as a chemical intermediate in formation pathways leading to interstellar complex organic molecules. However, gas-phase molecular formation mechanisms in cold regions of the ISM are poorly understood. Aims. We computationally determine the most favored gas-phase molecular formation mechanisms at local thermodynamic equilibrium conditions that can produce the detected amounts of H 2 CO in diffuse molecular clouds (DMCs), in dark, cold, and dense molecular clouds (DCDMCs), and in three regions of circumstellar envelopes of low-mass protostars (CELMPs). Methods. The potential energy surfaces, thermodynamic functions, and single-point energies for transition states were calculated at the CCSD(T)-F12/cc-pVTZ-F12 and MP2/aug-cc-pVDZ levels of theory and basis sets. Molecular thermodynamics and related partition functions were obtained by applying the Maxwell-Boltzmann quantum statistics theory from energies computed at CCSD(T)-F12/cc-pVTZ-F12 with corrections for zero-point energy. A literature review on detected abundances of reactants helped us to propose the most favorable formation routes. Results. The most probable reactions that produce H 2 CO in cold astrophysical regions are: 1 CH 2 + ⋅ 3 O 2 →more »1 H 2 CO + O⋅( 3 P) in DMCs, ⋅ 3 CH 2 + ⋅ 3 O 2 → 1 H 2 CO + ⋅O( 3 P) in DCDMCs, and ⋅CH 3 + ⋅O( 3 P) → 1 H 2 CO + ⋅H in region III, ⋅CH 3 +⋅O( 1 D) → 1 H 2 CO + ⋅H in region II, and 1 CH 2 + ⋅ 3 O 2 → 1 H 2 CO + ⋅O( 3 P) in region I belonging to CELMPs. Conclusions. Quantum chemical calculations suggest that the principal carbonaceous precursors of H 2 CO in cold regions for the gas-phase are CH 2 (a 1 A 1 ), and ⋅CH 2 (X 3 B 1 ) combined with ⋅O 2 ( 3 Σ g ) and ⋅CH 3 ( 2 A ” ) + ⋅O( 3 P) / O( 1 D). Reactions based on more complex reagents yield less effective thermodynamics in the gas-phase H 2 CO molecular formation.« less
  5. Abstract. A tethered-balloon system (TBS) has been developed and is beingoperated by Sandia National Laboratories (SNL) on behalf of the U.S.Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) UserFacility in order to collect in situ atmospheric measurements withinmixed-phase Arctic clouds. Periodic tethered-balloon flights have beenconducted since 2015 within restricted airspace at ARM's Advanced MobileFacility 3 (AMF3) in Oliktok Point, Alaska, as part of the AALCO (AerialAssessment of Liquid in Clouds at Oliktok), ERASMUS (Evaluation of RoutineAtmospheric Sounding Measurements using Unmanned Systems), and POPEYE(Profiling at Oliktok Point to Enhance YOPP Experiments) field campaigns. Thetethered-balloon system uses helium-filled 34 m3 helikites and 79 and104 m3 aerostats to suspend instrumentation that is used to measureaerosol particle size distributions, temperature, horizontal wind, pressure,relative humidity, turbulence, and cloud particle properties and tocalibrate ground-based remote sensing instruments. Supercooled liquid water content (SLWC) sondes using the vibrating-wireprinciple, developed by Anasphere Inc., were operated at Oliktok Point atmultiple altitudes on the TBS within mixed-phase clouds for over 200 h.Sonde-collected SLWC data were compared with liquid water content derivedfrom a microwave radiometer, Ka-band ARM zenith radar, and ceilometer at the AMF3, as well as liquid water content derived from AMF3 radiosonde flights. The in situ data collected by the Anasphere sensors were also comparedmore »with data collected simultaneously by an alternative SLWC sensor developed at the University of Reading, UK; both vibrating-wire instruments were typically observed to shed their ice quickly upon exiting the cloud or reaching maximum ice loading. Temperature sensing measurements distributed with fiber optic tethered balloons were also compared with AMF3 radiosonde temperature measurements. Combined, the results indicate that TBS-distributedtemperature sensing and supercooled liquid water measurements are inreasonably good agreement with remote sensing and radiosonde-basedmeasurements of both properties. From these measurements and sensorevaluations, tethered-balloon flights are shown to offer an effective methodof collecting data to inform and constrain numerical models, calibrate andvalidate remote sensing instruments, and characterize the flight environmentof unmanned aircraft, circumventing the difficulties of in-cloud unmanned aircraft flights such as limited flight time and in-flight icing.« less