Source author record

Isaac Pérez Castillo

Isaac Pérez Castillo 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

27works
21topics
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

27 published item(s)

preprint2022arXiv

Roadmap for Optical Tweezers

Optical tweezers are tools made of light that enable contactless pushing, trapping, and manipulation of objects ranging from atoms to space light sails. Since the pioneering work by Arthur Ashkin in the 1970s, optical tweezers have evolved into sophisticated instruments and have been employed in a broad range of applications in life sciences, physics, and engineering. These include accurate force and torque measurement at the femtonewton level, microrheology of complex fluids, single micro- and nanoparticle spectroscopy, single-cell analysis, and statistical-physics experiments. This roadmap provides insights into current investigations involving optical forces and optical tweezers from their theoretical foundations to designs and setups. It also offers perspectives for applications to a wide range of research fields, from biophysics to space exploration.

preprint2022arXiv

Spectral properties of the generalized diluted Wishart ensemble

The celebrated Marčenko-Pastur law, that considers the asymptotic spectral density of random covariance matrices, has found a great number of applications in physics, biology, economics, engineering, among others. Here, using techniques from statistical mechanics of spin glasses, we derive simple formulas concerning the spectral density of generalized diluted Wishart matrices. These are defined as $\bm{F}\equiv \frac{1}{2d}\left( \bm{X}\bm{Y}^T+ \bm{Y}\bm{X}^T\right)$, where $\bm{X}$ and $\bm{Y}$ are diluted $N\times P$ rectangular matrices, whose entries correspond to the links of doubly-weighted random bipartite Poissonian graphs following the distribution $P(x_{i}^μ,y_{i}^μ)=\frac{d}{N}\varrho(x_{i}^μ,y_{i}^μ)+\left(1-\frac{d}{N}\right)δ_{x_{i}^μ,0}δ_{y_{i}^μ,0}$, with the probability density $\varrho(x,y)$ controlling the correlation between the matrices entries of $\bm{X}$ and $\bm{Y}$. Our results cover several interesting cases by varying the parameters of the matrix ensemble, namely, the dilution of the graph $d$, the rectangularity of the matrices $α=N/P$, and the degree of correlation of the matrix entries via the density $\varrho(x,y)$. Finally, we compare our findings to numerical diagonalisation showing excellent agreement.

preprint2020arXiv

Optical Tweezers: A Comprehensive Tutorial from Calibration to Applications

Since their invention in 1986 by Arthur Ashkin and colleagues, optical tweezers have become an essential tool in several fields of physics, spectroscopy, biology, nanotechnology, and thermodynamics. In this Tutorial, we provide a primer on how to calibrate optical tweezers and how to use them for advanced applications. After a brief general introduction on optical tweezers, we focus on describing and comparing the various available calibration techniques. Then, we discuss some cutting-edge applications of optical tweezers in a liquid medium, namely to study single-molecule and single-cell mechanics, microrheology, colloidal interactions, statistical physics, and transport phenomena. Finally, we consider optical tweezers in vacuum, where the absence of a viscous medium offers vastly different dynamics and presents new challenges. We conclude with some perspectives for the field and the future application of optical tweezers. This Tutorial provides both a step-by-step guide ideal for non-specialists entering the field and a comprehensive manual of advanced techniques useful for expert practitioners. All the examples are complemented by the sample data and software necessary to reproduce them.

preprint2020arXiv

The exact two-spinon longitudinal dynamical structure factor of the anisotropic XXZ model

