Source author record

Christopher R. Iacovella

Christopher R. Iacovella appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2020arXiv

Towards Molecular Simulations that are Transparent, Reproducible, Usable By Others, and Extensible (TRUE)

Systems composed of soft matter (e.g., liquids, polymers, foams, gels, colloids, and most biological materials) are ubiquitous in science and engineering, but molecular simulations of such systems pose particular computational challenges, requiring time and/or ensemble-averaged data to be collected over long simulation trajectories for property evaluation. Performing a molecular simulation of a soft matter system involves multiple steps, which have traditionally been performed by researchers in a "bespoke" fashion, resulting in many published soft matter simulations not being reproducible based on the information provided in the publications. To address the issue of reproducibility and to provide tools for computational screening, we have been developing the open-source Molecular Simulation and Design Framework (MoSDeF) software suite. In this paper, we propose a set of principles to create Transparent, Reproducible, Usable by others, and Extensible (TRUE) molecular simulations. MoSDeF facilitates the publication and dissemination of TRUE simulations by automating many of the critical steps in molecular simulation, thus enhancing their reproducibility. We provide several examples of TRUE molecular simulations: All of the steps involved in creating, running and extracting properties from the simulations are distributed on open-source platforms (within MoSDeF and on GitHub), thus meeting the definition of TRUE simulations.

preprint2015arXiv

Development of a coarse-grained water forcefield via multistate iterative Boltzmann inversion

A coarse-grained water model is developed using multistate iterative Boltzmann inversion. Following previous work, the k-means algorithm is used to dynamically map multiple water molecules to a single coarse-grained bead, allowing the use of structure-based coarse-graining methods. The model is derived to match the bulk and interfacial properties of liquid water and improves upon previous work that used single state iterative Boltzmann inversion. The model accurately reproduces the density and structural correlations of water at 305 K and 1.0 atm, stability of a liquid droplet at 305 K, and shows little tendency to crystallize at physiological conditions. This work also illustrates several advantages of using multistate iterative Boltzmann inversion for deriving generally applicable coarse-grained forcefields.

preprint2014arXiv

Derivation of coarse-grained potentials via multistate iterative Boltzmann inversion

In this work, an extension to the standard iterative Boltzmann inversion (IBI) method to derive coarse-grained potentials is proposed. It is shown that the inclusion of target data from multiple states yields a less state-dependent potential, and is thus better suited to simulate systems over a range of thermodynamic states than the standard IBI method. The inclusion of target data from multiple states forces the algorithm to sample regions of potential phase space that match the radial distribution function at multiple state points, thus producing a derived potential that is more representative of the underlying potential interactions. It is shown that the algorithm is able to converge to the true potential for a system where the underlying potential is known. It is also shown that potentials derived via the proposed method better predict the behavior of n-alkane chains than those derived via the standard method. Additionally, through the examination of alkane monolayers, it is shown that the relative weight given to each state in the fitting procedure can impact bulk system properties, allowing the potentials to be further tuned in order to match the properties of reference atomistic and/or experimental systems.

preprint2013arXiv

Atomistic Simulations of Highly Conductive Molecular Transport Junctions Under Realistic Conditions

We report state-of-the-art atomistic simulations combined with high-fidelity conductance calculations to probe structure-conductance relationships in Au-benzenedithiolate(BDT)-Au junctions under elongation. Our results demonstrate that large increases in conductance are associated with the formation of monatomic chains (MACs) of Au atoms directly connected to BDT. An analysis of the electronic structure of the simulated junctions reveals that enhancement in the s-like states in Au MACs causes the increases in conductance. Other structures also result in increased conductance but are too short-lived to be detected in experiment, while MACs remain stable for long simulation times. Examinations of thermally evolved junctions with and without MACs show negligible overlap between conductance histograms, indicating that the increase in conductance is related to this unique structural change and not thermal fluctuation. These results, which provide an excellent explanation for a recently observed anomalous experimental result [Bruot et al., Nature Nanotech., 2012, 7, 35-40], should aid in the development of mechanically responsive molecular electronic devices.

preprint2013arXiv

Probing the Statistical Validity of the Ductile-to-Brittle Transition in Metallic Nanowires Using GPU Computing

