BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Heiko Rieger

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softquant-phBiological Physicscond-mat.dis-nncond-mat.quant-gasSubcellular Processescond-mat.str-elphysics.comp-phphysics.flu-dyncond-matcond-mat.otherhep-latphysics.chem-ph

Catalog footprint

What is connected

34works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

34 published item(s)

preprint2026arXiv

Generalized Finite-time Optimal Control Framework in Stochastic Thermodynamics

Optimal processes in stochastic thermodynamics are a frontier for understanding the control and design of non-equilibrium systems, with broad practical applications in biology, chemistry, and nanoscale/mesoscale systems. Optimal mass transport theory and thermodynamic geometry have emerged as optimal control methodology, but they are based on slow-driving and close to equilibrium assumptions. An optimal control framework in stochastic thermodynamics for finite time driving is still elusive. Therefore, we solve in this paper an optimal control problem for changing the control parameters of a discrete-state far-from-equilibrium process from an initial to a final value in finite-time. Optimal driving protocols are derived that minimize the total finite-time dissipation cost for the driving process. Our framework reveals that discontinuous endpoint jumps are a generic, model-independent physical mechanism that minimizes the optimal driving entropy production, whose importance is further amplified for far-from-equilibrium systems. The thermodynamic and dynamic physical interpretation and understanding of discontinuous endpoint jumps is formulated. An exact mapping between the finite-time to slow driving optimal control formulation is elucidated, developing the state-of-the-art of optimal mass transport theory and thermodynamic geometry, which has been the current paradigm for studying optimal processes in stochastic thermodynamics that relies on slow driving assumptions. Our framework opens up a plethora of applications to the thermodynamically efficient control of a far-from-equilibrium system in finite-time, which opens up a way to their efficient design principles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Polar flocks with discretized directions: the active clock model approaching the Vicsek model

We consider the off-lattice two-dimensional $q$-state active clock model (ACM) as a natural discretization of the Vicsek model (VM) describing flocking. The ACM consists of particles able to move in the plane in a discrete set of $q$ equidistant angular directions, as in the active Potts model (APM), with an alignment interaction inspired by the ferromagnetic equilibrium clock model. We find that for a small number of directions, the flocking transition of the ACM has the same phenomenology as the APM, including macrophase separation and reorientation transition. For a larger number of directions, the flocking transition in the ACM becomes equivalent to the one of the VM and displays microphase separation and only transverse bands, i.e. no re-orientation transition. Concomitantly also the transition of the $q\to\infty$ limit of the ACM, the active XY model (AXYM), is in the same universality class as the VM. We also construct a coarse-grained hydrodynamic description for the ACM and AXYM akin to the VM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Suppression of discontinuous phase transitions by particle diffusion

We investigate the phase transitions of the $q$-state Brownian Potts model in two dimensions (2d) comprising Potts spins that diffuse like Brownian particles and interact ferromagnetically with other spins within a fixed distance. With extensive Monte Carlo simulations we find a continuous phase transition from a paramagnetic to a ferromagnetic phase even for $q>4$. This is in sharp contrast to the existence of a discontinuous phase transition in the equilibrium $q$-state Potts model in 2d with $q>4$. We present detailed numerical evidence for a continuous phase transition and argue that diffusion generated dynamical positional disorder suppresses phase coexistence leading to a continuous transition.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stochastic model of T Cell repolarization during target elimination (II)