Inelastic neutron scattering experiments are commonly used to unveil how excitations on Heisenberg spin models play a role in dynamical correlations functions. For a certain class of materials, like CsCoCl$_3$ or CsCoBr$_3$ salts, it turns out that their magnetic properties are fairly well approximated by quasi-one dimensional XXZ models, which enjoy the property of quantum integrability. In these instances, one can in principle use their underlying algebraic structure to describe very precisely how excitations, the so-called spinons, participate in dynamical correlations functions. Even though the available theories (either algebraic Bethe ansatz or quantum group approach) provide all the needed physical quantities such as form factors, complete set of eigenstates and spectrum, it is typically a rather daunting task, however, to obtain sufficiently simple analytical expressions for computing Dynamical Structure Factors (DSFs), valuable, e.g., for parameters estimation based on experimental data. This is particularly the case for the longitudinal DSF of the XXZ model, which has eluded a formal mathematical treatment thus far. Using the quantum group approach, we present here an exact and simple expression of the 2-spinon longitudinal DSF and show our results to be consistent with the expected sum rules and the isotropic and Ising antiferromagnet limiting cases.

preprint2020arXiv

Using posterior predictive distributions to analyse epidemic models: COVID-19 in Mexico City

Epidemiological models contain a set of parameters that must be adjusted based on available observations. Once a model has been calibrated, it can be used as a forecasting tool to make predictions and to evaluate contingency plans. It is customary to employ only point estimators for such predictions. However, some models may fit the same data reasonably well for a broad range of parameter values, and this flexibility means that predictions stemming from such models will vary widely, depending on the particular parameter values employed within the range that give a good fit. When data are poor or incomplete, model uncertainty widens further. A way to circumvent this problem is to use Bayesian statistics to incorporate observations and use the full range of parameter estimates contained in the posterior distribution to adjust for uncertainties in model predictions. Specifically, given the epidemiological model and a probability distribution for observations, we use the posterior distribution of model parameters to generate all possible epidemiological curves via the posterior predictive distribution. From the envelope of all curves one can extract the worst-case scenario and study the impact of implementing contingency plans according to this assessment. We apply this approach to the potential evolution of COVID-19 in Mexico City and assess whether contingency plans are being successful and whether the epidemiological curve has flattened.

preprint2019arXiv

Phase transitions in atypical systems induced by a condensation transition on graphs

Random graphs undergo structural phase transitions that are crucial for dynamical processes and cooperative behavior of models defined on graphs. In this work we investigate the impact of a first-order structural transition on the thermodynamics of the Ising model defined on Erdös-Rényi random graphs, as well as on the eigenvalue distribution of the adjacency matrix of the same graphical model. The structural transition in question yields graph samples exhibiting condensation, characterized by a large number of nodes having degrees in a narrow interval. We show that this condensation transition induces distinct thermodynamic first-order transitions between the paramagnetic and the ferromagnetic phases of the Ising model. The condensation transition also leads to an abrupt change in the global eigenvalue statistics of the adjacency matrix, which renders the second moment of the eigenvalue distribution discontinuous. As a side result, we derive the critical line determining the percolation transition in Erdös-Rényi graph samples that feature condensation of degrees.

preprint2016arXiv

Large deviation function for the number of eigenvalues of sparse random graphs inside an interval

We present a general method to obtain the exact rate function $Ψ_{[a,b]}(k)$ controlling the large deviation probability $\text{Prob}[\mathcal{I}_N[a,b]=kN] \asymp e^{-NΨ_{[a,b]}(k)}$ that a $N \times N$ sparse random matrix has $\mathcal{I}_N[a,b]=kN$ eigenvalues inside the interval $[a,b]$. The method is applied to study the eigenvalue statistics in two distinct examples: (i) the shifted index number of eigenvalues for an ensemble of Erdös-Rényi graphs and (ii) the number of eigenvalues within a bounded region of the spectrum for the Anderson model on regular random graphs. A salient feature of the rate function in both cases is that, unlike rotationally invariant random matrices, it is asymmetric with respect to its minimum. The asymmetric character depends on the disorder in a way that is compatible with the distinct eigenvalue statistics corresponding to localized and delocalized eigenstates. The results also show that the number variance $σ_{N}^{2}$ for the Anderson model on a regular graph scales as $σ_{N}^{2} \propto N$ ($N \gg 1$) for any nonzero disorder, which is consistent with the absence of level-repulsion in the extended phase. Our theoretical findings are thoroughly compared to numerical diagonalization in both cases, showing a reasonable good agreement.

