Search

A one-dimensional charge density wave (CDW) is driven to slide by a dc electric field, carrying an electric current. In an additional ac field with frequency ${\omega}_{\mathrm{ex}}$, it is known that the sliding CDW can be synchronized to $\omega_{\mathrm{ex}}$, leading to the occurrence of Shapiro steps in the $I$-$V$ characteristics. Motivated by a recent experiment where ac fields with two frequencies $\omega_{\mathrm{ex}}$ and $\omega_{\mathrm{ex}}^{\prime}$ are simultaneously applied, we theoretically investigate the effects of frequency mixing on the Shapiro-step formation. Based on the Fukuyama-Lee-Rice model, we show that in addition to the main steps induced by $\omega_{\mathrm{ex}}$, satellite steps characterized by $\omega_{\mathrm{ex}}^{\prime}$ emerge. It is also found that with increasing the ac-field strength for $\omega_{\mathrm{ex}}^{\prime}$, each step width first exhibits a damped oscillation as in the one-frequency case, and then, exhibits a non-monotonic behavior. The origin of these behaviors and the relevance to the associated experiment are also discussed., Comment: 6 pages, 4 figures + supplemental material

In two-dimensional superconductors with a Rashba-type spin-orbit coupling, it is known that an in-plane magnetic field can induce a helical superconducting (SC) state with a phase modulation $e^{i {\bf q}\cdot {\bf r}}$. Here, we theoretically investigate the stability of a stripe order, a weight-biased superposition state composed of $+{\bf q}$ and $-{\bf q}$ modes taking the form of $\Delta_+ e^{i{\bf q}\cdot{\bf r}}+\Delta_- e^{-i{\bf q}\cdot{\bf r}}$ with $|\Delta_+|\neq|\Delta_-|\neq 0$, assuming that the spin-singlet pairing channel is dominant. Based on the Ginzburg-Landau theory, we show that for both s-wave and d-wave pairing symmetries, the stripe order can appear in the high-field and low-temperature region inside the helical phase and that the transition between the helical and stripe phases is of second order. It is noteworthy that for the d-wave pairing, the stability region of the stripe phase shrinks when the in-plane field is rotated from the nodal direction to the anti-nodal direction. It is also found that the nonreciprocity of the critical current, the so-called SC diode effect, emerges not only in the helical phase but also in the stripe phase, with no clear nonreciprocity anomaly at the helical-stripe transition due to its second-order nature., Comment: 11 pages, 8 figures

Effects of a bias electric current have been theoretically investigated in a spin-triplet superconducting ring in a magnetic field. Based on the Ginzburg-Landau theory, we show that the bias current can stabilize a half-quantum-flux (HQF) state via couplings to the Zeeman field and the dipole-type spin-orbit interaction, the latter becoming active when the field is tilted from the ring axis. The emergence of the HQF state is reflected as a field-induced half-quantum-shift in the Little-Parks (LP) oscillation in the critical current. Possible relevance to recent LP experiments is also discussed., Comment: 5 pages, 2 figures + supplemental material

We theoretically investigate the effects of surface acoustic waves (SAWs) on an electric-field-driven sliding motion of a one-dimensional charge density wave (CDW), which is initially pinned by impurities. By numerically analyzing an extended Fukuyama-Lee-Rice model, we show that a mechanical vibration of the SAW, which, in the model, is assumed to affect the CDW via the pinning site in the form of temporally oscillating pinning parameters, induces Shapiro steps with self-similarity, i.e., the devil's staircase, in the current-voltage characteristics. It is also found that when the SAW acts as the vibration in the pinning strength, the mechanism of the mode locking (harmonic and subharmonic responses) leading to the occurrence of the Shapiro steps is modified, and as a result, the fractal dimension and parameter dependence of the SAW-induced staircase can be considerably different from those for the conventional ac-electric-field-induced one. This suggests that an unconventional type of fractal phenomena can emerge in the SAW-induced CDW dynamics., Comment: 6 pages, 3 figures, + supplemental material