We perform a large-scale statistical analysis (> 2000 independent simulations) of the elongation and rupture of gold nanowires, probing the validity and scope of the recently proposed ductile-to-brittle transition that occurs with increasing nanowire length [Wu et. al., Nano Lett., 12, 910-914 (2012)]. To facilitate a high-throughput simulation approach, we implement the second-moment approximation to the tight-binding (TB-SMA) potential within HOOMD-Blue, a molecular dynamics package which runs on massively parallel graphics processing units (GPUs). In a statistical sense, we find that the nanowires obey the ductile-to-brittle model quite well; however, we observe several unexpected features from the simulations that build on our understanding of the ductile-to-brittle transition. First, occasional failure behavior is observed that qualitatively differs from that predicted by the model prediction; this is attributed to stochastic thermal motion of the Au atoms and occurs at temperatures as low as 10 K. In addition, we also find that the ductile-to-brittle model, which was developed using classical dislocation theory, holds for nanowires as small as 3 nm in diameter. Finally, we demonstrate that the nanowire critical length is higher at 298 K relative to 10 K, a result that is not predicted by the ductile-to-brittle model. These results offer practical design strategies for adjusting nanowire failure and structure, and also demonstrate that GPU computing is an excellent tool for studies requiring a large number independent trajectories in order to fully characterize a system's behavior.

preprint2013arXiv

Structural Origins of Conductance Fluctuations in Gold-Thiolate Molecular Transport Junctions

We report detailed atomistic simulations combined with high-fidelity conductance calculations to probe the structural origins of conductance fluctuations in thermally evolving Au-benzene-1,4-dithiolate-Au junctions. We compare the behavior of structurally ideal junctions (electrodes with flat surfaces) to structurally realistic, experimentally representative junctions resulting from break junction simulations. The enhanced mobility of metal atoms in structurally realistic junctions results in significant changes to the magnitude and origin of the conductance fluctuations. Fluctuations are larger by a factor of 2-3 in realistic junctions compared to ideal junctions. Moreover, in junctions with highly deformed electrodes, the conductance fluctuations arise primarily from changes in the Au geometry, in contrast to results for junctions with non-deformed electrodes, where the conductance fluctuations are dominated by changes in the molecule geometry. These results provide important guidance to experimentalists developing strategies to control molecular conductance for device applications, and also to theoreticians invoking simplified structural models of junctions to predict their behavior.

preprint2012arXiv

Large-Scale Atomistic Simulations of Environmental Effects on the Formation and Properties of Molecular Junctions

Using an updated simulation tool, we examine molecular junctions comprised of benzene-1,4-dithiolate bonded between gold nanotips, focusing on the importance of environmental factors and inter-electrode distance on the formation and structure of bridged molecules. We investigate the complex relationship between monolayer density and tip separation, finding that the formation of multi-molecule junctions is favored at low monolayer density, while single-molecule junctions are favored at high density. We demonstrate that tip geometry and monolayer interactions, two factors that are often neglected in simulation, affect the bonding geometry and tilt angle of bridged molecules. We further show that the structures of bridged molecules at 298 and 77 K are similar.

preprint2011arXiv

Self Assembly of Soft Matter Quasicrystals and Their Approximants

The surprising recent discoveries of quasicrystals and their approximants in soft matter systems poses the intriguing possibility that these structures can be realized in a broad range of nano- and micro-scale assemblies. It has been theorized that soft matter quasicrystals and approximants are largely entropically stabilized, but the thermodynamic mechanism underlying their formation remains elusive. Here, we use computer simulation and free energy calculations to demonstrate a simple design heuristic for assembling quasicrystals and approximants in soft matter systems. Our study builds on previous simulation studies of the self-assembly of dodecagonal quasicrystals and approximants in minimal systems of spherical particles with complex, highly-specific interaction potentials. We demonstrate an alternative entropy-based approach for assembling dodecagonal quasicrystals and approximants based solely on particle functionalization and shape, thereby recasting the interaction-potential-based assembly strategy in terms of simpler-to-achieve bonded and excluded-volume interactions. Here, spherical building blocks are functionalized with mobile surface entities to encourage the formation of structures with low surface contact area, including non-close-packed and polytetrahedral structures. The building blocks also possess shape polydispersity, where a subset of the building blocks deviate from the ideal spherical shape, discouraging the formation of close-packed crystals. We show that three different model systems with both of these features -- mobile surface entities and shape polydispersity -- consistently assemble quasicrystals and/or approximants. We argue that this design strategy can be widely exploited to assemble quasicrystals and approximants on the nano- and micro- scales. In addition, our results further elucidate the formation of soft matter quasicrystals in experiment.

