Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Near-critical confined fluids and Ising films: Density-matrix renormalization-group study

Anna Maciołek, A. Drzewiński, Robert Evans

Published Oct 30, 2001• 25 citations
Publisher

Abstract

Two-dimensional Ising strips subject to identical surface fields h(1)=h(2) > or =0 are studied for temperatures above and below the bulk critical temperature T(c) and a range of bulk fields h by means of the density-matrix renormalization-group method. In the case of nonvanishing surface fields, the near-critical behavior of the solvation force f(solv), total adsorption Gamma, inverse longitudinal correlation length xi(parallel)( -1) and specific heat C(H) is strongly influenced by the (pseudo) capillary condensation that occurs below T(c). We obtain scaling functions of f(solv), Gamma, and xi(parallel)(-1). C(H) exhibits a weakly rounded singularity on crossing the pseudocoexistence line. We contrast these results with those for the case of free boundaries where, for temperatures slightly below T(c), f(solv) and C(H) exhibit a sharp extremum away from h=0. Our results have direct repercussions for the properties of near-critical Ising films in three dimensions and we argue that the long-ranged solvation (Casimir) force in confined fluids should be more attractive in the neighborhood of the capillary critical point than exactly at the bulk critical point.

DOI
10.1103/physreve.64.056137

Related papers

Effect of bulk magnetic field on critical Ising films

A. Drzewiński, Anna Maciołek, A. Ciach

Two-dimensional Ising films $L\ifmmode\times\else\texttimes\fi{}\ensuremath{\infty}$ in a nonvanishing bulk magnetic field H are studied at the bulk critical temperature ${T}_{c}$ for two choices of surface fields (a) ${H}_{1}{=H}_{L}=0$ (ordinary transition), and (b) ${H}_{1}{=H}_{L}=\ensuremath{\infty}$ (normal transition) by the density-matrix renormalization-group method. Universal scaling functions for magnetization profiles, the excess magnetization $\ensuremath{\Gamma},$ the longitudinal correlation length ${\ensuremath{\xi}}_{\ensuremath{\Vert}},$ and for the analog of the solvation force ${f}_{\mathrm{solv}}$ are found and discussed. When ${H}_{1}=0$ the scaling function for ${f}_{\mathrm{solv}}$ has two symmetric minima at $y=\mathrm{sgn}(H)L|H{|}^{\ensuremath{\nu}/\ensuremath{\Delta}}\ensuremath{\approx}\ifmmode\pm\else\textpm\fi{}1$ with an amplitude at the minimum about 3.8 times the value at $H=0,$ the Casimir amplitude. For the normal transition the scaling function for ${f}_{\mathrm{solv}}$ has a single minimum near the continuation of the pseudocoexistence (capillary condensation) line, with an amplitude about 100 times the Casimir amplitude.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•May 1, 2000•19 citations

Solvation forces in Ising films with long-range boundary fields: density-matrix renormalization-group study

A. Drzewiński, Anna Maciołek, Artur Barasiński

Using the quasi-exact density-matrix renormalization-group method we calculate the solvation forces in two-dimensional Ising films of thickness L subject to identical algebraically decaying boundary fields with various decay exponents p. At the bulk critical point the solvation force acquires a universal contribution which is long-ranged in L due to the critical fluctuations, a phenomenon known as the critical Casimir effect. For p = 2, 3 and 50, we study the scaling behaviour of the solvation force along the pseudo-phase coexistence and along the critical and sub-critical isotherms.

arXiv:1102.1203•Molecular Physics•Mar 30, 2011•4 citations

Interplay of complete wetting, critical adsorption, and capillary condensation

A. Drzewiński, Anna Maciołek, Artur Barasiński +1 more

The excess adsorption Gamma in two-dimensional Ising strips (infinityxL), subject to identical boundary fields at both one-dimensional surfaces decaying in the orthogonal direction j as -h1j(-p), is studied for various values of p and along various thermodynamic paths below the bulk critical point by means of the density-matrix renormalization-group method. The crossover behavior between the complete-wetting and critical-adsorption regimes, occurring in semi-infinite systems, is strongly influenced by confinement effects. Along isotherms T=const the asymptotic power-law dependences on the external bulk field, which characterize these two regimes, are pre-empted by capillary condensation. Along the pseudo-first-order phase-coexistence line of the strips, which varies with temperature, we find a broad crossover regime in which both the thickness of the wetting film and Gamma increase as functions of the reduced temperature tau but do not follow any power law. Above the wetting temperature the order-parameter profiles are not slablike but exhibit wide interfacial variations and pronounced tails.