The stability of a miniature skyrmion crystal (SkX) with only a small number of spins in the magnetic unit cell has been theoretically investigated in $J_1$-$J_3$ antiferromagnets on the breathing kagome lattice with a single-ion anisotropy $D$ at zero field. It is found by means of Monte Carlo simulations that due to the breathing bond-alternation, a zero-field triple-${\bf Q}$ miniature SkX can be stabilized not only in the specific case of $D=0$ [K. Aoyama and H. Kawamura, Phys. Rev. B 105, L100407 (2022)] but also in more general situations with easy-axis ($D<0$) and easy-plane ($D>0$) anisotropies which favor triple-${\bf Q}$ collinear and noncoplanar states, respectively. Since the SkX and anti-SkX each having positive or negative chirality are energetically degenerate, the topological Hall effect of alternative sign is possible at zero field. It is also found that reflecting the chiral degeneracy, the collinear and coplanar phases preempting the SkX phase possess random domain structures consisting of positive- and negative-chirality clusters., Comment: 5 pages, 3 figures

We investigate spin and thermal transport near the N\'{e}el transition temperature $T_N$ in three dimensions, by numerically analyzing the classical antiferromagnetic $XXZ$ model on the cubic lattice, where in the model, the anisotropy of the exchange interaction $\Delta=J_z/J_x$ plays a role to control the universality class of the transition. It is found by means of the hybrid Monte-Carlo and spin-dynamics simulations that in the $XY$ and Heisenberg cases of $\Delta \leq 1$, the longitudinal spin conductivity $\sigma^s_{\mu\mu}$ exhibits a divergent enhancement on cooling toward $T_N$, while not in the Ising case of $\Delta>1$. In all the three cases, the temperature dependence of the thermal conductivity $\kappa_{\mu\mu}$ is featureless at $T_N$, being consistent with experimental results. The divergent enhancement of $\sigma^s_{\mu\mu}$ toward $T_N$ is attributed to the spin-current relaxation time which gets longer toward $T_N$, showing a power-law divergence characteristic of critical phenomena. It is also found that in contrast to the $XY$ case where the divergence in $\sigma^s_{\mu\mu}$ is rapidly suppressed below $T_N$, $\sigma^s_{\mu\mu}$ likely remains divergent even below $T_N$ in the Heisenberg case, which might experimentally be observed in the ideally isotropic antiferromagnet RbMnF$_3$., Comment: 17 pages, 9 figures. arXiv admin note: text overlap with arXiv:1908.06630

We theoretically investigate a $J_1$-$J_3$ classical Heisenberg model on the breathing pyrochlore lattice, where the nearest-neighbor (NN) exchange interactions for small and large tetrahedra, $J_1$ and $J_1'$, take different values due to the breathing bond-alternation and $J_3$ is the third NN antiferromagnetic interaction along the bond direction. It is found by means of Monte Carlo simulations that for large $J_3$, a hedgehog lattice, a three-dimensional periodic array of magnetic monopoles and antimonopoles, emerges in the form of a quadruple-${\bf Q}$ state characterized by the ordering vector of ${\bf Q}=(\pm\frac{1}{2},\pm\frac{1}{2},\pm\frac{1}{2})$, being irrespective of the signs of $J_1$ and/or $J_1'$ as long as $J_1\neq J_1'$. It is also found that in an applied magnetic field, there appear six quadruple-${\bf Q}$ states depending on the values of $J_1$ and $J_1'$, among which three phases including the in-field hedgehog-lattice state exhibit nonzero total chirality $\mbox{\boldmath $\chi$}^{\rm T}$ associated with the anomalous Hall effect of chirality origin. In the remaining two chiral phases, which are realized in the presence of ferromagnetic $J_1$ and/or $J_1'$, the spin structure is not topologically nontrivial, in spite of the fact that $\mbox{\boldmath $\chi$}^{\rm T} \neq 0$. The role of the topological objects of the monopoles in $\mbox{\boldmath $\chi$}^{\rm T}$ is also discussed., Comment: 21 pages, 20 figures

We theoretically investigate the critical bias current $j_c$ for a superconducting (SC) ring in a magnetic field. Based on the Ginzburg-Landau theory, we show that $j_c$ exhibits a Little-Parks (LP) oscillation as a function of the magnetic flux passing through the ring, similarly to the LP oscillation in the SC transition temperature. It is also found that for a spin-triplet SC ring, the ${\bf d}$-vector rotates to yield a larger $j_c$, forming a texture along the circumference of the ring. Experimental implications of our result are discussed., Comment: 6 pages, 3 figures

