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 A

Exact and approximate differential renormalization-group generators. II. The equation of state

Jeffrey F. Nicoll, T. S. Chang

Published Jun 1, 1978• 35 citations
Publisher

Abstract

We apply an approximate form of the exact one-particle-irreducible renormalization-group generator to the calculation of the equation of state. Several approaches are explored. (i) Global nonlinear trajectories for the Hamiltonian parameters are exploited to give nonlinear crossover equations of state. Logarithms of the reduced temperature $t\ensuremath{\equiv}\frac{(T\ensuremath{-}{T}_{c})}{{T}_{c}}$ are automatically exponentiated to give power-law behavior. (ii) The form and asymptotic properties of the equation of state are described for those systems whose complete nonlinear trajectories cannot be explicitly obtained. (iii) Nonlinear trajectories for the irreducible Green's functions are solved to give a fully exponentiated equation of state containing no logarithmic terms. (iv) Simple critical (noncrossover) equations of state are obtained by iteration of the generator with linear and quasilinear trajectories. (v) Operator nonlinear equations are solved to give the crossover equation of state for an arbitrary order $\ensuremath{\theta}$ critical system, both in an $\ensuremath{\epsilon}$ expansion and at the borderline dimension. (vi) Combination of the Green's functions are formed into Green's eigenfunctions or operators and their nonlinear trajectories used to calculate fully exponentiated equations of state for order $\ensuremath{\theta}$ Ising systems. Systems considered include the usual Wilson-Fisher Hamiltonian [methods (i)-(iv)], the Sak-model compressible ferromagnet [(i)-(ii)], the $\mathrm{nm}$-hypercubical model (of which the dilute quenched random ferromagnet forms the $n\ensuremath{\rightarrow}0$ limit) [(ii)], and the isotropic $n$-component order $\ensuremath{\theta}$ ferromagnet [(iv)-(vi)]. The tricritical ($\ensuremath{\theta}=3$) case is given explicitly as an example of the general formalism of methods (iv)-(vi). Throughout, we use a bare critical propagator which is an arbitrary generalized homogeneous function of the components of the wave vector. This allows us to simultaneously describe ordinary critical systems and anisotropic Lifschitz points as well as certain structural and spin-reorientation phase transitions.

DOI
10.1103/physreva.17.2083

Related papers

Exact and approximate differential renormalization-group generators

Jeffrey F. Nicoll, T. S. Chang, Eugene Stanley

In dealing with critical phenomena of complex systems that simulate realistic materials, the full structure of the renormalization group is often unnecessarily cumbersome. For approximate calculations and for systems with special properties, specialized generators are simpler to apply. We derive several such exact and approximate differential generators and solve a number of interesting practical problems to illustrate this approach: (i) We derive a new approximate differential generator based on the Wilson incomplete-integration generator. Using this generator we calculate for an $n$-component spin system the eigenvalues (critical-point exponents) associated with a critical point of arbitrary order $\mathcal{O}$ and "propagator exponent" $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}$ to first order in the expansion parameter ${\ensuremath{\epsilon}}_{\mathcal{O}}(\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}})=d+\mathcal{O}(\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}\ensuremath{-}d)$; this extends previous work for $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=2$, $\mathcal{O}$ arbitrary; $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}\ensuremath{\le}2$, $\mathcal{O}=2$ (long-rang forces); and $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=4$, $\mathcal{O}=2$ (the "Lifshitz point"). Our results agree with those obtainable using an approximate generator based on the Wegner-Houghton equation. The cases $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=2L$ ($L$ a positive integer \ensuremath{\ge}2) describe the onset of helical ordering for which $|\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}|\ensuremath{\sim}{(\ensuremath{-}p)}^{{\ensuremath{\beta}}_{k}}$, where ${\ensuremath{\beta}}_{k}=\frac{1}{2(L\ensuremath{-}1)}+o({\ensuremath{\epsilon}}_{\mathcal{O}}^{2}(2L))$ and $p$ parametrizes the hypersurface of critical points. For $p>0$, the ordered phase is uniform; for $p<0$ there is spiral order. The point $p=0$, at which such nonuniform ordering commences, we term a generalized Lifshitz point of Lifschitz character $L$. (ii) We consider the full Wilson and Wegner-Houghton generators in the paired spin-momenta limit and the $n\ensuremath{\rightarrow}\ensuremath{\infty}$ limit for even-order critical points. These limiting generators are identical for both full generators. This demonstrates that at least in these cases the Wilson and Wegner-Houghton generators agree exactly, without recourse to perturbation theory. These simple exact generators should provide "anchors" for calculations of exponents for higher-order critical points. (iii) We derive approximate generators which are suitable for compressible magnetic systems and more general systems with constraints for which the spin momenta are grouped in any arbitrary manner. We apply this to the case of a simple compressible magnet model and obtain the exact renormalization-group trajectories to order $\ensuremath{\epsilon}$ with $\ensuremath{\epsilon}\ensuremath{\equiv}{\ensuremath{\epsilon}}_{2}(2)=4\ensuremath{-}d$.