arXiv:0808.2966•Physical Review E•Apr 28, 2009•28 citations

Critical fluctuations of a confined binary mixture

Małgorzata Zubaszewska, Andrej Gendiar, Tomasz Masłowski

Exploiting the mapping between a binary mixture and the Ising model we have analyzed the critical fluctuations by means of the density-matrix renormalization group technique. The calculations have been carried out for a two-dimensional Ising strip subject to equal strong surface fields. It was found that the critical Casimir force displays significantly different behavior on opposite sides of the capillary condensation line, especially below the critical temperature. It can be concluded that in a real binary mixtures the most attractive force appears at temperatures near T_C and at reservoir compositions slightly away from the critical composition.

arXiv:1410.2015•Oct 8, 2014•0 citations

Crossover between ordinary and normal transitions in the presence of a bulk field

Anna Maciołek, A. Drzewiński, A. Ciach

We investigate two-dimensional Ising films at the critical temperature T(c) and nonzero bulk magnetic field h using the density-matrix renormalization-group method. The crossover between ordinary (h(1)=0) and normal (h(1)=infinity) transitions corresponding to finite values of the surface fields h(1)=h(2), is studied. The structure and the solvation force f(solv) as a function of h, crucially depend on the value of h(1). Scaling functions for f(solv) and the longitudinal correlation length are given and discussed.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jul 25, 2001•6 citations

Crossover between ordinary and normal transitions in two dimensional critical Ising films

Anna Maciołek, A. Ciach, A. Drzewiński