We show that classical Heisenberg antiferromagnets on the breathing kagome lattice can be a platform to realize a zero-field topological order of the scalar spin chirality which can be viewed as a miniature skyrmion crystal (SkX) of discrete form with a small number of spins in its magnetic unit cell. In the model, a third nearest-neighbor (NN) antiferromagnetic interaction along the bond direction $J_3$ and the breathing bond-alternation characterized by the ratio of the NN interaction for large triangles to that for small ones, $J_1'/J_1$, are essential. It is found by means of Monte Carlo simulations that a commensurate triple-${\bf Q}$ state appearing for relatively strong $J_3$ at zero field is the noncoplanar state with the SkX structure in the breathing case of $J_1'/J_1 \neq 1$, while in the uniform case of $J_1'/J_1 =1$, it is a collinear state favored by thermal fluctuations. The origin of this chiral order and experimental implications of our result are also discussed., Comment: 6 pages 3 figures + 5 pages 5 figures

We theoretically investigate spin-lattice coupling (SLC) effects on the in-field ordering properties of classical Heisenberg antiferromagnets on the breathing pyrochlore lattice. Here, we use the two possible simplified models describing the effect of local lattice distortions on the spin ordering via the SLC, the bond-phonon and site-phonon models. It is found by means of Monte Carlo simulations that in both models, the $\frac{1}{2}$ plateau shows up in the magnetization curve being relatively robust against the breathing bond-alternation, although magnetic long-range orders (LRO's) are realized only in the site-phonon model. In the bond-phonon model, additional further neighbor interactions are necessary to induce a magnetic LRO. In the site-phonon model, it is also found that in addition to the low-field, middle-field 1/2-plateau, and high-field phases appearing on both the uniform and breathing pyrochlore lattices, various types of unconventional phases which can be viewed as LRO's in units of tetrahedron are induced by the breathing bond-alternation just below the 1/2 plateau and the saturation field. The occurrence of these tetrahedron-based orders could be attributed to the nature characteristic of the breathing pyrochlore lattice, i.e., the existence of the nonequivalent small and large tetrahedra. Experimental implications of our result are also discussed., Comment: 18 pages, 13 figures

The hedgehog lattice, a three-dimensional periodic array of magnetic monopoles and antimonopoles, is known to be realized in the presence of the Dzyaloshinskii-Moriya (DM) interaction. Here, we demonstrate by means of Monte Carlo simulations that the hedgehog lattice is induced by not the DM interaction but frustration in classical Heisenberg antiferromagnets on the breathing pyrochlore lattice. In the model, the breathing bond-alternation is characterized by the ratio of the nearest-neighbor (NN) antiferromagnetic exchange interaction for large tetrahedra to that for small ones, J1'/J1. A quadruple-q state with the ordering vector of q=(\pm 1/2,\pm 1/2,\pm 1/2), which is realized for a large third NN antiferromagnetic interaction along the bond direction J3, turns out to become the hedgehog-lattice state in the breathing case of J1'/J1 <1, while in the uniform case of J1'/J1 =1, it is a collinear state favored by thermal fluctuations. It is also found that in a magnetic field, the structure of the (1/2,1/2,1/2) hedgehog lattice is changed from cubic to tetragonal, resulting in a nonzero net spin chirality which in a metalic system, should yield a characteristic topological Hall effect., Comment: 16 pages, 12 figures, accepted in PRB

We report a negative resistance, namely, a voltage drop along the opposite direction of a current flow, in the superconducting gap of NbSe$_2$ thin films under the irradiation of surface acoustic waves (SAWs). The amplitude of the negative resistance becomes larger by increasing the SAW power and decreasing temperature. As one possible scenario, we propose that soliton-antisoliton pairs in the charge density wave of NbSe$_2$ modulated by the SAW serve as a time-dependent capacitance in the superconducting state, leading to the dc negative resistance. The present experimental result would provide a previously unexplored way to examine nonequilibrium manipulation of the superconductivity., Comment: 25 pages, 4 figures

Transport properties of the classical antiferromagnetic XXZ model on the square lattice have been theoretically investigated, putting emphasis on how the occurrence of a phase transition is reflected in spin and thermal transports. As is well known, the anisotropy of the exchange interaction $\Delta\equiv J_z/J_x$ plays a role to control the universality class of the transition of the model, i.e., either a second-order transition at $T_N$ into a magnetically ordered state or the Kosterlitz-Thouless (KT) transition at $T_{KT}$, which respectively occur for the Ising-type ($\Delta >1$) and $XY$-type ($\Delta <1$) anisotropies, while for the isotropic Heisenberg case of $\Delta=1$, a phase transition does not occur at any finite temperature. It is found by means of the hybrid Monte-Carlo and spin-dynamics simulations that the spin current probes the difference in the ordering properties, while the thermal current does not. For the $XY$-type anisotropy, the longitudinal spin-current conductivity $\sigma^s_{xx}$ ($=\sigma^s_{yy}$) exhibits a divergence at $T_{KT}$ of the exponential form, $\sigma^s_{xx} \propto \exp\big[ B/\sqrt{T/T_{KT}-1 }\, \big]$ with $B={\cal O}(1)$, while for the Ising-type anisotropy, the temperature dependence of $\sigma^s_{xx}$ is almost monotonic without showing a clear anomaly at $T_{N}$ and such a monotonic behavior is also the case in the Heisenberg-type spin system. The significant enhancement of $\sigma^s_{xx}$ at $T_{KT}$ is found to be due to the exponential rapid growth of the spin-current-relaxation time toward $T_{KT}$, which can be understood as a manifestation of the topological nature of a vortex whose lifetime is expected to get longer toward $T_{KT}$. Possible experimental platforms for the spin-transport phenomena associated with the KT topological transition are discussed., Comment: 20 pages, 15 figures

We theoretically investigate effects of lattice distortions on the spin ordering of classical Heisenberg antiferromagnets on the breathing pyrochlore lattice. In the model, local lattice distortions originating from the site vibration are taken into account to yield effective spin interactions via the spin-lattice coupling (SLC). The breathing lattice-alternation is characterized by the ratio of the nearest-neighbor interaction for large tetrahedra to that for small ones, J'/J. It is found by Monte Carlo simulations that the system exhibits a first-order transition into four different types of collinear magnetic ordered states. In the uniform case (J'/J=1), the state realized at stronger SLC is cubic-symmetric characterized by the (1/2,1/2,1/2) magnetic Bragg peaks, while the one at weaker SLC is tetragonal-symmetric characterized by the (1,1,0) ones. With increasing the breathing distortion (decreasing J'/J), the ordering pattern of the (1/2,1/2,1/2) state becomes non-cubic with its magnetic Bragg reflections almost unchanged, while the (1,1,0) state is robust. The non-cubic state peculiar to the breathing pyrochlores is further categorized into two types. We demonstrate that these two non-cubic orderings result in the massive degeneracy of the ground state and that the associated residual entropy per spin takes characteristic values of k_B ln(4)/16 and k_B ln(6)/16. Experimental implications of the results are discussed., Comment: 14 pages, 9 figures, 1 table

Surface scattering effects in thin films of the superfluid 3He B phase have been theoretically investigated, with an emphasis on the stability of the stripe order with spontaneous broken translational symmetry in the film plane and quasiparticle excitations in this spatially inhomogeneous phase. Based on the Ginzburg-Landau theory in the weak coupling limit, we have shown that the stripe order, which was originally discussed for a film with two specular surfaces, can be stable in a film with one specular and one diffusive surfaces which should correspond to superfluid 3He on a substrate. It is also found by numerically solving the Eilenberger equation that due to the stripe structure, a midgap state distinct from the surface Andreev bound state emerges and its signature is reflected in the local density of states., Comment: 8 pages, 5 figures

Effects of local lattice distortions on the spin ordering are investigated for the antiferromagnetic classical Heisenberg model on the pyrochlore lattice. It is found by Monte Carlo simulations that the spin-lattice coupling (SLC) originating from site phonons induces a first-order transition into two different types of collinear magnetic ordered states. The state realized at stronger SLC is cubic symmetric characterized by the magnetic (1/2,1/2,1/2) Bragg peaks, while that at weaker SLC is tetragonal symmetric characterized by the (1,1,0) ones, each accompanied by the commensurate local lattice distortions. Experimental implications to chromium spinels are discussed., Comment: 5 pages, 4 figures

We theoretically investigate pairing states of the spin-triplet p-wave superfluid 3He confined in narrow cylinders. The surface-induced distortion and the multiple internal degrees of freedom of the order parameter lead to the occurrence of a stripe structure along the cylinder axis in the superfluid 3He-B phase. We show that in sufficiently small cylinders with an anisotropic surface scattering, the stripe order with broken translational symmetry may be stabilized as the lowest energy state. Periodic spin-current textures caused in this inhomogeneous superfluid phase are also discussed., Comment: 5 pages, 3 figures, accepted for publication in PRB (R)

Domains in non-centrosymmetric materials represent regions of different crystal structure and spin-orbit coupling. Twin boundaries separating such domains display unusual properties in non-centrosymmetric superconductors (NCS), where magneto-electric effects influence the local lower and upper critical magnetic fields. As a model system, we investigate NCS with tetragonal crystal structure and Rashba spin-orbit coupling (RSOC), and with twin boundaries parallel to their basal planes. There, we report that there are two types of such twin boundaries which separate domains of opposite RSOC. In a magnetic field parallel to the basal plane, magneto-electric coupling between the spin polarization and supercurrents induces an effective magnetic field at these twin boundaries. We show this leads to unusual effects in such superconductors, and in particular to the modification of the upper and lower critical fields, in ways that depend on the type of twin boundary, as analyzed in detail, both analytically and numerically. Experimental implications of these effects are discussed., Comment: 10 pages, 6 figures

The effect of spatial inhomogeneity on the properties of a two-dimensional non-centrosymmetric superconductor in an in-plane magnetic field is investigated, as it can be realized in LaAlO_3-SrTiO_3 interfaces. We demonstrate that the spatial variation of Rashba spin-orbit coupling (RSOC) yields a local magnetic flux pattern due to the field-induced inhomogeneous helical phase. For sufficiently strong fields, vortices can nucleate at inhomogeneities of the RSOC., Comment: 5 pages, 3 figures, accepted for publication in PRL

Effects of the strong Pauli-paramagnetic pair-breaking (PPB) on the vortex lattice in d-wave superconductors are theoretically studied by putting emphasis on consequences of the PPB-induced antiferromagnetic (AFM) ordering in the spatial modulation in the vortex lattice. It is shown that the PPB-induced AFM fluctuation in the superconducting state leads to an enhancement of the vortex lattice form factor which is a measure of spatial variations of the internal magnetic field and that the enhancement becomes more remarkable as an AFM instability is approached. It is also demonstrated that the PPB-induced AFM ordering is assisted by the vortex-lattice modulation, and thus, that the resulting AFM order is spatially modulated, while it is not localized in the vortex cores but coexistent with the nonvanishing superconducting order parameter. These results are discussed in connection with two phenomena observed in CeCoIn5, the anomalous field dependence of the vortex lattice form factor and the AFM order appearing inside the high-field and low-temperature superconducting phase., Comment: 26 pages, 13 figures

Core structures of a single vortex in A-like and B-like phases of superfluid 3He in uniaxially compressed and stretched aerogels are studied by numerically solving Ginzburg-Landau equations derived microscopically. It is found that, although any uniaxial deformation leads to a wider A-like phase with the axial pairing in the pressure-temperature phase diagram, the vortex core states in the two phases in aerogel depend highly on the type of deformation. In a compressed aerogel, the first-order vortex core transition (VCT) previously seen in the bulk B phase appears at any pressure in the B-like phase while no strange vortex core is expected in the corresponding A-like phase. By contrast, in a stretched aerogel, the VCT in the B-like phase is lost while another VCT is expected to occur between a nonunitary core and a polar one in the A-like phase. Experimental search for these results is hoped to understand correlation between superfluid 3He and aerogel structure., Comment: 7 pages, 6 figures Text was changed. Resubmitted version

Antiferromagnetic (AFM) (or spin-density wave) quantum critical fluctuation enhanced just below H_c2(0) have been often observed in d-wave superconductors with a strong Pauli paramagnetic depairing (PD) including CeCoIn_5. It is shown here that such a tendency of field-induced AFM ordering is a consequence of strong PD and should appear particularly in superconductors with a gap node along the AFM modulation. Two phenomena seen in CeCoIn_5, the anomalous vortex lattice form factor and the AFM order in the Fulde-Ferrell-Larkin-Ovchinnikov state, are explained based on this peculiar PD effect., Comment: 4 pages, 3 figures Title and text were changed. References are added. Resubmitted version

The high field phase diagram and magnetic properties of CeCoIn_5 below H_c2(0) are examined from the picture regarding the high field and low temperature (HFLT) phase as a possible Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) vortex lattice. Crucial roles of antiferromagnetic (AFM) fluctuations enhanced close to H_c2(0) are stressed. The FFLO vortex lattice with a longitudinal modulation parallel to the field is stabilized compared with those with lateral modulations as a consequence of the presence of AFM fluctuations. Further, an unusual field-induced enhancement of the flux distribution is argued to be a consequence of interplay between the paramagnetic depairing and AFM fluctuation, both of which are enhanced with increasing field., Comment: 2 pages, 2 figures, M2S-IX, Tokyo (2009) conference proceedings, Available online in Physica C