preprint2016arXiv

When does inequality freeze an economy?

Inequality and its consequences are the subject of intense recent debate. Using a simplified model of the economy, we address the relation between inequality and liquidity, the latter understood as the frequency of economic exchanges. Assuming a Pareto distribution of wealth for the agents, that is consistent with empirical findings, we find an inverse relation between wealth inequality and overall liquidity. We show that an increase in the inequality of wealth results in an even sharper concentration of the liquid financial resources. This leads to a congestion of the flow of goods and the arrest of the economy when the Pareto exponent reaches one.

preprint2015arXiv

A generalised Airy distribution function for the accumulated area swept by $N$ vicious Brownian paths

In this work exact expressions for the distribution function of the accumulated area swept by reunions and meanders of $N$ vicious Brownian particles up to time $T$ are derived. The results are expressed in terms of a generalised Airy distribution function, containing the Vandermonde determinant of the Airy roots. By mapping the problem to an Random Matrix Theory ensemble we are able to perform Monte Carlo simulations finding perfect agreement with the theoretical results.

preprint2015arXiv

A unified approach for large deviations of bulk and extreme eigenvalues of the Wishart ensemble

Within the framework of the Coulomb fluid picture, we present a unified approach to derive the large deviations of bulk and extreme eigenvalues of large Wishart matrices. By analysing the statistics of the shifted index number we are able to derive a rate function $Ψ(c, x)$ depending on two variables: the fraction $c$ of eigenvalues to the left of an infinite energetic barrier at position $x$. For a fixed value of $c$, the rate function gives the large deviations of the bulk eigenvalues. In particular, in the limits $c\to 0$ or $c\to 1$ it is possible to extract the left and right deviations of the smallest and largest eigenvalues, respectively. Alternatively, for a fixed value $x$ of the barrier, the rate function provides the large deviations of the shifted index number. All our analytical findings are compared with Metropolis Monte Carlo simulations, obtaining excellent agreement.

preprint2015arXiv

Large deviations of the shifted index number in the Gaussian ensemble

We show that, using the Coulomb fluid approach, we are able to derive a rate function $Ψ(c,x)$ of two variables that captures: (i) the large deviations of bulk eigenvalues; (ii) the large deviations of extreme eigenvalues (both left and right large deviations); (iii) the statistics of the fraction $c$ of eigenvalues to the left of a position $x$. Thus, $Ψ(c,x)$ explains the full order statistics of the eigenvalues of large random Gaussian matrices as well as the statistics of the shifted index number. All our analytical findings are thoroughly compared with Monte Carlo simulations, obtaining excellent agreement. A summary of preliminary results was already presented in [22] in the context of one-dimensional trapped spinless fermions in a harmonic potential.

preprint2015arXiv

Minimal cooling speed for glass transition in a simple solvable energy landscape model

The minimal cooling speed required to form a glass is obtained for a simple solvable energy landscape model. The model, made from a two-level system modified to include the topology of the energy landscape, is able to capture either a glass transition or a crystallization depending on cooling rate. In this setup, the minimal cooling speed to achieve glass formation is then found to be related with the relaxation time and with the thermal history. In particular, we obtain that the thermal history encodes small fluctuations around the equilibrium population which are exponentially amplified near the glass transition, which mathematically corresponds to the boundary layer of the master equation. Finally, to verify our analytical results, a kinetic Monte-Carlo simulation was implemented.

preprint2015arXiv

On the passage properties of the gradual capture of a diffusive particle in the presence of drift

We investigate a stochastic process consisting of a two-dimensional particle with anisotropic diffusion, mortality rate and a drift velocity, in the presence of an absorbing boundary. After the particle has encountered the boundary, the process is restarted with updated values of its diffusion constants and drift velocity. We then derive the joint probability that, after $M$ encounters, the particle is absorbed at a point of the boundary at a given time and show that, under certain conditions, the eventual hitting probability corresponds to a Bessel distribution. In the context of the man-mosquitoes problem, the mosquito is able to gradually capture the man, after which the mosquito follows a diffusion process with no drift. Our results are compared thoroughly with Monte Carlo simulations showing excellent agreement.

