The critical behavior at the special surface transition and crossover bevavior from special to ordinary surface transition in semi-infinite n-component anisotropic cubic models are investigated by applying the field theoretic approach directly in d=3 dimensions up to the two-loop approximation. The crossover behavior for random semi-infinite Ising-like system, which is the nontrivial particular case of the cubic model in the limit
n→0, is also investigated. The numerical estimates of the resulting two-loop series expansions for the critical exponents of the special surface transition, surface crossover critical exponent
Φ and the surface critical exponents of the layer,
α1, and local specific heats,
α11, are computed by means of Pade and Pade-Borel resummation techniques. For
nnc the cubic fixed point is stable, where
nc is the marginal spin dimensionality of the cubic model. The obtained results indicate that the surface critical behavior of semi-infinite systems with cubic anisotropy is characterized by new set of surface critical exponents for
n>nc.
The crossover behavior of the susceptibilities ${\ensuremath{\chi}}^{\ensuremath{\alpha}\ensuremath{\alpha}}$ of spin systems with anisotropic exchange coupling is discussed on the basis of an "extended" scaling hypothesis which gives ${\ensuremath{\chi}}^{\ensuremath{\alpha}\ensuremath{\alpha}}(T,g)\ensuremath{\approx}A{t}^{\ensuremath{-}\ensuremath{\gamma}}{X}_{\ensuremath{\alpha}}(\frac{Bg}{{t}^{\ensuremath{\varphi}}})$, where
g is the anisotropy parameter, $t=\frac{(T\ensuremath{-}{T}_{c0})}{{T}_{c0}}$, with
Tc0 and $\ensuremath{\gamma}$ being the iso tropic critical temperature, and susceptibility exponent, while
A and
B are model-dependent amplitudes. Analysis of high-temperature series expansions for ${\ensuremath{\chi}}^{\ensuremath{\alpha}\ensuremath{\alpha}}(T,g)$ as polynomials in
g for rhombic and axial anisotropy in the fcc, bcc, and simple cubic classical Heisenberg (
n=3) and classical
XY and planar-spin models (
N=2), verifies the scaling with crossover exponents $\ensuremath{\varphi}=1.25\ifmmode\pm\else\textpm\fi{}0.015$ (
n=3) and $\ensuremath{\varphi}=1.175\ifmmode\pm\else\textpm\fi{}0.015$ (
n=2). The universality of the scaling functions ${X}_{\ensuremath{\alpha}}(x)$ is demonstrated both for small
x and in the anisotropic limit $T\ensuremath{\rightarrow}{T}_{c}(g)$, where $x\ensuremath{\rightarrow}\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{x}$; for
n=3, accurate representations are constructed in the form ${X}_{\ensuremath{\alpha}}(x)\ensuremath{\simeq}{P}_{\ensuremath{\alpha}}(\frac{x}{\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{x}}){(1\ensuremath{-}\frac{x}{\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{x}})}^{\ensuremath{-}\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{\ensuremath{\gamma}}}$, where $\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{\ensuremath{\gamma}}$ is the anisotropic susceptibility exponent, while ${P}_{\ensuremath{\alpha}}(z)$ varies smoothly from ${P}_{\ensuremath{\alpha}}(0)=1$ to ${P}_{\ensuremath{\alpha}}(1)\ensuremath{\simeq}1.10$.