Researcher profile

A. Carpio

A. Carpio contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
16topics
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

15 published item(s)

preprint2023arXiv

Scale free chaos in the confined Vicsek flocking model

The Vicsek model encompasses the paradigm of active dry matter. Motivated by collective behavior of insects in swarms, we have studied finite size effects and criticality in the three dimensional, harmonically confined Vicsek model. We have discovered a phase transition that exists for appropriate noise and small confinement strength. On the critical line of confinement versus noise, swarms are in a state of scale-free chaos characterized by minimal correlation time, correlation length proportional to swarm size and topological data analysis. The critical line separates dispersed single clusters from confined multicluster swarms. Scale-free chaotic swarms occupy a compact region of space and comprise a recognizable `condensed' nucleus and particles leaving and entering it. Susceptibility, correlation length, dynamic correlation function and largest Lyapunov exponent obey power laws. The critical line and a narrow criticality region close to it move simultaneously to zero confinement strength for infinitely many particles. At the end of the first chaotic window of confinement, there is another phase transition to infinitely dense clusters of finite size that may be termed flocking black holes.

preprint2021arXiv

Uncertainty quantification in covid-19 spread: lockdown effects

We develop a Bayesian inference framework to quantify uncertainties in epidemiological models. We use SEIJR and SIJR models involving populations of susceptible, exposed, infective, diagnosed, dead and recovered individuals to infer from covid-19 data rate constants, as well as their variations in response to lockdown measures. To account for confinement, we distinguish two susceptible populations at different risk: confined and unconfined. We show that transmission and recovery rates within them vary in response to facts. A key unknown to predict the evolution of the epidemic is the fraction of the population affected by the virus, including asymptomatic subjects. Our study tracks its time evolution with quantified uncertainty from available official data from the onset of the epidemic, limited, however, by the data quality. We exemplify the technique with data from Spain, country in which late drastic lockdowns were enforced for months. In late actions and in the absence of other measures, spread is delayed but not stopped unless a large enough fraction of the population is confined until the asymptomatic population is depleted. To some extent, confinement could be replaced by strong distancing through masks in adequate circumstances.

preprint2016arXiv

Noninvasive imaging of three dimensional micro and nanostructures by topological methods

We present topological derivative and energy based procedures for the imaging of micro and nanostructures using one beam of visible light of a single wavelength. Objects with diameters as small as 10 nm can be located, and their position tracked with nanometer precision. Multiple objects distributed either on planes perpendicular to the incidence direction or along axial lines in the incidence direction are distinguishable. More precisely, the shape and size of plane sections perpendicular to the incidence direction can be clearly determined, even for asymmetric and non-convex scatterers. Axial resolution improves as the size of the objects decreases. Initial reconstructions may proceed by glueing together 2D horizontal slices between axial peaks or by locating objects at 3D peaks of topological energies, depending on the effective wavenumber. Below a threshold size, topological derivative based iterative schemes improve initial predictions of the location, size and shape of objects by postprocessing fixed measured data. For larger sizes, tracking the peaks of topological energy fields that average information from additional incident light beams seems to be more effective.

preprint2015arXiv

Differential growth of wrinkled biofilms

Biofilms are antibiotic-resistant bacterial aggregates that grow on moist surfaces and can trigger hospital-acquired infections. They provide a classical example in biology where the dynamics of cellular communities may be observed and studied. Gene expression regulates cell division and differentiation, which affect the biofilm architecture. Mechanical and chemical processes shape the resulting structure. We gain insight into the interplay between cellular and mechanical processes during biofilm development on air-agar interfaces by means of a hybrid model. Cellular behavior is governed by stochastic rules informed by a cascade of concentration fields for nutrients, waste and autoinducers. Cellular differentiation and death alter the structure and the mechanical properties of the biofilm, which is deformed according to Foppl-Von Karman equations informed by cellular processes and the interaction with the substratum. Stiffness gradients due to growth and swelling produce wrinkle branching. We are able to reproduce wrinkled structures often formed by biofilms on air-agar interfaces, as well as spatial distributions of differentiated cells commonly observed with B. subtilis.

preprint2015arXiv

Measuring strain and rotation fields at the dislocation core in graphene

Strain fields, dislocations and defects may be used to control electronic properties of graphene. By using advanced imaging techniques with high-resolution transmission electron microscopes, we have measured the strain and rotation fields about dislocations in monolayer graphene with single-atom sensitivity. These fields differ qualitatively from those given by conventional linear elasticity. However, atom positions calculated from two dimensional (2D) discrete elasticity and three dimensional discrete periodized Föppl-von Kármán equations (dpFvKEs) yield fields close to experiments when determined by geometric phase analysis. 2D theories produce symmetric fields whereas those from experiments exhibit asymmetries. Numerical solutions of dpFvKEs provide strain and rotation fields of dislocation dipoles and pairs that also exhibit asymmetries and, compared with experiments, may yield information on out-of-plane displacements of atoms. While discrete theories need to be solved numerically, analytical formulas for strains and rotation about dislocations can be obtained from 2D Mindlin's hyperstress theory. These formulas are very useful for fitting experimental data and provide a template to ascertain the importance of nonlinear and nonplanar effects. Measuring the parameters of this theory, we find two characteristic lengths between three and four times the lattice spacings that control dilatation and rotation about a dislocation. At larger distances from the dislocation core, the elastic fields decay to those of conventional elasticity. Our results may be relevant for strain engineering in graphene and other 2D materials of current interest.

preprint2013arXiv

Sawtooth patterns in biomolecules force-extension curves: an equilibrium-statistical-mechanics theory

We analyze the force-extension curve for a general class of systems, which are described at the mesoscopic level by a free energy depending on the extension of its components. Similarly to what is done in real experiments, the total length of the system is the controlled parameter. This imposes a global constraint in the minimization procedure leading to the equilibrium values of the extensions. As a consequence, the force-extension curve has multiple branches in a certain range of forces. The stability of these branches is governed by the free energy: there are a series of first-order phase transitions at certain values of the total length, in which the free energy itself is continuous but its first derivative, the force, has a finite jump. This behavior is completely similar to the one observed in real experiments with biomolecules like proteins, and other complex systems.

preprint2012arXiv

Model of ripples in graphene

We propose a model of ripples in suspended graphene sheets based on plate equations that are made discrete with the periodicity of the honeycomb lattice and then periodized. In addition, the equation for the displacements with respect to the planar configuration contains a double-well site potential, a nonlinear friction and a multiplicative white noise term satisfying the fluctuation-dissipation theorem. The nonlinear friction terms agree with those proposed by Eichler et al [Nature Nanotech. {\bf 6}, 339 (2011)] to explain their experiments with a graphene resonator. The site double-well potential indicates that the carbon atoms at each lattice point have equal probability to move upward or downward off-plane. For the considered parameter values, the relaxation time due to friction is much larger than the periods of membrane vibrations and the noise is quite small. Then ripples with no preferred orientation appear as long-lived metastable states for any temperature. Numerical solutions confirm this picture.

preprint2012arXiv

Ripples in a graphene membrane coupled to Glauber spins

We propose a theory of ripples in suspended graphene sheets based on two-dimensional elasticity equations that are made discrete on the honeycomb lattice and then periodized. At each point carbon atoms are coupled to Ising spins whose values indicate the atoms local trend to move vertically off-plane. The Ising spins are in contact with a thermal bath and evolve according to Glauber dynamics. In the limit of slow spin flip compared to membrane vibrations, ripples with no preferred orientation appear as long-lived metastable states for any temperature. Numerical solutions confirm this picture.

preprint2012arXiv

Spin-oscillator model for DNA/RNA unzipping by mechanical force

We model unzipping of DNA/RNA molecules subject to an external force by a spin-oscillator system. The system comprises a macroscopic degree of freedom, represented by a one-dimensional oscillator, and internal degrees of freedom, represented by Glauber spins with nearest-neighbor interaction and a coupling constant proportional to the oscillator position. At a critical value $F_c$ of an applied external force $F$, the oscillator rest position (order parameter) changes abruptly and the system undergoes a first-order phase transition. When the external force is cycled at different rates, the extension given by the oscillator position exhibits a hysteresis cycle at high loading rates whereas it moves reversibly over the equilibrium force-extension curve at very low loading rates. Under constant force, the logarithm of the residence time at the stable and metastable oscillator rest position is proportional to $(F-F_c)$ as in an Arrhenius law.

preprint2012arXiv

Strain and rotation fields of dislocations in graphene

Strain and rotation fields of dislocations in monolayer graphene have been mapped in a recent experiment. These fields are finite everywhere and differ from those given by linear elasticity which does not consider rotation explicitly and predicts infinite rotation and strains at the dislocation point. A derivative regularization of two-dimensional linear elasticity fixes these shortcomings. The theory adds rotation, dislocation and residual strain energies to the usual elastic energy. There are two extra material constants that determine the size of the dislocation core and are determined from experimental data. These findings are useful for studies of dislocations in graphene and for analyzing incipient plasticity in two dimensional crystals.

preprint2011arXiv

Statics and dynamics of a harmonic oscillator coupled to a one-dimensional Ising system

We investigate an oscillator linearly coupled with a one-dimensional Ising system. The coupling gives rise to drastic changes both in the oscillator statics and dynamics. Firstly, there appears a second order phase transition, with the oscillator stable rest position as its order parameter. Secondly, for fast spins, the oscillator dynamics is described by an effective equation with a nonlinear friction term that drives the oscillator towards the stable equilibrium state.

preprint2011arXiv

Theory of defect dynamics in graphene: defect groupings and their stability

We use our theory of periodized discrete elasticity to characterize defects in graphene as the cores of dislocations or groups of dislocations. Earlier numerical implementations of the theory predicted some of the simpler defect groupings observed in subsequent Transmission Electron Microscope experiments. Here we derive the more complicated defect groupings of three or four defect pairs from our theory, show that they correspond to the cores of two pairs of dislocation dipoles and ascertain their stability.

preprint2010arXiv

Nonequilibrium dynamics of a fast oscillator coupled to Glauber spins

A fast harmonic oscillator is linearly coupled with a system of Ising spins that are in contact with a thermal bath, and evolve under a slow Glauber dynamics at dimensionless temperature $θ$. The spins have a coupling constant proportional to the oscillator position. The oscillator-spin interaction produces a second order phase transition at $θ=1$ with the oscillator position as its order parameter: the equilibrium position is zero for $θ>1$ and non-zero for $θ< 1$. For $θ<1$, the dynamics of this system is quite different from relaxation to equilibrium. For most initial conditions, the oscillator position performs modulated oscillations about one of the stable equilibrium positions with a long relaxation time. For random initial conditions and a sufficiently large spin system, the unstable zero position of the oscillator is stabilized after a relaxation time proportional to $θ$. If the spin system is smaller, the situation is the same until the oscillator position is close to zero, then it crosses over to a neighborhood of a stable equilibrium position about which keeps oscillating for an exponentially long relaxation time. These results of stochastic simulations are predicted by modulation equations obtained from a multiple scale analysis of macroscopic equations.

preprint2010arXiv

Phase transitions in a mechanical system coupled to Glauber spins

A harmonic oscillator linearly coupled with a linear chain of Ising spins is investigated. The $N$ spins in the chain interact with their nearest neighbours with a coupling constant proportional to the oscillator position and to $N^{-1/2}$, are in contact with a thermal bath at temperature $T$, and evolve under Glauber dynamics. The oscillator position is a stochastic process due to the oscillator-spin interaction which produces drastic changes in the equilibrium behaviour and the dynamics of the oscillator. Firstly, there is a second order phase transition at a critical temperature $T_c$ whose order parameter is the oscillator stable rest position: this position is zero above $T_c$ and different from zero below $T_c$. This transition appears because the oscillator moves in an effective potential equal to the harmonic term plus the free energy of the spin system at fixed oscillator position. Secondly, assuming fast spin relaxation (compared to the oscillator natural period), the oscillator dynamical behaviour is described by an effective equation containing a nonlinear friction term that drives the oscillator towards the stable equilibrium state of the effective potential. The analytical results are compared with numerical simulation throughout the paper.

preprint2006arXiv

Kinetics of helium bubble formation in nuclear materials

The formation and growth of helium bubbles due to self-irradiation in plutonium has been modelled by a discrete kinetic equations for the number densities of bubbles having $k$ atoms. Analysis of these equations shows that the bubble size distribution function can be approximated by a composite of: (i) the solution of partial differential equations describing the continuum limit of the theory but corrected to take into account the effects of discreteness, and (ii) a local expansion about the advancing leading edge of the distribution function in size space. Both approximations contribute to the memory term in a close integrodifferential equation for the monomer concentration of single helium atoms. The present boundary layer theory for discrete equations is compared to the numerical solution of the full kinetic model and to previous approximation of Schaldach and Wolfer involving a truncated system of moment equations.