Core structures of a single vortex in A-like and B-like phases of superfluid 3He in uniaxially compressed and stretched aerogels are studied by numerically solving Ginzburg-Landau equations derived microscopically. It is found that, although the A-like phase with axial pairing is extended by any uniaxial deformation, the vortex core states in the two phases depend highly on the type of deformation in aerogel: In a compressed aerogel, the first order vortex core transition (VCT) previously seen in the bulk B phase appears at any pressure in B-like phase, while no strange vortex core is expected in A-like phase. By contrast, in a stretched aerogel, the VCT in B-like phase is lost, while another VCT is expected to occur between a nonunitary core and the polar one in A-like phase. An experimental search for these results is hoped to understand correlation between superfluid $^3$He and aerogel structures., Comment: This manuscript is withdrawn because the corresponding extended version is available as arXiv:1007.4045

Recent theoretical and experimental studies of superfluid $^3$He in aerogels with a global anisotropy, e.g., due to an external stress, have definitely shown that the A-like phase with an equal spin pairing (ESP) in such aerogel samples is in the ABM (or, axial) pairing state. In this paper, the A-like phase of superfluid $^3$He in globally {\it isotropic} aerogel is studied in details by assuming a weakly disordered system in which singular topological defects are absent. Through calculation of the free energy, a disordered ABM state is found to be the best candidate of the pairing state of the globally isotropic A-like phase. Further, it is found through a one-loop renormalization group calculation that the coreless continuous vortices (or, vortex-skyrmions) are irrelevant to the long-distance behavior of the disorder-induced textures, and that the superfluidity is maintained in spite of lack of the conventional off-diagonal long range order. Therefore, the globally isotropic A-like phase at weak disorder is, like in the case with a global stretched anisotropy, a superfluid glass with the ABM pairing., Comment: Revised version accepted for publication in Phys.Rev.B