Cytotoxic T lymphocytes (T cells) and natural killer cells form a tight contact, the immunological synapse (IS), with target cells, where they release their lytic granules containing perforin/granzyme and cytokine containing vesicles. During this process the cell repolarizes and moves the microtubule organizing center (MTOC) towards the IS. In the first part of our work we developed a a computational model for the molecular-motor-driven motion of the MT cytoskeleton during T cell polarization and analyzed effects of cortical sliding and capture-shrinkage mechanisms. Here we use this model to analyze the dynamics of the MTOC repositioning in situations in which a) the IS is in an arbitrary position with respect to the initial position of the MTOC and b) the T cell has two IS at two arbitrary positions. In the case of one IS, we found that the initial position determines which mechanism is dominant and that the time of repositioning does not rise monotonously with the MTOC-IS distance. In the case of two IS, we observe several scenarios that have also been reported experimentally: the MTOC alternates stochastically (but with a well defined average transition time) between the two IS; it wiggles in between the two IS without transiting to one of the two; or it is at some point pulled to one of the two IS and stays there. Our model allows to predict which scenario emerges in dependency of the mechanisms in action and the number of dyneins present. We report that the presence of capture-shrinkage mechanism at least in one IS is necessary to assure the transitions in every cell configuration. Moreover, the frequency of transitions does not decrease with the distance between the two IS and is the highest when both mechanisms are present in both IS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Flocking and reorientation transition in the 4-state active Potts model

We study the active 4-state Potts model (APM) on the square lattice in which active particles have four internal states corresponding to the four directions of motion. A local alignment rule inspired by the ferromagnetic 4-state Potts model and self-propulsion via biased diffusion according to the internal particle states leads to flocking at high densities and low noise. We compute the phase diagram of the APM and explore the flocking dynamics in the region, in which the high-density (liquid) phase coexists with the low-density (gas) phase and forms a fluctuating band of coherently moving particles. As a function of the particle self-propulsion velocity, a novel reorientation transition of the phase-separated profiles from transversal to longitudinal band motion is revealed, which is absent in the Vicsek model and the active Ising model. We further construct a coarse-grained hydrodynamic description of the model which validates the results for the microscopic model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Haldane Insulator in the 1D Nearest-Neighbor Extended Bose-Hubbard Model with Cavity-Mediated Long-Range Interactions

In the one-dimensional Bose-Hubbard model with on-site and nearest neighbor interactions, a gapped phase characterized by an exotic non-local order parameter emerges, the Haldane insulator. Bose-Hubbard models with cavity-mediated global range interactions display phase diagrams, which are very similar to those with nearest neighbor repulsive interactions, but the Haldane phase remains elusive there. Here we study the one-dimensional Bose-Hubbard model with nearest-neighbor and cavity-mediated global-range interactions and scrutinize the existence of a Haldane Insulator phase. With the help of extensive quantum Monte-Carlo simulations we find that in the Bose-Hubbard model with only cavity-mediated global-range interactions no Haldane phase exists. For a combination of both interactions, the Haldane Insulator phase shrinks rapidly with increasing strength of the cavity-mediated global-range interactions. Thus, in spite of the otherwise very similar behavior the mean-field like cavity-mediated interactions strongly suppress the non-local order favored by nearest neighbor repulsion in some regions of the phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic model of T Cell repolarization during target elimination (I)

Cytotoxic T lymphocytes (T) and natural killer (NK) cells are the main cytotoxic killer cells of the human body to eliminate pathogen-infected or tumorigenic cells (i.e. target cells). Once a NK or T cell has identified a target cell, they form a tight contact zone, the immunological synapse (IS). One then observes repolarization of the cell involving the rotation of the microtubule (MT) cytoskeleton and a movement of the microtubule organizing center (MTOC) to a position that is just underneath the plasma membrane at the IS. Concomitantly a massive relocation of organelles attached to MTs is observed, including the Golgi apparatus, lytic granules and mitochondria. Since the mechanism of this relocation is still elusive we devise a theoretical model for the molecular motor driven motion of the MT cytoskeleton confined between membrane and nucleus. We analyze scenarios currently discussed in the literature, the cortical sliding and the capture-shrinkage mechanisms, and compare quantitative predictions about the spatio-temporal evolution of MTOC position and MT cytoskeleton morphology with experiments. The model predicts the experimentally observed biphasic nature of the process due to an interplay between MT cytoskeleton geometry and motor forces and confirms the dominance of the capture-shrinkage over the cortical sliding mechanism when MTOC and IS are initially diametrically opposed. We also find that the two mechanisms act synergistically, reducing the resources necessary for repositioning. The localization of dyneins in the pSMAC facilitates their interaction with the MTs. Our model also opens a way to infer details of the dynein distribution from the experimentally observed features of the MT cytoskeleton dynamics. In a subsequent publication, we will address the issue of general initial configurations and situations in which the T cell established two immunological synapses.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Capillary action in scalar active matter