preprint2014arXiv

Cavity approach for modeling and fitting polymer stretching

The mechanical properties of molecules are today captured by single molecule manipulation experiments, so that polymer features are tested at a nanometric scale. Yet devising mathematical models to get further insight beyond the commonly studied force--elongation relation is typically hard. Here we draw from techniques developed in the context of disordered systems to solve models for single and double--stranded DNA stretching in the limit of a long polymeric chain. Since we directly derive the marginals for the molecule local orientation, our approach allows us to readily calculate the experimental elongation as well as other observables at wish. As an example, we evaluate the correlation length as a function of the stretching force. Furthermore, we are able to fit successfully our solution to real experimental data. Although the model is admittedly phenomenological, our findings are very sound. For single--stranded DNA our solution yields the correct (monomer) scale and, yet more importantly, the right persistence length of the molecule. In the double--stranded case, our model reproduces the well-known overstretching transition and correctly captures the ratio between native DNA and overstretched DNA. Also in this case the model yields a persistence length in good agreement with consensus, and it gives interesting insights into the bending stiffness of the native and overstretched molecule, respectively.

preprint2014arXiv

Phase transitions in the condition number distribution of Gaussian random matrices

We study the statistics of the condition number $κ=λ_{\mathrm{max}}/λ_{\mathrm{min}}$ (the ratio between largest and smallest squared singular values) of $N\times M$ Gaussian random matrices. Using a Coulomb fluid technique, we derive analytically and for large $N$ the cumulative $\mathcal{P}[κ<x]$ and tail-cumulative $\mathcal{P}[κ>x]$ distributions of $κ$. We find that these distributions decay as $\mathcal{P}[κ<x]\approx\exp\left(-βN^2 Φ_{-}(x)\right)$ and $\mathcal{P}[κ>x]\approx\exp\left(-βN Φ_{+}(x)\right)$, where $β$ is the Dyson index of the ensemble. The left and right rate functions $Φ_{\pm}(x)$ are independent of $β$ and calculated exactly for any choice of the rectangularity parameter $α=M/N-1>0$. Interestingly, they show a weak non-analytic behavior at their minimum $\langleκ\rangle$ (corresponding to the average condition number), a direct consequence of a phase transition in the associated Coulomb fluid problem. Matching the behavior of the rate functions around $\langleκ\rangle$, we determine exactly the scale of typical fluctuations $\sim\mathcal{O}(N^{-2/3})$ and the tails of the limiting distribution of $κ$. The analytical results are in excellent agreement with numerical simulations.

preprint2014arXiv

Spectral density of products of Wishart dilute random matrices. Part I: the dense case

In this work we study the spectral density of products of Wishart diluted random matrices of the form $X(1)\cdots X(M)(X(1)\cdots X(M))^T$ using the Edwards-Jones trick to map this problem into a system of interacting particles with random couplings on a multipartite graph. We apply the cavity method to obtain recursive relations in typical instances from which to obtain the spectral density. As this problem is fairly rich, we start by reporting in part I a lengthy analysis for the case of dense matrices. Here we derive that the spectral density is a solution of a polynomial equation of degree $M+1$ and obtain exact expressions of it for $M=1$, $2$ and $3$. For general $M$, we are able to find the exact expression of the spectral density only when all the matrices $X(t)$ for $t=1,\ldots, M$ are square. We also make some observations for general $M$, based admittedly on some weak numerical evidence, which we expect to be correct.

preprint2014arXiv

Spectral order statistics of Gaussian random matrices: large deviations for trapped fermions and associated phase transitions