Physical Review A•Mar 1, 1976•78 citations

Nonlinear Crossover between Critical and Tricritical Behavior

Jeffrey F. Nicoll, T. S. Chang, Eugene Stanley

Recently approximate solutions to renormalization-group equations were used to calculate the double-power "tricritical" scaling equation of state. We show that this equation can be simply calculated using our exact renormalization-group trajectories. The equation of state is expressed explicitly in terms of nonlinear scaling fields.

Physical Review Letters•Jan 12, 1976•67 citations

Q -state Potts model by Wilson's exact renormalization-group equation

Kathie E. Neuman, Eberhard K. Riedel, Shunichi Muto

Critical properties of the $Q$-state Potts model for dimensions $3\ensuremath{\le}d\ensuremath{\le}6$ are calculated by means of Wilson's exact momentum-space renormalization-group equation. The scaling-field method of Golner and Riedel is used to approximate the functional differential equation by a set of 11 ordinary coupled differential equations. For $d=4\ensuremath{-}\ensuremath{\epsilon}$, lines of critical and tricritical Potts fixed points are found as functions of $Q$ that annihilate as $Q$ approaches a critical value ${Q}_{c}=\frac{2+{\ensuremath{\epsilon}}^{2}}{a+O({\ensuremath{\epsilon}}^{3})}$. For $Q>{Q}_{c}$, the Potts transition is first order. Along these fixed lines the critical and tricritical exponents (upper and lower sign, respectively) are to leading order: $\frac{1}{\ensuremath{\nu}}=2\ensuremath{-}\frac{1}{6}[\ensuremath{\epsilon}\ifmmode\pm\else\textpm\fi{}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}]$, $\frac{\ensuremath{\varphi}}{\ensuremath{\nu}}=\ensuremath{\mp}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}$, and $\ensuremath{\eta}=\frac{{[\ensuremath{\epsilon}\ifmmode\pm\else\textpm\fi{}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}]}^{2}}{216}+b\ensuremath{\delta}$, where $\ensuremath{\epsilon}=4\ensuremath{-}d$, $\ensuremath{\delta}=Q\ensuremath{-}2$, and $\ensuremath{\delta}\ensuremath{\le}{\ensuremath{\delta}}_{c}=\frac{{\ensuremath{\epsilon}}^{2}}{a}+O({\ensuremath{\epsilon}}^{3})$. While the form of the $\ensuremath{\epsilon}$ and $\ensuremath{\delta}$ dependences is exact, the coefficients $a$ and $b$ cannot be obtained systematically by $\ensuremath{\epsilon}$ expansion, since the upper critical dimensionality of the Potts model is six when $Q\ensuremath{\ne}2$. In our truncation, $a=6.52$ and $b=0.065$. The results have been extended to dimensions $3.4\ensuremath{\le}d\ensuremath{\lesssim}4$ by solving the renormalization-group equations numerically. The percolation limit of the Potts model, $Q=1$, is also investigated and the critical exponents ${\ensuremath{\nu}}^{P},{\ensuremath{\varphi}}^{P}$, and ${\ensuremath{\eta}}^{P}$ determined as functions of dimension for $3\ensuremath{\le}d\ensuremath{\le}6$.

Physical review. B, Condensed matter•Jan 1, 1984•33 citations

Renormalization and Statistical Mechanics in Many-Particle Systems. II. Statistical Perturbation Method

Duane C. Wallace