We study the capacity of active matter to rise in thin tubes against gravity and other related phenomena, like, wetting of vertical plates and spontaneous imbibition, where a wetting liquid is drawn into a porous medium. This capillary action or capillarity is well known in classical fluids and originates from attractive interactions between the liquid molecules and the container walls, and from the attraction of the liquid molecules among each other. We observe capillarity in a minimal model for scalar active matter with purely repulsive interactions, where an effective attraction emerges due to slowdown during collisions between active particles and between active particles and walls. Simulations indicate that the capillary rise in thin tubes is approximately proportional to the active sedimentation length $λ$ and that the wetting height of a vertical plate grows superlinear with $λ$. In a disordered porous medium the imbibition height scales as $\langle h\rangle\proptoλϕ_m$, where $ϕ_m$ is its packing fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Computational models for active matter

A variety of computational models have been developed to describe active matter at different length and time scales. The diversity of the methods and the challenges in modeling active matter---ranging from molecular motors and cytoskeletal filaments over artificial and biological swimmers on microscopic to groups of animals on macroscopic scales---mainly originate from their out-of-equilibrium character, multiscale nature, nonlinearity, and multibody interactions. In the present review, various modeling approaches and numerical techniques are addressed, compared, and differentiated to illuminate the innovations and current challenges in understanding active matter. The complexity increases from minimal microscopic models of dry active matter toward microscopic models of active matter in fluids. Complementary, coarse-grained descriptions and continuum models are elucidated. Microscopic details are often relevant and strongly affect collective behaviors, which implies that the selection of a proper level of modeling is a delicate choice, with simple models emphasizing universal properties and detailed models capturing specific features. Finally, current approaches to further advance the existing models and techniques to cope with real-world applications, such as complex media and biological environments, are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Reentrant Random Quantum Ising Antiferromagnet

We consider the quantum Ising chain with uniformly distributed random antiferromagnetic couplings $(1 \le J_i \le 2)$ and uniformly distributed random transverse fields ($Γ_0 \le Γ_i \le 2Γ_0$) in the presence of a homogeneous longitudinal field, $h$. Using different numerical techniques (DMRG, combinatorial optimisation and strong disorder RG methods) we explore the phase diagram, which consists of an ordered and a disordered phase. At one end of the transition line ($h=0,Γ_0=1$) there is an infinite disorder quantum fixed point, while at the other end ($h=2,Γ_0=0$) there is a classical random first-order transition point. Close to this fixed point, for $h>2$ and $Γ_0>0$ there is a reentrant ordered phase, which is the result of quantum fluctuations by means of an order through disorder phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical analysis of homogeneous and inhomogeneous intermittent search strategies

A random search is a stochastic process representing the random motion of a particle (denoted as the searcher) that is terminated when it reaches (detects) a target particle or area the first time. In intermittent search the random motion alternates between two or more motility modes, one of which is non-detecting. An example is the slow diffusive motion as the detecting mode and fast, directed ballistic motion as the non-detecting mode, which can lead to much faster detection than a purely diffusive search. The transition rate between the diffusive and the ballistic mode (and back) together with the probability distribution of directions for the ballistic motion defines a search strategy. If these transition rates and/or probability distributions depend on the spatial coordinates within the search domain it is a spatially inhomogeneous search strategy, if both are constant, it is a homogeneous one. Here we study the efficiency, measured in terms of the mean first-passage time, of spatially homogeneous and inhomogeneous search strategies for three paradigmatic search problems: 1) the narrow escape problem, where the searcher has to find a small area on the boundary of the search domain, 2) reaction kinetics, which involves the detection of an immobile target in the interior of a search domain, and 3) the reaction-escape problem, where the searcher first needs to find a diffusive target before it can escape through a narrow region on the boundary. Using families of spatially inhomogeneous search strategies, partially motivated by the spatial organization of the cytoskeleton in living cells with a centrosome, we show that they can be made almost always more efficient than homogeneous strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimality of spatially inhomogeneous search strategies