preprint2011arXiv

The Influence of Molecular Adsorption on Elongating Gold Nanowires

Using molecular dynamics simulations, we study the impact of physisorbing adsorbates on the structural and mechanical evolution of gold nanowires (AuNWs) undergoing elongation. We used various adsorbate models in our simulations, with each model giving rise to a different surface coverage and mobility of the adsorbed phase. We find that the local structure and mobility of the adsorbed phase remains relatively uniform across all segments of an elongating AuNW, except for the thinning region of the wire where the high mobility of Au atoms disrupts the monolayer structure, giving rise to higher solvent mobility. We analyzed the AuNW trajectories by measuring the ductile elongation of the wires and detecting the presence of characteristic structural motifs that appeared during elongation. Our findings indicate that adsorbates facilitate the formation of high-energy structural motifs and lead to significantly higher ductile elongations. In particular, our simulations result in a large number of monatomic chains and helical structures possessing mechanical stability in excess of what we observe in vacuum. Conversely, we find that a molecular species that interacts weakly (i.e., does not adsorb) with AuNWs worsens the mechanical stability of monatomic chains.

preprint2011arXiv

The Role of Polytetrahedral Structures in the Elongation and Rupture of Gold Nanowires

We report comprehensive high-accuracy molecular dynamics simulations using the ReaxFF forcefield to explore the structural changes that occur as Au nanowires are elongated, establishing trends as a function of both temperature and nanowire diameter. Our simulations and subsequent quantitative structural analysis reveal that polytetrahedral structures (e.g., icosahedra) form within the "amorphous" neck regions, most prominently for systems with small diameter at high temperature. We demonstrate that the formation of polytetrahedra diminishes the conductance quantization as compared to systems without this structural motif. We demonstrate that use of the ReaxFF forcefield, fitted to high-accuracy first principles calculations of Au, combines the accuracy of quantum calculations with the speed of semi-empirical methods.

preprint2010arXiv

Harmonic Order Parameters for Characterizing Complex Particle Morphologies

Order parameters based on spherical harmonics and Fourier coefficients already play a significant role in condensed matter research in the context of systems of spherical or point particles. Here, we extend these types of order parameter to more complex shapes, such as those encountered in nanoscale self-assembly applications. To do so, we build on a powerful set of techniques that originate in the computer science field of "shape matching." We demonstrate how shape matching techniques can be applied to identify unknown structures and create highly-specialized \textit{ad hoc} order parameters. Additionally, we investigate the special symmetry properties of harmonic descriptors, and demonstrate how they can be exploited to provide optimal solutions to certain classes of problems. Our techniques can be applied to particle systems in general, both simulated and experimental, provided the particle positions are known.

preprint2010arXiv

Pair Interaction Potentials of Colloids by Extrapolation of Confocal Microscopy Measurements of Collective Structure

A method for measuring the pair interaction potential between colloidal particles by extrapolation measurement of collective structure to infinite dilution is presented and explored using simulation and experiment. The method is particularly well suited to systems in which the colloid is fluorescent and refractive index matched with the solvent. The method involves characterizing the potential of mean force between colloidal particles in suspension by measurement of the radial distribution function using 3D direct visualization. The potentials of mean force are extrapolated to infinite dilution to yield an estimate of the pair interaction potential, $U(r)$. We use Monte Carlo (MC) simulation to test and establish our methodology as well as to explore the effects of polydispersity on the accuracy. We use poly-12-hydroxystearic acid-stabilized poly(methyl methacrylate) (PHSA-PMMA) particles dispersed in the solvent dioctyl phthalate (DOP) to test the method and assess its accuracy for three different repulsive systems for which the range has been manipulated by addition of electrolyte.