We investigate two dimensional critical Ising films of width L with surface fields H(1)=H(L) in the crossover between ordinary (H(1)=0) and normal (H(1)=infinity) transitions. Using exact transfer-matrix diagonalization and density matrix renormalization-group (DMRG) methods, we calculate magnetization profiles m(z), the excess magnetization Gamma, and the analog of the solvation force f(solv) as functions of H1 for several L. Scaling functions of the above quantities deviate substantially from their asymptotic forms at fixed points for a broad region of the scaling variable LH21 approximately L/l(1), where l(1) is the length induced by the surface field H1. The scaling function for /f(solv)/ has a deep minimum near LH(2)(1)=1, which is about one order of magnitude smaller than its value at both fixed points (the "Casimir" amplitude). For weak H1 (l(1)>L) the magnetization profile has a maximum at the center of the film, and f(solv) decays much faster than L-2. For stronger H1 (1<l(1)<L), the magnetization has two maxima at a distance approximately l(1) from the walls, and the solvation force decays much slower than L-2. For L>>l(1) the solvation force decays according to the universal power law f(solv) approximately L(-2). The results of the approximate DMRG method show remarkable agreement with the exact ones.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Sep 1, 1999•28 citations
Two-dimensional Ising strips subject to identical surface fields h(1)=h(2) > or =0 are studied for temperatures above and below the bulk critical temperature T(c) and a range of bulk fields h by means of the density-matrix renormalization-group method. In the case of nonvanishing surface fields, the near-critical behavior of the solvation force f(solv), total adsorption Gamma, inverse longitudinal correlation length xi(parallel)( -1) and specific heat C(H) is strongly influenced by the (pseudo) capillary condensation that occurs below T(c). We obtain scaling functions of f(solv), Gamma, and xi(parallel)(-1). C(H) exhibits a weakly rounded singularity on crossing the pseudocoexistence line. We contrast these results with those for the case of free boundaries where, for temperatures slightly below T(c), f(solv) and C(H) exhibit a sharp extremum away from h=0. Our results have direct repercussions for the properties of near-critical Ising films in three dimensions and we argue that the long-ranged solvation (Casimir) force in confined fluids should be more attractive in the neighborhood of the capillary critical point than exactly at the bulk critical point.
Two-dimensional Ising films $L\ifmmode\times\else\texttimes\fi{}\ensuremath{\infty}$ in a nonvanishing bulk magnetic field H are studied at the bulk critical temperature Tc{T}_{c}Tc​ for two choices of surface fields (a) H1=HL=0{H}_{1}{=H}_{L}=0H1​=HL​=0 (ordinary transition), and (b) ${H}_{1}{=H}_{L}=\ensuremath{\infty}$ (normal transition) by the density-matrix renormalization-group method. Universal scaling functions for magnetization profiles, the excess magnetization $\ensuremath{\Gamma},$ the longitudinal correlation length ${\ensuremath{\xi}}_{\ensuremath{\Vert}},$ and for the analog of the solvation force fsolv{f}_{\mathrm{solv}}fsolv​ are found and discussed. When H1=0{H}_{1}=0H1​=0 the scaling function for fsolv{f}_{\mathrm{solv}}fsolv​ has two symmetric minima at $y=\mathrm{sgn}(H)L|H{|}^{\ensuremath{\nu}/\ensuremath{\Delta}}\ensuremath{\approx}\ifmmode\pm\else\textpm\fi{}1$ with an amplitude at the minimum about 3.8 times the value at H=0,H=0,H=0, the Casimir amplitude. For the normal transition the scaling function for fsolv{f}_{\mathrm{solv}}fsolv​ has a single minimum near the continuation of the pseudocoexistence (capillary condensation) line, with an amplitude about 100 times the Casimir amplitude.
Using the quasi-exact density-matrix renormalization-group method we calculate the solvation forces in two-dimensional Ising films of thickness L subject to identical algebraically decaying boundary fields with various decay exponents p. At the bulk critical point the solvation force acquires a universal contribution which is long-ranged in L due to the critical fluctuations, a phenomenon known as the critical Casimir effect. For p = 2, 3 and 50, we study the scaling behaviour of the solvation force along the pseudo-phase coexistence and along the critical and sub-critical isotherms.
The excess adsorption Gamma in two-dimensional Ising strips (infinityxL), subject to identical boundary fields at both one-dimensional surfaces decaying in the orthogonal direction j as -h1j(-p), is studied for various values of p and along various thermodynamic paths below the bulk critical point by means of the density-matrix renormalization-group method. The crossover behavior between the complete-wetting and critical-adsorption regimes, occurring in semi-infinite systems, is strongly influenced by confinement effects. Along isotherms T=const the asymptotic power-law dependences on the external bulk field, which characterize these two regimes, are pre-empted by capillary condensation. Along the pseudo-first-order phase-coexistence line of the strips, which varies with temperature, we find a broad crossover regime in which both the thickness of the wetting film and Gamma increase as functions of the reduced temperature tau but do not follow any power law. Above the wetting temperature the order-parameter profiles are not slablike but exhibit wide interfacial variations and pronounced tails.
Exploiting the mapping between a binary mixture and the Ising model we have analyzed the critical fluctuations by means of the density-matrix renormalization group technique. The calculations have been carried out for a two-dimensional Ising strip subject to equal strong surface fields. It was found that the critical Casimir force displays significantly different behavior on opposite sides of the capillary condensation line, especially below the critical temperature. It can be concluded that in a real binary mixtures the most attractive force appears at temperatures near T_C and at reservoir compositions slightly away from the critical composition.
We investigate two-dimensional Ising films at the critical temperature T(c) and nonzero bulk magnetic field h using the density-matrix renormalization-group method. The crossover between ordinary (h(1)=0) and normal (h(1)=infinity) transitions corresponding to finite values of the surface fields h(1)=h(2), is studied. The structure and the solvation force f(solv) as a function of h, crucially depend on the value of h(1). Scaling functions for f(solv) and the longitudinal correlation length are given and discussed.
We investigate two dimensional critical Ising films of width L with surface fields H(1)=H(L) in the crossover between ordinary (H(1)=0) and normal (H(1)=infinity) transitions. Using exact transfer-matrix diagonalization and density matrix renormalization-group (DMRG) methods, we calculate magnetization profiles m(z), the excess magnetization Gamma, and the analog of the solvation force f(solv) as functions of H1 for several L. Scaling functions of the above quantities deviate substantially from their asymptotic forms at fixed points for a broad region of the scaling variable LH21 approximately L/l(1), where l(1) is the length induced by the surface field H1. The scaling function for /f(solv)/ has a deep minimum near LH(2)(1)=1, which is about one order of magnitude smaller than its value at both fixed points (the "Casimir" amplitude). For weak H1 (l(1)>L) the magnetization profile has a maximum at the center of the film, and f(solv) decays much faster than L-2. For stronger H1 (1<l(1)<L), the magnetization has two maxima at a distance approximately l(1) from the walls, and the solvation force decays much slower than L-2. For L>>l(1) the solvation force decays according to the universal power law f(solv) approximately L(-2). The results of the approximate DMRG method show remarkable agreement with the exact ones.