We consider random search processes alternating stochastically between diffusion and ballistic motion, in which the distribution function of ballistic motion directions varies from point to point in space. The specific space dependence of the directional distribution together with the switching rates between the two modes of motion establishes a spatially inhomogeneous search strategy. We show that the mean first passage times for several standard search problems - narrow escape, reaction partner finding, reaction-escape - can be minimized with a directional distribution that is reminiscent of the spatial organization of the cytoskeleton filaments of cells with a centrosome: radial ballistic transport from center to periphery and back, and ballistic transport in random directions within a concentric shell of thickness $Δ_{\rm opt}$ along the domain boundary. The results suggest that living cells realize efficient search strategies for various intracellular transport problems economically through a spatial cytoskeleton organization that involves radial microtubules in the central region and only a narrow actin cortex rather than a cell body filled with randomly oriented actin filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum relaxation and finite size effects in the XY chain in a transverse field after global quenches

We consider global quenches in the quantum XY chain in a transverse field and study the nonequilibrium relaxation of the magnetization and the correlation function as well as the entanglement entropy in finite systems. For quenches in the ordered phase, the exact results are well described by a semiclassical theory (SCT) in terms of ballistically moving quasiparticles. In the thermodynamic limit the SCT is exact for the entanglement entropy and its modified version following the method of Calabrese, Essler and Fagotti: arXiv:1204.3911 is exact for the magnetization and the correlation function, too. The stationary correlation function is shown to be described by a generalized Gibbs ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Run-and-pause dynamics of cytoskeletal motor proteins

Cytoskeletal motor proteins are involved in major intracellular transport processes which are vital for maintaining appropriate cellular function. The motor exhibits distinct states of motility: active motion along filaments, and effectively stationary phase in which it detaches from the filaments and performs passive diffusion in the vicinity of the detachment point due to cytoplasmic crowding. The transition rates between motion and pause phases are asymmetric in general, and considerably affected by changes in environmental conditions which influences the efficiency of cargo delivery to specific targets. By considering the motion of molecular motor on a single filament as well as a dynamic filamentous network, we present an analytical model for the dynamics of self-propelled particles which undergo frequent pause phases. The interplay between motor processivity, structural properties of filamentous network, and transition rates between the two states of motility drastically changes the dynamics: multiple transitions between different types of anomalous diffusive dynamics occur and the crossover time to the asymptotic diffusive or ballistic motion varies by several orders of magnitude. We map out the phase diagrams in the space of transition rates, and address the role of initial conditions of motion on the resulting dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spatial Organization of the Cytoskeleton enhances Cargo Delivery to Specific Target Areas on the Plasma Membrane of Spherical Cells

Intracellular transport is vital for the proper functioning and survival of a cell. Cargo (proteins, vesicles, organelles, etc.) is transferred from its place of creation to its target locations via molecular motor assisted transport along cytoskeletal filaments. The transport efficiency is strongly affected by the spatial organization of the cytoskeleton, which constitutes an inhomogeneous, complex network. In cells with a centrosome microtubules grow radially from the central microtubule organizing center towards the cell periphery whereas actin filaments form a dense meshwork, the actin cortex, underneath the cell membrane with a broad range of orientations. The emerging ballistic motion along filaments is frequently interrupted due to constricting intersection nodes or cycles of detachment and reattachment processes in the crowded cytoplasm. In order to investigate the efficiency of search strategies established by the cell's specific spatial organization of the cytoskeleton we formulate a random velocity model with intermittent arrest states. With extensive computer simulations we analyze the dependence of the mean first passage times for narrow escape problems on the structural characteristics of the cytoskeleton, the motor properties and the fraction of time spent in each state. We find that an inhomogeneous architecture with a small width of the actin cortex constitutes an efficient intracellular search strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ultracold Bosons with cavity-mediated long-range interactions: A local mean field analysis of the phase diagram