A general nonperturbation problem, described by the Hamiltonian $\mathcal{H}$, is considered. Approximate creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are defined as operators which satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}+{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where ${\ensuremath{\omega}}_{i}$ are the creation energies, and the remainder operators ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small in the sense that statistical averages involving the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small. A zeroth approximation is to neglect the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$. In this case, a straight-forward derivation provides general relations between statistical averages such as $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}〉$ and $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}\ifmmode\pm\else\textpm\fi{}\ensuremath{\Omega}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}〉$, where $\ensuremath{\Omega}$ is any operator. The usual boson and fermion occupation numbers are a trivial result of these zeroth-order relations. Corrections to the general zeroth-order relations arise from the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$; such corrections are considered in the spirit of perturbation theory. An explicit first-order correction is obtained for the averages $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\theta}}_{i}〉$, and this correction is related to a self-consistent energy-renormalization procedure. The present method is similar to the method of thermodynamic Green's functions in that the derivations given here do not require evaluation of any state vectors or of the partition function. The general results are applied to the Heisenberg ferromagnet problem. At low temperatures, the zeroth-order relations give the Bloch spin-wave results, and an approximate evaluation of the first-order contribution gives the leading term in Dyson's ${T}^{4}$ correction to the spontaneous magnetization. For arbitrary temperatures, the zeroth-order relations give the Tyablikov equations for $〈{S}^{z}〉$, and Callen's results are obtained from the first-order corrections. These examples illustrate the simplicity of calculating statistical averages of functions of the ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, ${\ensuremath{\theta}}_{i}$ operators by the present method.

Physical Review•Dec 2, 1966•21 citations

Renormalization group for first-order phase transitions: Equation of state of the compressible Ising magnet

J. Bruno, Josef Sak

Thermodynamic functions of an $n$-component magnetic system coupled to an isotropic elastic medium in $d=4\ensuremath{-}\ensuremath{\epsilon}$ dimensions are calculated to leading order in $\ensuremath{\epsilon}$. Solutions to renormalization-group recursion relations correct to $O(\ensuremath{\epsilon})$ along with trajectory-integral methods are employed in these calculations. For $n>~2$ the finite compressibility modifies the critical-point behavior producing new corrections to scaling in the thermodynamic functions which are given. In the $n=1$ case renormalization-group predictions of a first-order transition are verified. Thermodynamic functions are given and analyzed in this case. In particular the transition temperature and discontinuity in order parameter are determined and the phase diagram is constructed.

Physical review. B, Condensed matter•Oct 1, 1980•42 citations

Application of the renormalization group to phase transitions in disordered systems

G. Grinstein, Alan H. Luther

The critical behavior of spin systems with quenched disorder is studied by renormalization-group methods. For the randomly dilute $m$-vector model, the $n=0$ limit is used to construct a translationally invariant effective Hamiltonian which describes the original disordered system. This Hamiltonian is analyzed in the $\ensuremath{\epsilon}$ expansion to order ${\ensuremath{\epsilon}}^{2}$. Sharp second-order phase transitions with exponents which do not depend continuously on impurity concentration are predicted. For $m>{m}_{c}\ensuremath{\equiv}4\ensuremath{-}4\ensuremath{\epsilon}+O({\ensuremath{\epsilon}}^{2})$ the isotropic $m$-component fixed point, which characterizes the critical behavior of the pure system, is stable. For $m<{m}_{c}$, a new random fixed point becomes stable. The exponents corresponding to this fixed point are $\ensuremath{\eta}=[\frac{(5{m}^{2}\ensuremath{-}8m)}{256{(m\ensuremath{-}1)}^{2}}]{\ensuremath{\epsilon}}^{2}+O({\ensuremath{\epsilon}}^{3})$, $\ensuremath{\nu}=\frac{1}{2}+[\frac{3m}{32(m\ensuremath{-}1)}]\ensuremath{\epsilon}+[\frac{m(127{m}^{2}\ensuremath{-}572m\ensuremath{-}32)}{4096{(m\ensuremath{-}1)}^{3}}]{\ensuremath{\epsilon}}^{2}+O({\ensuremath{\epsilon}}^{3})$ for $m\ensuremath{\ne}1$, and $\ensuremath{\eta}=\ensuremath{-}\frac{\ensuremath{\epsilon}}{106}+O({\ensuremath{\epsilon}}^{\frac{3}{2}})$, $\ensuremath{\nu}=\frac{1}{2}+\frac{{(\frac{6\ensuremath{\epsilon}}{53})}^{\frac{1}{2}}}{4}+O(\ensuremath{\epsilon})$ for $m=1$. More general random systems are qualitatively discussed from the effective-Hamiltonian viewpoint.