We compute the full order statistics of a one-dimensional gas of fermions in a harmonic trap at zero temperature, including its large deviation tails. The problem amounts to computing the probability distribution of the $k$th smallest eigenvalue $λ_{(k)}$ of a large dimensional Gaussian random matrix. We find that this probability behaves for large $N$ as $\mathcal{P}[λ_{(k)}=x]\approx \exp\left(-βN^2 ψ(k/N,x)\right)$, where $β$ is the Dyson index of the ensemble. The rate function $ψ(c,x)$, computed explicitly as a function of $x$ in terms of the intensive label $c=k/N$, has a quadratic behavior modulated by a weak logarithmic singularity at its minimum. This is shown to be related to phase transitions in the associated Coulomb gas problem. The connection with statistics of extreme eigenvalues of random matrices is also elucidated.

preprint2013arXiv

A Novel Methodology to Estimate Metabolic Flux Distributions in Constraint-Based Models

Quite generally, constraint-based metabolic flux analysis describes the space of viable flux configurations for a metabolic network as a high-dimensional polytope defined by the linear constraints that enforce the balancing of production and consumption fluxes for each chemical species in the system. In some cases, the complexity of the solution space can be reduced by performing an additional optimization, while in other cases, knowing the range of variability of fluxes over the polytope provides a sufficient characterization of the allowed configurations. There are cases, however, in which the thorough information encoded in the individual distributions of viable fluxes over the polytope is required. Obtaining such distributions is known to be a highly challenging computational task when the dimensionality of the polytope is sufficiently large, and the problem of developing cost-effective {\it ad hoc} algorithms has recently seen a major surge of interest. Here, we propose a method that allows us to perform the required computation heuristically in a time scaling {\it linearly} with the number of reactions in the network, overcoming some limitations of similar techniques employed in recent years. As a case study, we apply it to the analysis of the human red blood cell metabolic network, whose solution space can be sampled by different exact techniques, like Hit-and-Run Monte Carlo (scaling roughly like the third power of the system size). Remarkably accurate estimates for the true distributions of viable reaction fluxes are obtained, suggesting that, although further improvements are desirable, our method enhances our ability to analyze the space of allowed configurations for large biochemical reaction networks.

preprint2013arXiv

Reunion probabilities of $N$ one-dimensional random walkers with mixed boundary conditions

In this work we extend the results of the reunion probability of $N$ one-dimensional random walkers to include mixed boundary conditions between their trajectories. The level of the mixture is controlled by a parameter $c$, which can be varied from $c=0$ (independent walkers) to $c\to\infty$ (vicious walkers). The expressions are derived by using Quantum Mechanics formalism (QMf) which allows us to map this problem into a Lieb-Liniger gas (LLg) of $N$ one-dimensional particles. We use Bethe ansatz and Gaudin's conjecture to obtain the normalized wave-functions and use this information to construct the propagator. As it is well-known, depending on the boundary conditions imposed at the endpoints of a line segment, the statistics of the maximum heights of the reunited trajectories have some connections with different ensembles in Random Matrix Theory (RMT). Here we seek to extend those results and consider four models: absorbing, periodic, reflecting, and mixed. In all four cases, the probability that the maximum height is less or equal than $L$ takes the form $F_N(L)=A_N\sum_{k\inΩ_{B}}\int Dz e^{-\sum_{j=1}^Nk_j^2+G_N(k)-\sum_{j,\ell=1}^N z_jV_{j\ell}(k)\overline{z}_\ell}$, where $A_N$ is a normalization constant, $G_N(k)$ and $V_{j\ell}(k)$ depend on the type of boundary condition, and $Ω_{B}$ is the solution set of quasi-momenta $k$ obeying the Bethe equations for that particular boundary condition.

preprint2012arXiv

A weighted belief-propagation algorithm to estimate volume-related properties of random polytopes

In this work we introduce a novel weighted message-passing algorithm based on the cavity method to estimate volume-related properties of random polytopes, properties which are relevant in various research fields ranging from metabolic networks, to neural networks, to compressed sensing. Unlike the usual approach consisting in approximating the real-valued cavity marginal distributions by a few parameters, we propose an algorithm to faithfully represent the entire marginal distribution. We explain various alternatives to implement the algorithm and benchmark the theoretical findings by showing concrete applications to random polytopes. The results obtained with our approach are found to be in very good agreement with the estimates produced by the Hit-and-Run algorithm, known to produce uniform sampling.