Ultracold bosonic atoms in optical lattices self-organize into a variety of structural and quantum phases when placed into a single-mode cavity and pumped by a laser. Cavity optomechanical effects induce an atom density modulation at the cavity-mode wave length that competes with the optical lattice arrangement. Simultaneously short-range interactions via particle hopping promote superfluid order, such that a variety of structural and quantum coherent phases can occur. We analyze the emerging phase diagram in two dimensions by means of an extended Bose-Hubbard model using a local mean field approach combined with a superfluid cluster analysis. For commensurate ratios of the cavity and external lattice wave lengths the Mott insulator-superfluid transition is modified by the appearance of charge density wave and supersolid phases, at which the atomic density supports the buildup of a cavity field. For incommensurate ratios, the optomechanical forces induce the formation of Bose-glass and superglass phases, namely non-superfluid and superfluid phases, respectively, displaying quasi-periodic density modulations, which in addition can exhibit structural and superfluid stripe formation. The onset of such structures is constrained by the onsite interaction and is favourable at fractional densities. Experimental observables are identified and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Persistent-random-walk approach to anomalous transport of self-propelled particles

The motion of self-propelled particles is modeled as a persistent random walk. An analytical framework is developed that allows the derivation of exact expressions for the time evolution of arbitrary moments of the persistent walk's displacement. It is shown that the interplay of step length and turning angle distributions and self-propulsion produces various signs of anomalous diffusion at short time scales and asymptotically a normal diffusion behavior with a broad range of diffusion coefficients. The crossover from the anomalous short time behavior to the asymptotic diffusion regime is studied and the parameter dependencies of the crossover time are discussed. Higher moments of the displacement distribution are calculated and analytical expressions for the time evolution of the skewness and the kurtosis of the distribution are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anomalous Diffusion of Self-Propelled Particles in Directed Random Environments

We theoretically study the transport properties of self-propelled particles on complex structures, such as motor proteins on filament networks. A general master equation formalism is developed to investigate the persistent motion of individual random walkers, which enables us to identify the contributions of key parameters: the motor processivity, and the anisotropy and heterogeneity of the underlying network. We prove the existence of different dynamical regimes of anomalous motion, and that the crossover times between these regimes as well as the asymptotic diffusion coefficient can be increased by several orders of magnitude within biologically relevant control parameter ranges. In terms of motion in continuous space, the interplay between stepping strategy and persistency of the walker is established as a source of anomalous diffusion at short and intermediate time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice Model for Spontaneous Imbibition in Porous Media: The Role of Effective Tension and Universality Class

Recently anomalous scaling properties of front broadening during spontaneous imbibition of water in Vycor glass, a nano-porous medium, were reported: the mean height and the width of the propagating front increase with time $t$ both proportional to $t^{1/2}$. Here we propose a simple lattice imbibition model and elucidate quantitatively how the correlation range of the hydrostatic pressure and the disorder strength of the pore radii affect the scaling properties of the imbibition front. We introduce an effective tension of liquid across neighboring pores, which depends on the aspect ratio of each pore, and show that it leads to a dynamical crossover: both the mean height and the roughness grow faster in the presence of tension in the intermediate-time regime but eventually saturate in the long-time regime. The universality class of the long-time behavior is discussed by examining the associated scaling exponents and their relation to directed percolation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-equilibrium quantum relaxation across a localization-delocalization transition