Effects of impurity scatterings on the strong coupling (SC) contribution, stabilizing the ABM (axial) pairing state, to the quartic term of the Ginzburg-Landau (GL) free energy of superfluid $^3$He are theoretically studied to examine recent observations suggestive of an anomalously small SC effect in superfluid $^3$He in aerogels. To study the SC corrections, two approaches are used. One is based on a perturbation in the short-range repulsive interaction, and the other is a phenomenological approach used previously for the bulk liquid by Sauls and Serene [Phys.Rev.B 24, 183 (1981)]. It is found that the impurity scattering favors the BW pairing state and shrinks the region of the ABM pairing state in the T-P phase diagram. In the phenomenological approach, the resulting shrinkage of the ABM region is especially substantial and, if assuming an anisotropy over a large scale in aerogel, leads to justifying the phase diagrams determined experimentally., Comment: 19 pages, 9 figures, Accepted for publication in Phys. Rev. B

Stable pairing states of superfluid 3He in aerogel are examined in the case with a global uniaxial anisotropy which may be created by applying an uniaxial stress to the aerogel. Due to such a global anisotropy, the stability region of an ABM pairing state becomes wider. In an uniaxially stretched aerogel, the pure polar pairing state with a horizontal line node is predicted to occur, as a 3d superfluid phase, over a measurable width just below the superfluid transition at T_c(P). A possible relevance of the present results to the case with no global anisotropy is also discussed., Comment: 4 pages, 4 figures

The equal-spin pairing (ESP) state, the so-called A-like phase, of superfluid $^3$He in aerogels is studied theoretically in the Ginzburg-Landau (GL) region by examining thermodynamics, and the resulting equilibrium phase diagram is mapped out. We find that, among the ABM, planar, and robust pairing states, the ABM state with presumably quasi long-ranged superfluid order is the best candidate of the A-like phase with a strange lowering of the polycritical point (PCP) observed experimentally., Comment: 4 pages, 1 figure, one reference added, accepted for publication in Phys. Rev. B