Researcher profile

Juan J. de Pablo

Juan J. de Pablo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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.