We consider the one-dimensional $XX$-model in a quasi-periodic transverse-field described by the Harper potential, which is equivalent to a tight-binding model of spinless fermions with a quasi-periodic chemical potential. For weak transverse field (chemical potential), $h<h_c$, the excitations (fermions) are delocalized, but become localized for $h>h_c$. We study the non-equilibrium relaxation of the system by applying two protocols: a sudden change of $h$ (quench dynamics) and a slow change of $h$ in time (adiabatic dynamics). For a quench into the delocalized (localized) phase, the entanglement entropy grows linearly (saturates) and the order parameter decreases exponentially (has a finite limiting value). For a critical quench the entropy increases algebraically with time, whereas the order parameter decreases with a stretched-exponential. The density of defects after an adiabatic field change through the critical point is shown to scale with a power of the rate of field change and a scaling relation for the exponent is derived.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The quantum phase transition and correlations in the multi-spin-boson model

We consider multiple non-interacting quantum mechanical two-level systems coupled to a common bosonic bath and study its quantum phase transition with Monte Carlo simulations using a continuous imaginary time cluster algorithm. The common bath induces an effective ferromagnetic interaction between the otherwise independent two-level systems, which can be quantified by an effective interaction strength. For degenerate energy levels above a critical value of the bath coupling strength $α$ all two-level systems freeze into the same state and the critical value $α_c$ decreases asymptotically as $1/N$ with increasing $N$. For a finite number, $N$, of two-level systems the quantum phase transition (at zero temperature) is in the same universality class as the single spin-boson model, in the limit $N\to\infty$ the system shows mean-field critical behavior independent of the power of the spectral function of the bosonic bath. We also study the influence of a spatial separation of the spins in a bath of bosonic modes with linear dispersion relation on the location and characteristics of the phase transition as well as on correlations between the two-level systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum phases of incommensurate optical lattices due to cavity backaction

Ultracold bosonic atoms are confined by an optical lattice inside an optical resonator and interact with a cavity mode, whose wave length is incommensurate with the spatial periodicity of the confining potential. We predict that the intracavity photon number can be significantly different from zero when the atoms are driven by a transverse laser whose intensity exceeds a threshold value and whose frequency is suitably detuned from the cavity and the atomic transition frequency. In this parameter regime the atoms form clusters in which they emit in phase into the cavity. The clusters are phase locked, thereby maximizing the intracavity photon number. These predictions are based on a Bose-Hubbard model, whose derivation is here reported in detail. The Bose-Hubbard Hamiltonian has coefficients which are due to the cavity field and depend on the atomic density at all lattice sites. The corresponding phase diagram is evaluated using Quantum Monte Carlo simulations in one-dimension and mean-field calculations in two dimensions. Where the intracavity photon number is large, the ground state of the atomic gas lacks superfluidity and possesses finite compressibility, typical of a Bose-glass.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling theory for spontaneous imbibition in random networks of elongated pores

We present a scaling theory for the long time behavior of spontaneous imbibition in porous media consisting of interconnected pores with a large length-to-width ratio. At pore junctions the meniscus propagation in one or more branches can come to a halt when the Laplace pressure of the meniscus exceeds the hydrostatic pressure within the junction. We derive the scaling relations for the emerging arrest time distribution and show that the average front width is proportional to the height, yielding a roughness exponent of exactly beta=1/2 and explaining recent experimental results for nano-porous Vycor glass (NVG). Extensive simulations of a pore network model confirm these predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Superfluid clusters, percolation and phase transitions in the disordered, two dimensional Bose-Hubbard model