Physical review. B, Solid state•Feb 1, 1976•529 citations
We apply an approximate form of the exact one-particle-irreducible renormalization-group generator to the calculation of the equation of state. Several approaches are explored. (i) Global nonlinear trajectories for the Hamiltonian parameters are exploited to give nonlinear crossover equations of state. Logarithms of the reduced temperature t\ensuremath≡(T\ensuremath−Tc)Tct\ensuremath{\equiv}\frac{(T\ensuremath{-}{T}_{c})}{{T}_{c}}t\ensuremath≡Tc​(T\ensuremath−Tc​)​ are automatically exponentiated to give power-law behavior. (ii) The form and asymptotic properties of the equation of state are described for those systems whose complete nonlinear trajectories cannot be explicitly obtained. (iii) Nonlinear trajectories for the irreducible Green's functions are solved to give a fully exponentiated equation of state containing no logarithmic terms. (iv) Simple critical (noncrossover) equations of state are obtained by iteration of the generator with linear and quasilinear trajectories. (v) Operator nonlinear equations are solved to give the crossover equation of state for an arbitrary order \ensuremathθ\ensuremath{\theta}\ensuremathθ critical system, both in an \ensuremathϵ\ensuremath{\epsilon}\ensuremathϵ expansion and at the borderline dimension. (vi) Combination of the Green's functions are formed into Green's eigenfunctions or operators and their nonlinear trajectories used to calculate fully exponentiated equations of state for order \ensuremathθ\ensuremath{\theta}\ensuremathθ Ising systems. Systems considered include the usual Wilson-Fisher Hamiltonian [methods (i)-(iv)], the Sak-model compressible ferromagnet [(i)-(ii)], the nm\mathrm{nm}nm-hypercubical model (of which the dilute quenched random ferromagnet forms the n\ensuremath→0n\ensuremath{\rightarrow}0n\ensuremath→0 limit) [(ii)], and the isotropic nnn-component order \ensuremathθ\ensuremath{\theta}\ensuremathθ ferromagnet [(iv)-(vi)]. The tricritical (\ensuremathθ=3\ensuremath{\theta}=3\ensuremathθ=3) case is given explicitly as an example of the general formalism of methods (iv)-(vi). Throughout, we use a bare critical propagator which is an arbitrary generalized homogeneous function of the components of the wave vector. This allows us to simultaneously describe ordinary critical systems and anisotropic Lifschitz points as well as certain structural and spin-reorientation phase transitions.
In dealing with critical phenomena of complex systems that simulate realistic materials, the full structure of the renormalization group is often unnecessarily cumbersome. For approximate calculations and for systems with special properties, specialized generators are simpler to apply. We derive several such exact and approximate differential generators and solve a number of interesting practical problems to illustrate this approach: (i) We derive a new approximate differential generator based on the Wilson incomplete-integration generator. Using this generator we calculate for an nnn-component spin system the eigenvalues (critical-point exponents) associated with a critical point of arbitrary order O\mathcal{O}O and "propagator exponent" $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}$ to first order in the expansion parameter ${\ensuremath{\epsilon}}_{\mathcal{O}}(\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}})=d+\mathcal{O}(\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}\ensuremath{-}d)$; this extends previous work for $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=2$, O\mathcal{O}O arbitrary; $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}\ensuremath{\le}2$, O=2\mathcal{O}=2O=2 (long-rang forces); and $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=4$, O=2\mathcal{O}=2O=2 (the "Lifshitz point"). Our results agree with those obtainable using an approximate generator based on the Wegner-Houghton equation. The cases $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\sigma}}=2L$ (LLL a positive integer \ensuremath{\ge}2) describe the onset of helical ordering for which $|\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}|\ensuremath{\sim}{(\ensuremath{-}p)}^{{\ensuremath{\beta}}_{k}}$, where ${\ensuremath{\beta}}_{k}=\frac{1}{2(L\ensuremath{-}1)}+o({\ensuremath{\epsilon}}_{\mathcal{O}}^{2}(2L))$ and ppp parametrizes the hypersurface of critical points. For p>0p>0p>0, the ordered phase is uniform; for p<0p<0p<0 there is spiral order. The point p=0p=0p=0, at which such nonuniform ordering commences, we term a generalized Lifshitz point of Lifschitz character LLL. (ii) We consider the full Wilson and Wegner-Houghton generators in the paired spin-momenta limit and the $n\ensuremath{\rightarrow}\ensuremath{\infty}$ limit for even-order critical points. These limiting generators are identical for both full generators. This demonstrates that at least in these cases the Wilson and Wegner-Houghton generators agree exactly, without recourse to perturbation theory. These simple exact generators should provide "anchors" for calculations of exponents for higher-order critical points. (iii) We derive approximate generators which are suitable for compressible magnetic systems and more general systems with constraints for which the spin momenta are grouped in any arbitrary manner. We apply this to the case of a simple compressible magnet model and obtain the exact renormalization-group trajectories to order $\ensuremath{\epsilon}$ with $\ensuremath{\epsilon}\ensuremath{\equiv}{\ensuremath{\epsilon}}_{2}(2)=4\ensuremath{-}d$.
Critical properties of the QQQ-state Potts model for dimensions $3\ensuremath{\le}d\ensuremath{\le}6$ are calculated by means of Wilson's exact momentum-space renormalization-group equation. The scaling-field method of Golner and Riedel is used to approximate the functional differential equation by a set of 11 ordinary coupled differential equations. For $d=4\ensuremath{-}\ensuremath{\epsilon}$, lines of critical and tricritical Potts fixed points are found as functions of QQQ that annihilate as QQQ approaches a critical value ${Q}_{c}=\frac{2+{\ensuremath{\epsilon}}^{2}}{a+O({\ensuremath{\epsilon}}^{3})}$. For Q>QcQ>{Q}_{c}Q>Qc​, the Potts transition is first order. Along these fixed lines the critical and tricritical exponents (upper and lower sign, respectively) are to leading order: $\frac{1}{\ensuremath{\nu}}=2\ensuremath{-}\frac{1}{6}[\ensuremath{\epsilon}\ifmmode\pm\else\textpm\fi{}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}]$, $\frac{\ensuremath{\varphi}}{\ensuremath{\nu}}=\ensuremath{\mp}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}$, and $\ensuremath{\eta}=\frac{{[\ensuremath{\epsilon}\ifmmode\pm\else\textpm\fi{}{({\ensuremath{\epsilon}}^{2}\ensuremath{-}a\ensuremath{\delta})}^{\frac{1}{2}}]}^{2}}{216}+b\ensuremath{\delta}$, where $\ensuremath{\epsilon}=4\ensuremath{-}d$, $\ensuremath{\delta}=Q\ensuremath{-}2$, and $\ensuremath{\delta}\ensuremath{\le}{\ensuremath{\delta}}_{c}=\frac{{\ensuremath{\epsilon}}^{2}}{a}+O({\ensuremath{\epsilon}}^{3})$. While the form of the $\ensuremath{\epsilon}$ and $\ensuremath{\delta}$ dependences is exact, the coefficients aaa and bbb cannot be obtained systematically by $\ensuremath{\epsilon}$ expansion, since the upper critical dimensionality of the Potts model is six when $Q\ensuremath{\ne}2$. In our truncation, a=6.52a=6.52a=6.52 and b=0.065b=0.065b=0.065. The results have been extended to dimensions $3.4\ensuremath{\le}d\ensuremath{\lesssim}4$ by solving the renormalization-group equations numerically. The percolation limit of the Potts model, Q=1Q=1Q=1, is also investigated and the critical exponents ${\ensuremath{\nu}}^{P},{\ensuremath{\varphi}}^{P}$, and ${\ensuremath{\eta}}^{P}$ determined as functions of dimension for $3\ensuremath{\le}d\ensuremath{\le}6$.
A general nonperturbation problem, described by the Hamiltonian H\mathcal{H}H, is considered. Approximate creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are defined as operators which satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}+{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where ${\ensuremath{\omega}}_{i}$ are the creation energies, and the remainder operators ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small in the sense that statistical averages involving the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small. A zeroth approximation is to neglect the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$. In this case, a straight-forward derivation provides general relations between statistical averages such as $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}〉$ and $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}\ifmmode\pm\else\textpm\fi{}\ensuremath{\Omega}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}〉$, where $\ensuremath{\Omega}$ is any operator. The usual boson and fermion occupation numbers are a trivial result of these zeroth-order relations. Corrections to the general zeroth-order relations arise from the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$; such corrections are considered in the spirit of perturbation theory. An explicit first-order correction is obtained for the averages $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\theta}}_{i}〉$, and this correction is related to a self-consistent energy-renormalization procedure. The present method is similar to the method of thermodynamic Green's functions in that the derivations given here do not require evaluation of any state vectors or of the partition function. The general results are applied to the Heisenberg ferromagnet problem. At low temperatures, the zeroth-order relations give the Bloch spin-wave results, and an approximate evaluation of the first-order contribution gives the leading term in Dyson's T4{T}^{4}T4 correction to the spontaneous magnetization. For arbitrary temperatures, the zeroth-order relations give the Tyablikov equations for 〈Sz〉〈{S}^{z}〉〈Sz〉, and Callen's results are obtained from the first-order corrections. These examples illustrate the simplicity of calculating statistical averages of functions of the ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, ${\ensuremath{\theta}}_{i}$ operators by the present method.
Thermodynamic functions of an nnn-component magnetic system coupled to an isotropic elastic medium in $d=4\ensuremath{-}\ensuremath{\epsilon}$ dimensions are calculated to leading order in $\ensuremath{\epsilon}$. Solutions to renormalization-group recursion relations correct to $O(\ensuremath{\epsilon})$ along with trajectory-integral methods are employed in these calculations. For n> 2n>~2n> 2 the finite compressibility modifies the critical-point behavior producing new corrections to scaling in the thermodynamic functions which are given. In the n=1n=1n=1 case renormalization-group predictions of a first-order transition are verified. Thermodynamic functions are given and analyzed in this case. In particular the transition temperature and discontinuity in order parameter are determined and the phase diagram is constructed.
The critical behavior of spin systems with quenched disorder is studied by renormalization-group methods. For the randomly dilute mmm-vector model, the n=0n=0n=0 limit is used to construct a translationally invariant effective Hamiltonian which describes the original disordered system. This Hamiltonian is analyzed in the $\ensuremath{\epsilon}$ expansion to order ${\ensuremath{\epsilon}}^{2}$. Sharp second-order phase transitions with exponents which do not depend continuously on impurity concentration are predicted. For $m>{m}_{c}\ensuremath{\equiv}4\ensuremath{-}4\ensuremath{\epsilon}+O({\ensuremath{\epsilon}}^{2})$ the isotropic mmm-component fixed point, which characterizes the critical behavior of the pure system, is stable. For m<mcm<{m}_{c}m<mc​, a new random fixed point becomes stable. The exponents corresponding to this fixed point are $\ensuremath{\eta}=[\frac{(5{m}^{2}\ensuremath{-}8m)}{256{(m\ensuremath{-}1)}^{2}}]{\ensuremath{\epsilon}}^{2}+O({\ensuremath{\epsilon}}^{3})$, $\ensuremath{\nu}=\frac{1}{2}+[\frac{3m}{32(m\ensuremath{-}1)}]\ensuremath{\epsilon}+[\frac{m(127{m}^{2}\ensuremath{-}572m\ensuremath{-}32)}{4096{(m\ensuremath{-}1)}^{3}}]{\ensuremath{\epsilon}}^{2}+O({\ensuremath{\epsilon}}^{3})$ for $m\ensuremath{\ne}1$, and $\ensuremath{\eta}=\ensuremath{-}\frac{\ensuremath{\epsilon}}{106}+O({\ensuremath{\epsilon}}^{\frac{3}{2}})$, $\ensuremath{\nu}=\frac{1}{2}+\frac{{(\frac{6\ensuremath{\epsilon}}{53})}^{\frac{1}{2}}}{4}+O(\ensuremath{\epsilon})$ for m=1m=1m=1. More general random systems are qualitatively discussed from the effective-Hamiltonian viewpoint.