Source author record

Juan J. de Pablo

Juan J. de Pablo 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

23works
9topics
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

23 published item(s)

preprint2022arXiv

Current rectification by nanoparticles in bipolar nanopores

Bipolar nanochannels comprising two domains of positively and negatively charged walls along the pore axis are known to rectify current when exposed to an electric potential bias. We find that addition of charged nanoparticles can increase rectification considerably, by approximately one order of magnitude. Two bipolar channel geometries are considered here; their behavior is examined at rest and under the influence of a negative bias and a positive bias, respectively. We do so by relying on a molecular-level model of the electrolyte solution in the channels. The large increase in current rectification can be explained by the inherent electric field that charged nanoparticles generate within the channel. This effect is found to be largely dependent on the pore's geometry, its charge distribution, and the sign of the nanoparticles' charge, thereby offering new opportunities for design of engineered nanopore membrane-nanoparticle systems for energy storage.

preprint2022arXiv

Ionic transport in electrostatic Janus Membranes. An explicit solvent molecular dynamic simulation

Janus --or two-sided, charged membranes offer promise as ionic current rectifiers. In such systems, pores consisting of two regions of opposite charge can be used to generate a current from a gradient in salinity. The efficiency of \textcolor{black}{nanoscale} Janus pores increases dramatically as their diameter becomes smaller. However, little is known about the underlying transport processes\textcolor{black}{, particularly under experimentally accessible conditions}. In this work, we examine the molecular basis for rectification in Janus nanopores using an applied electric field. Molecular simulations with explicit water and ions are used to examine the structure and dynamics of all molecular species in aqueous electrolyte solutions. \textcolor{black}{For several macroscopic observables, the results of such simulations are consistent with experimental observations on asymmetric membranes. Our analysis reveals a number of previously unknown features, including a pronounced local re-orientation of water molecules in the pores, and a segregation of ionic species that has not been anticipated by previously reported continuum analyses of Janus pores. Using these insights, a model is proposed for ionic current rectification in which electric leakage at pore entrance controls net transport.

preprint2022arXiv

Sobolev Sampling of Free Energy Landscapes

A family of fast sampling methods is introduced here for molecular simulations of systems having rugged free energy landscapes. The methods represent a generalization of a strategy consisting of adjusting a model for the free energy as a function of one- or more collective variables as a simulation proceeds. Such a model is gradually built as a system evolves through phase space from both the frequency of visits to distinct states and generalized force estimates corresponding to such states. A common feature of the methods is that the underlying functional models and their gradients are easily expressed in terms of the same parameters, thereby providing faster and more effective fitting of the model from simulation data than other available sampling techniques. They also eliminate the need to train simultaneously separate neural networks, while retaining the advantage of generating smooth and continuous functional estimates that enable biasing outside the support grid. Implementation of the methods is relatively simple and, more importantly, they are found to provide gains of up to several orders of magnitude in computational efficiency over existing approaches.

preprint2022arXiv

Systematic Modification of Functionality in Disordered Elastic Networks Through Free Energy Surface Tailoring

Advances in manufacturing and characterization of complex molecular systems have created a need for new methods for design at molecular length scales. Emerging approaches are increasingly relying on the use of Artificial Intelligence (AI), and the training of AI models on large data libraries. This paradigm shift has led to successful applications, but shortcomings related to interpretability and generalizability continue to pose challenges. Here, we explore an alternative paradigm in which AI is combined with physics-based considerations for molecular and materials engineering. Specifically, collective variables, akin to those used in enhanced sampled simulations, are constructed using a machine learning (ML) model trained on data gathered from a single system. Through the ML-constructed collective variables, it becomes possible to identify critical molecular interactions in the system of interest, the modulation of which enables a systematic tailoring of the system's free energy landscape. To explore the efficacy of the proposed approach, we use it to engineer allosteric regulation, and uniaxial strain fluctuations in a complex disordered elastic network. Its successful application in these two cases provides insights regarding how functionality is governed in systems characterized by extensive connectivity, and points to its potential for design of complex molecular systems.

preprint2022arXiv

Transformation between elastic dipoles, quadrupoles, octupoles and hexadecapoles driven by surfactant self-assembly in nematic emulsion

Emulsions comprising isotropic fluid drops within a nematic host are of interest for applications ranging from biodetection to smart windows, which rely on changes of molecular alignment structures around the drops in response to chemical, thermal, electric and other stimuli. We show that absorption or desorption of trace amounts of common surfactants can drive continuous transformations of elastic multipoles induced by the droplets within the uniformly aligned nematic host. Out-of-equilibrium dynamics of director structures emerge from a controlled self-assembly or desorption of different surfactants at the drop-nematic interfaces, with ensuing forward and reverse transformations between elastic dipoles, quadrupoles, octupoles and hexadecapoles. We characterize inter-transformations of droplet-induced surface and bulk defects, probe elastic pair interactions and discuss emergent prospects for fundamental science and applications of the reconfigurable nematic emulsions.

preprint2021arXiv

Stability and molecular pathways to the formation of spin defects in silicon carbide

Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate--the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon ($V_{Si}$) to carbon monovacancies ($V_{C}$), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of $V_{Si}$ relative to $V_{C}$ favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.

preprint2021arXiv

Stress-activated Constraints in Dense Suspension Rheology

Dispersing small particles in a liquid can produce surprising behaviors when the solids fraction becomes large: rapid shearing drives these systems out of equilibrium and can lead to dramatic increases in viscosity (shear-thickening) or even solidification (shear jamming). These phenomena occur above a characteristic onset stress when particles are forced into frictional contact. Here we show via simulations how this can be understood within a framework that abstracts details of the forces acting at particle-particle contacts into general stress-activated constraints on relative particle movement. We find that focusing on just two constraints, affecting sliding and rolling at contact, can reproduce the experimentally observed shear thickening behavior quantitatively, despite widely different particle properties, surface chemistries, and suspending fluids. Within this framework parameters such as coefficients of sliding and rolling friction can each be viewed as a proxy for one or more forces of different physical or chemical origin, while the parameter magnitudes indicate the relative importance of the associated constraint. In this way, a new link is established that connects features observable in macroscale rheological measurements to classes of constraints arising from micro- or nano-scale properties.

preprint2020arXiv

Molecular Characterization of Ebselen Binding Activity to SARS-CoV-2 Main Protease

The Coronavirus Disease (COVID-19) pandemic caused by the SARS-coronavirus 2 (SARS-CoV-2) urgently calls for the design of drugs directed against this new virus. Given its essential role in proteolytic processing, the main protease Mpro has been identified as an attractive candidate for drugs against SARS-CoV-2 and similar coronaviruses. Recent high-throughput screening studies have identified a set of existing, small-molecule drugs as potent Mpro inhibitors. Amongst these, Ebselen (2-Phenyl-1,2-benzoselenazol-3-one), a glutathione peroxidase mimetic seleno-organic compound, is particularly attractive. Recent experiments suggest that its effectiveness is higher than that of other molecules that also act at the enzyme catalytic site. By relying on extensive simulations with all-atom models, in this study we examine at a molecular level the potential of Ebselen to decrease Mpro catalytic activity. Our results indicate that Ebselen exhibits a distinct affinity for the catalytic site cavity of Mpro. In addition, our molecular models reveal a second, previously unkown binding site for Ebselen in the dimerization region localized between the II and III domains of the protein. A detailed analysis of the free energy of binding indicates that the affinity of Ebselen to this second binding site is in fact significantly larger than that to the catalytic site. A strain analysis indicates that Ebselen bound between the II-III domains exerts a pronounced allosteric effect that regulates catalytic site access through surface loop interactions, and induces a displacement and reconfiguration of water hotspots, including the catalytic water, which could interfere with normal enzymatic function. Taken together, these findings provide a framework for the future design of more potent and specific Mpro inhibitors, based on the Ebselen scaffold, that could lead to new therapeutic strategies for COVID-19.

preprint2020arXiv

Sculpting stable structures in pure liquids

Pure liquids in thermodynamic equilibrium are structurally homogeneous. In liquid crystals, flow and light pulses are used to create reconfigurable domains with polar order. Moreover, through careful engineering of concerted microfluidic flows and localized opto-thermal fields, it is possible to achieve complete control over the nucleation, growth, and shape of such domains. Experiments, theory, and simulations indicate that the resulting structures can be stabilized indefinitely, provided the liquids are maintained in a controlled non-equilibrium state. The resulting sculpted liquids could find applications in microfluidic devices for selective encapsulation of solutes and particles into optically active compartments that interact with external stimuli.

preprint2020arXiv

Shape induced segregation and anomalous particle transport under spherical confinement

Colloid or nanoparticle mobility under confinement is of central importance to a wide range of physical and biological processes. Here, we introduce a minimal model of particles in a hydrodynamic continuum to examine how particle shape and concentration affect the transport of particles in spherical confinement. Specifically, an immersed boundary-General geometry Ewald-like approach is adopted to simulate the dynamics of spheres and cylinders under the influence of short-and long-range fluctuating hydrodynamic interactions with appropriate non-slip conditions at the confining walls. An efficient $\it{O(N)}$ parallel finite element algorithm is used, thereby allowing simulations at high concentrations, while a Chebyshev polynomial approximation is implemented in order to satisfy the fluctuation-dissipation theorem. A concentration-dependent anomalous diffusion is observed for suspended particles. It is found that introducing cylinders in a background of spheres, i.e. particles with a simple degree of anisotropy, has a pronounced influence on the structure and dynamics of the particles. First, increasing the fraction of cylinders induces a particle segregation effect, where spheres are pushed towards the wall and cylinders remain near the center of the cavity. This segregation leads to lower mobility for the spheres relative to that encountered in a system of pure spheres at the same volume fraction. Second, the diffusive-to-anomalous transition and the degree of anomaly--quantified by the power-law exponent in the mean square displacement vs. time relation-both increase as the fraction of cylinders becomes larger. These findings are of relevance for studies of diffusion in the cytoplasm, where proteins exhibit a distribution of size and shapes that could lead to some of the effects identified in the simulations reported here.

preprint2020arXiv

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: spheres and cylinders

The structure and dynamics of confined suspensions of particles of arbitrary shape is of interest in multiple disciplines, from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an Immersed-Boundary General geometry Ewald-like method to capture lubrication and long-range hydrodynamics, and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than $10\%$, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion towards the walls, and a sub-diffusive regime --caused by crowding-- in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.

preprint2016arXiv

Controlling mechanical properties of a polymeric glass

In this work we use molecular simulations to examine methods of controlling mechanical properties of polymeric glass materials such as elastic moduli, mechanical heterogeneity as well as their glass transition temperature. We study filled and unfilled polymers and examine the effect of particle size, volume fraction and polymer-particle interactions. We identify a relationship between mobility and dynamic heterogeneity with elastic moduli and glass transition temperature.

preprint2016arXiv

Endohedral confinement of a DNA dodecamer onto pristine carbon nanotubes and the stability of the canonical B form

Although carbon nanotubes are potential candidates for DNA encapsulation and subsequent delivery of biological payloads to living cells, the thermodynamical spontaneity of DNA encapsulation under physiological conditions is still a matter of debate. Using enhanced sampling techniques, we show for the first time that, given a sufficiently large carbon nanotube, the confinement of a double-stranded DNA segment, 5'-D(*CP*GP*CP*GP*AP*AP*TP*TP*CP*GP*CP*G)-3', is thermodynamically favourable under physiological environments (134 mM, 310 K, 1 bar), leading to DNA-nanotube hybrids with lower free energy than the unconfined biomolecule. A diameter threshold of 3 nm is established below which encapsulation is inhibited. The confined DNA segment maintains its translational mobility and exhibits the main geometrical features of the canonical B form. To accommodate itself within the nanopore, the DNA end-to-end length increases from 3.85 nm up to approximately 4.1 nm, due to a 0.3 nm elastic expansion of the strand termini. The canonical Watson-Crick H-bond network is essentially conserved throughout encapsulation, showing that the contact between the DNA segment and the hydrophobic carbon walls results in minor rearrangements of the nucleotides H-bonding. The results obtained here are paramount to the usage of carbon nanotubes as encapsulation media for next generation drug delivery technologies.

preprint2016arXiv

Inherent Structure Energy is a Good Indicator of Molecular Mobility in Glasses

Glasses produced via physical vapor deposition can display greater kinetic stability and lower enthalpy than glasses prepared by liquid cooling. While the reduced enthalpy has often been used as a measure of the stability, it is not obvious whether dynamic measures of stability provide the same view. Here, we study dynamics in vapor-deposited and liquid-cooled glass films using molecular simulations of a bead-spring polymer model as well as a Lennard-Jones binary mixture in two and three dimensions. We confirm that the dynamics in vapor-deposited glasses is indeed slower than in ordinary glasses. We further show that the inherent structure energy is a good reporter of local dynamics, and that aged systems and glasses prepared by cooling at progressively slower rates exhibit the same behavior as vapor-deposited materials when they both have the same inherent structure energy. These findings suggest that the stability inferred from measurements of the energy is also manifested in dynamic observables, and they strengthen the view that vapor deposition processes provide an effective strategy for creation of stable glasses.

preprint2016arXiv

Mesoscale Structure of Chiral Nematic Shells

There is considerable interest in understanding and controlling topological defects in nematic liquid crystals (LCs). Confinement, in the form of droplets, has been particularly effective in that regard. Here, we employ the Landau-de Gennes method to explore the geometrical frustration of nematic order in shell geometries, and focus on chiral materials. By varying the chirality and thickness in uniform shells, we construct a phase diagram that includes tetravalent structures, bipolar structures (BS), bent structures and radial spherical structures (RSS). It is found that, in uniform shells, the BS-to-RSS structural transition, in response to both chirality and shell geometry, is accompanied by an abrupt change of defect positions, implying a potential use for chiral nematic shells as sensors. Moreover, we investigate thickness heterogeneity in shells and demonstrate that non-chiral and chiral nematic shells exhibit distinct equilibrium positions of their inner core that are governed by shell chirality c.

preprint2016arXiv

Monte Carlo simulation study of diblock copolymer self assembly

A technique is presented which maps the parameters of a bead spring model, using the Flory Huggins theory, to a specific experimental system. By keeping only necessary details, for the description of these systems, the mapping procedure turns into an estimation of a few characteristic parameters. An asset of this technique is that it is simple to apply and captures the behavior of block copolymer phase separation. In our study this mapping technique is utilized in conjunction with a Monte Carlo (MC) algorithm to perform simulations on block copolymer systems. The microphase separation is investigated in the bulk and under confinement, on unpatterned and patterned surfaces.

preprint2016arXiv

Nanoscopic Characterization of DNA within Hydrophobic Pores:Thermodynamics and Kinetics

The energetic and transport properties of a double-stranded DNA dodecamer encapsulated in hydrophobic carbon nanotubes are probed employing two limiting nanotube diameters, D=4 nm and D=3 nm, corresponding to (51,0) and (40,0) zig-zag topologies, respectively. It is observed that the thermodynamically spontaneous encapsulation in the 4 nm nanopore (40 kJ/mol) is annihilated when the solid diameter narrows down to 3 nm, and that the confined DNA termini directly contact the hydrophobic walls with no solvent slab in-between. During the initial moments after confinement (2-3 ns),the biomolecule translocates along the nano pore's inner volume according to Fick's law (t) with a self-diffusion coefficient D=1.713 x 10-9m2/s, after which molecular diffusion assumes a single-file type mechanism (t1/2). As expected, diffusion is anisotropic, with the pore main axis as the preferred direction, but an in-depth analysis shows that the instantaneous velocity probabilities are essentially identical along the x, y and z directions. The 3D velocity histogram shows a maximum probability located at v=30.8 m/s, twice the observed velocity for a single-stranded three nucleotide DNA encapsulated in comparable armchair geometries (v=16.7 m/s, D=1.36-1.89 nm). Because precise physiological conditions (310 K, [NaCl] = 134 mM) are employed throughout, the present study establishes a landmark for the development of next generation in vivo drug delivery technologies based on carbon nanotubes as encapsulation agents.

preprint2016arXiv

Tension-dependent Free Energies of Nucleosome Unwrapping

Nucleosomes form the basic unit of compaction within eukaryotic genomes and their locations represent an important, yet poorly understood, mechanism of genetic regulation. Quantifying the strength of interactions within the nucleosome is a central problem in biophysics and is critical to understanding how nucleosome positions influence gene expression. By comparing to single-molecule experiments, we demonstrate that a coarse-grained molecular model of the nucleosome can reproduce key aspects of nucleosome unwrapping. Using detailed simulations of DNA and histone proteins, we calculate the tension-dependent free energy surface corresponding to the unwrapping process. The model reproduces quantitatively the forces required to unwrap the nucleosome, and reveals the role played by electrostatic interactions during this process. We then demonstrate that histone modifications and DNA sequence can have significant effects on the energies of nucleosome formation. Most notably, we show that histone tails are crucial for stabilizing the outer turn of nucleosomal DNA.

preprint2015arXiv

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged soft particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

preprint2015arXiv

Turning Statistical Physics Models Into Materials Design Engines

Despite the success statistical physics has enjoyed at predicting the properties of materials for given parameters, the inverse problem, identifying which material parameters produce given, desired properties, is only beginning to be addressed. Recently, several methods have emerged across disciplines that draw upon optimization and simulation to create computer programs that tailor material responses to specified behaviors. However, so far the methods developed either involve black-box techniques, in which the optimizer operates without explicit knowledge of the material's configuration space, or they require carefully tuned algorithms with applicability limited to a narrow subclass of materials. Here we introduce a formalism that can generate optimizers automatically by extending statistical mechanics into the realm of design. The strength of this new approach lies in its capability to transform statistical models that describe materials into optimizers to tailor them. By comparing against standard black-box optimization methods, we demonstrate how optimizers generated by this formalism can be faster and more effective, while remaining straightforward to implement. The scope of our approach includes new possibilities for solving a variety of complex optimization and design problems concerning materials both in and out of equilibrium.

preprint2014arXiv

Coarse-grained Modeling of DNA Curvature

Modeling of DNA-protein interactions is a complex process involving many important time and length scales. This can be facilitated through the use of coarse-grained models which reduce the number of degrees of freedom and allow efficient exploration of binding configurations. It is known that the local structure of DNA can significantly affect its protein-binding properties (i.e. intrinsic curvature in DNA-histone complexes). In a step towards comprehensive DNA-protein modeling, we expand the 3SPN.2 coarse-grained model to include intrinsic shape, and validate the refined model against experimental data including melting temperature, local flexibility, persistence length, and minor groove width profile.

preprint2013arXiv

Model Vapor-Deposited Glasses: Growth Front and Composition Effects

A growing body of experimental work indicates that physical vapor deposition provides an effective route for preparation of stable glasses, whose properties correspond in some cases to those expected for glasses that have been aged for thousands of years. In this work, model binary glasses are prepared in a process inspired by physical vapor deposition, in which particles are sequentially added to the free surface of a growing film. The resulting glasses are shown to be more stable than those prepared by gradual cooling from the liquid phase. However, it is also shown that the composition of the resulting glass, which is difficult to control in physical vapor deposition simulations of thin films, plays a significant role on the physical characteristics of the material. That composition dependence leads to a re- evaluation of previous results from simulations of thinner films than those considered here, where the equivalent age of the corresponding glasses was overestimated. The simulations presented in this work, which correspond to films that are approximately 38 molecular diameters thick, also enable analysis of the devitrification mechanism of vapor-deposited glasses. Consistent with experiments, it is found that this mechanism consists of a mobility front that propagates from the free interface into the interior of the films. Eliminating surface mobility eliminates this route of transformation into the supercooled liquid.

preprint2013arXiv

Optimizing directed self-assembled morphology

Directed assembly of block polymers is rapidly becoming a viable strategy for lithographic patterning of nanoscopic features. One of the key attributes of directed assembly is that an underlying chemical or topographic substrate pattern used to direct assembly need not exhibit a direct correspondence with the sought after block polymer morphology, and past work has largely relied on trial-and-error approaches to design appropriate patterns. In this work, a computational evolutionary strategy is proposed to solve this optimization problem. By combining the Cahn-Hilliard equation, which is used to find the equilibrium morphology, and the covariance-matrix evolutionary strategy, which is used to optimize the combined outcome of particular substrate-copolymer combinations, we arrive at an efficient method for design of substrates leading to non-trivial, desirable outcomes.