The Bose glass (BG) phase is the Griffiths region of the disordered Bose Hubbard model (BHM), characterized by finite, quasi-superfluid clusters within a Mott insulating background. We propose to utilize this characterization to identify the complete zero-temperature phase diagram of the disordered BHM in $d\ge2$ dimensions by analyzing the geometric properties of what we call superfluid (SF) clusters, which are defined to be clusters of sites with non-integer expectation values for the local boson occupation number. The Mott insulator (MI) phase then is the region in the phase diagram where no SF clusters exist, and the SF phase the region, where SF clusters percolate - the BG phase is in between: SF clusters exist, but do not percolate. This definition is particularly useful in the context of local mean field (LMF, or Gutzwiller-Ansatz) calculations, where we show that an identification of the phases on the basis of global quantities like the averaged SF order parameter and the compressibility are misleading. We apply the SF cluster analysis to the LMF ground states of the two dimensional disordered BHM to produce its phase diagram and find a) an excellent agreement with the phase diagram predicted on the basis of quantum Monte Carlo simulations for the commensurate density $n=1$, and b) large differences to stochastic mean field and other mean field predictions for fixed disorder strength. The relation of the percolation transition of the SF clusters with the onset of non-vanishing SF stiffness indicating the BG to SF transition is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anomalous Front Broadening During Spontaneous Imbibition in a Matrix with Elongated Pores

During spontaneous imbibition a wetting liquid is drawn into a porous medium by capillary forces. In systems with comparable pore length and diameter, such as paper and sand, the front of the propagating liquid forms a continuous interface. Sections of this interface advance in a highly correlated manner due to an effective surface tension, which restricts front broadening. Here we investigate water imbibition in a nanoporous glass (Vycor) in which the pores are much longer than they are wide. In this case, no continuous liquid-vapor interface with coalesced menisci can form. Anomalously fast imbibition front roughening is experimentally observed by neutron imaging.We propose a theoretical pore network model, whose structural details are adapted to the microscopic pore structure of Vycor glass, and show that it displays the same large scale roughening characteristics as observed in the experiment. The model predicts that menisci movements are uncorrelated. This indicates that despite the connectivity of the network the smoothening effect of surface tension on the imbibition front roughening is negligible. These results suggest a new universality class of imbibition behavior which is expected to occur in any matrix with elongated, interconnected pores of random radii.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bose-glass phases of ultracold atoms due to cavity backaction

We determine the quantum ground-state properties of ultracold bosonic atoms interacting with the mode of a high-finesse resonator. The atoms are confined by an external optical lattice, whose period is incommensurate with the cavity mode wave length, and are driven by a transverse laser, which is resonant with the cavity mode. While for pointlike atoms photon scattering into the cavity is suppressed, for sufficiently strong lasers quantum fluctuations can support the build-up of an intracavity field, which in turn amplifies quantum fluctuations. The dynamics is described by a Bose-Hubbard model where the coefficients due to the cavity field depend on the atomic density at all lattice sites. Quantum Monte Carlo simulations and mean-field calculations show that for large parameter regions cavity backaction forces the atoms into clusters with a checkerboard density distribution. Here, the ground state lacks superfluidity and possesses finite compressibility, typical of a Bose-glass. This system constitutes a novel setting where quantum fluctuations give rise to effects usually associated with disorder.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient kinetic Monte Carlo method for reaction-diffusion processes with spatially varying annihilation rates

We present an efficient Monte Carlo method to simulate reaction-diffusion processes with spatially varying particle annihilation or transformation rates as it occurs for instance in the context of motor-driven intracellular transport. Like Green's function reaction dynamics and first-passage time methods, our algorithm avoids small diffusive hops by propagating sufficiently distant particles in large hops to the boundaries of protective domains. Since for spatially varying annihilation or transformation rates the single particle diffusion propagator is not known analytically, we present an algorithm that generates efficiently either particle displacements or annihilations with the correct statistics, as we prove rigorously. The numerical efficiency of the algorithm is demonstrated with an illustrative example.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-equilibrium quantum dynamics after local quenches