preprint2012arXiv

Spectral Properties of the Jacobi Ensembles via the Coulomb Gas approach

Using the Coulomb gas method and standard methods of statistical physics, we compute analytically the joint cumulative probability distribution of the extreme eigenvalues of the Jacobi-MANOVA ensemble of random matrices, in the limit of large matrices. This allows us to derive the rate functions for the large fluctuations to the left and the right of the expected values of the smallest and largest eigenvalues analytically. Our findings are compared with some available known exact results as well as with numerical simulations finding good agreement.

preprint2011arXiv

A weighted message-passing algorithm to estimate volume-related properties of random polytopes

In this letter, we introduce a novel message-passing algorithm for a class of problems which can be mathematically understood as estimating volume-related properties of random polytopes. Unlike the usual approach consisting in approximating the real-valued cavity marginal distributions by a few parameters, we propose a weighted message-passing algorithm to deal with the entire function. Various alternatives of how to implement our approach are discussed and numerical results for random polytopes are compared with results using the Hit-and-Run algorithm.

preprint2011arXiv

The Storm and Nelson's model for polymer stretching revisited

The quite recent technological rise in molecular biology allowed single molecule manipulation experiments, where molecule stretching plays a primary role. In order to understand the experimental data, it is felt the urge of some physical and mathematical models to quantitatively express the mechanical properties of the observed molecules. In this paper we reconsider a simple phenomenological model which reproduces the behaviour of single and double stranded DNA under tension. The problem is easily solved via the cavity method both in the small forces range and in presence of overstretching transition, so that some properties such as bending stiffness and elasticity of DNA emerge in a very clear manner. Our theoretical findings are successfully fitted to real measurements and compared to Monte Carlo simulations, confirming the quality of the approach.

preprint2010arXiv

Large Deviations of the Smallest Eigenvalue of the Wishart-Laguerre Ensemble

We consider the large deviations of the smallest eigenvalue of the Wishart-Laguerre Ensemble. Using the Coulomb gas picture we obtain rate functions for the large fluctuations to the left and the right of the hard edge. Our findings are compared with known exact results for $β=1$ finding good agreement. We also consider the case of almost square matrices finding new universal rate functions describing large fluctuations.

preprint2009arXiv

Spectral density of random graphs with topological constraints

The spectral density of random graphs with topological constraints is analysed using the replica method. We consider graph ensembles featuring generalised degree-degree correlations, as well as those with a community structure. In each case an exact solution is found for the spectral density in the form of consistency equations depending on the statistical properties of the graph ensemble in question. We highlight the effect of these topological constraints on the resulting spectral density.

preprint2005arXiv

Cavity approach for real variables on diluted graphs and application to synchronization in small-world lattices

We study XY spin systems on small world lattices for a variety of graph structures, e.g. Poisson and scale-free, superimposed upon a one dimensional chain. In order to solve this model we extend the cavity method in the one pure-state approximation to deal with real-valued dynamical variables. We find that small-world architectures significantly enlarge the region in parameter space where synchronization occurs. We contrast the results of population dynamics performed on a truncated set of cavity fields with Monte Carlo simulations and find excellent agreement. Further, we investigate the appearance of replica symmetry breaking in the spin-glass phase by numerically analyzing the proliferation of pure states in the message passing equations.

preprint2005arXiv

On exact mappings between fermionic Ising spin glass and classical spin glass models

We present in this paper exact analytical expressions for the thermodynamical properties and Green functions of a certain family of fermionic Ising spin-glass models with Hubbard interaction, by noticing that their Hamiltonian is a function of the number operator only. The thermodynamical properties are mapped to the classical Ghatak-Sherrington spin-glass model while the the Density of States (DoS) is related to its joint spin-field distribution. We discuss the presence of the pseudogap in the DoS with the help of this mapping.