We study the quantum dynamics resulting from preparing a one-dimensional quantum system in the ground state of initially two decoupled parts which are then joined together (local quench). Specifically we focus on the transverse Ising chain and compute the time-dependence of the magnetization profile, m_l(t), and correlation functions at the critical point, in the ferromagnetically ordered phase and in the paramagnetic phase. At the critical point we find finite size scaling forms for the nonequilibrium magnetization and compare predictions of conformal field theory with our numerical results. In the ferromagnetic phase the magnetization profiles are well matched by our predictions from a quasi-classical calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum relaxation after a quench in systems with boundaries

We study the time-dependence of the magnetization profile, m_l(t), of a large finite open quantum Ising chain after a quench. We observe a cyclic variation, in which starting with an exponentially decreasing period the local magnetization arrives to a quasi-stationary regime, which is followed by an exponentially fast reconstruction period. The non-thermal behavior observed at near-surface sites turns over to thermal behavior for bulk sites. Besides the standard time- and length-scales a non-standard time-scale is identified in the reconstruction period.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Semi-classical theory for quantum quenches in finite transverse Ising chains

We present a quantitative semi-classical theory for the non-equilibrium dynamics of transverse Ising chains after quantum quenches, in particular sudden changes of the transverse field strength. We obtain accurate predictions for the quench dependent relaxation times and correlation lengths, and also about the recurrence times and quasiperiodicity of time dependent correlations in finite systems with open or periodic boundary conditions. We compare the quantiative predictions of our semi-classical theory (local magnetization, equal time bulk-bulk and surface-to-bulk correlations, and bulk autocorrelations) with the results from exact free fermion calculations and discuss the range of applicability of the semi-classical theory and possible generaliztaions and extensions.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The quantum phase transition in the sub-ohmic spin-boson model: Quantum Monte-Carlo study with a continuous imaginary time cluster algorithm

A continuous time cluster algorithm for two-level systems coupled to a dissipative bosonic bath is presented and applied to the sub-ohmic spin-Boson model. When the power s of the spectral function J(w) \propto w^s is smaller than 1/2, the critical exponents are found to be classical, mean-field like. Potential sources for the discrepancy with recent renormalization group predictions are traced back to the effect of a dangerously irrelevant variable.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Path integral Monte Carlo study of the interacting quantum double-well model: Quantum phase transition and phase diagram

The discrete time path integral Monte Carlo (PIMC) with a one-particle density matrix approximation is applied to study the quantum phase transition in the coupled double-well chain. To improve the convergence properties, the exact action for a single particle in a double well potential is used to construct the many-particle action. The algorithm is applied to the interacting quantum double-well chain for which the zero-temperature phase diagram is determined. The quantum phase transition is studied via finite-size scaling and the critical exponents are shown to be compatible with the classical two-dimensional (2D) Ising universality class -- not only in the order-disorder limit (deep potential wells) but also in the displacive regime (shallow potential wells).

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Super-Aging in two-dimensional random ferromagnets

We study the aging properties, in particular the two-time autocorrelations, of the two-dimensional randomly diluted Ising ferromagnet below the critical temperature via Monte-Carlo simulations. We find that the autocorrelation function displays additive aging $C(t,t_w)=C_{st}(t)+C_{ag}(t,t_w)$, where the stationary part $C_{st}$ decays algebraically. The aging part shows anomalous scaling $C_{ag}(t,t_w)={\cal C}(h(t)/h(t_w))$, where $h(u)$ is a non-homogeneous function excluding a $t/t_w$ scaling.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Aging in Disordered Systems

The dynamics of strongly disordered systems becomes extremely slow or glassy at low temperatures, which results in a characteristic aging scenario. This means that the outcome of measurements strongly depends on the history of the system within the glassy phase, even on macroscopic time scale like hours, weeks or years. This area of non-equilibrium dynamics in disordered systems became recently a major focus of research interest, in particular with respect to spin glasses and related systems. Here we give an overview on these activities with a focus on Monte Carlo studies.

0 citations0 reviewsNo code linkedOpen access