Antineutrino science in KamLAND

The primary goal of KamLAND is a search for the oscillation of \({\bar{\nu }}_\mathrm{e}\) ’s emitted from distant power reactors. The long baseline, typically 180 km, enables KamLAND to address the oscillation solution of the “solar neutrino problem” with \({\bar{\nu }}_{e} \) ’s under laboratory conditions. KamLAND found fewer reactor \({\bar{\nu }}_{e} \) events than expected from standard assumptions about \(\overline{\nu }_e\) propagation at more than 9 \(\sigma \) confidence level (C.L.). The observed energy spectrum disagrees with the expected spectral shape at more than 5 \(\sigma \) C.L., and prefers the distortion from neutrino oscillation effects. A three-flavor oscillation analysis of the data from KamLAND and KamLAND + solar neutrino experiments with CPT invariance, yields \(\Delta m_{21}^2 \) = [ \(7.54_{-0.18}^{+0.19} \times \) 10 \(^{-5}\) eV \(^{2}\) , \(7.53_{-0.18}^{+0.19} \times \) 10 \(^{-5}\) eV \(^{2}\) ], tan \(^{2}\theta _{12}\) = [ \(0.481_{-0.080}^{+0.092} \) , \(0.437_{-0.026}^{+0.029} \) ], and sin \(^{2}\theta _{13}\) = [ \(0.010_{-0.034}^{+0.033} \) , \(0.023_{-0.015}^{+0.015} \) ]. All solutions to the solar neutrino problem except for the large mixing angle region are excluded. KamLAND also demonstrated almost two cycles of the periodic feature expected from neutrino oscillation effects. KamLAND performed the first experimental study of antineutrinos from the Earth’s interior so-called geoneutrinos (geo \({\bar{\nu }}_{e} \) ’s), and succeeded in detecting geo \({\bar{\nu }}_{e} \) ’s produced by the decays of \(^{238}\) U and \(^{232}\) Th within the Earth. Assuming a chondritic Th/U mass ratio, we obtain \(116_{-27}^{+28} {\bar{\nu }}_{e}\) events from \(^{238}\) U and \(^{232}\) Th, corresponding a geo \({\bar{\nu }}_{e}\) flux of \(3.4_{-0.8}^{+0.8}\times \) 10 \(^{6}\) cm \(^{-2}\)  s \(^{-1}\) at the KamLAND location. We evaluate various bulk silicate Earth composition models using the observed geo \({\bar{\nu }}_{e} \) rate.     

Investigation of conductivity of adhesive layer including indium tin oxide particles for multi-junction solar cells

We report connection conductivity ( \(C_{\rm c}\) ) of adhesive which including \(\hbox {In}_2\hbox {O}_3\) – \(\hbox {SnO}_2\) (ITO) particles developed for fabrication of stacked-type-multi-junction solar cells. The commercial 20- \(\upmu \) m sized ITO particles were heated in vacuum at temperature ranging from 800 to 1,300  \(^{\circ }{\rm C}\) for 10 min to increase \(C_{\rm c}\) . 6.2 wt% ITO particles were dispersed in commercial Cemedine adhesive gel to form 100 samples structured with n-type Si/adhesive/n-type Si (n-Si sample) and p-type Si/adhesive/p-type Si (p-Si sample). Current density as a function of voltage (J–V) characteristics gave \(C_{\rm c}\) . It ranged from 4.3 to 1.0 S/cm \(^2\) for the n-Si sample with 800 \(^{\circ }{\rm C}\) heat-treated ITO particles. Its standard deviation was 0.59 S/cm \(^2\) . On the other hand, it ranged from 2.0 to 0.6 S/cm \(^2\) for the p-Si sample with 800  \(^{\circ }{\rm C}\) heat-treated ITO particles. Its standard deviation was 0.22 S/cm \(^2\) . The distribution of \(C_{\rm c}\) mainly resulted from contact efficiency of ITO particles to substrate. We theoretically estimated that present \(C_{\rm c}\) achieved a low loss of the power conversion efficiency ( \(E_{\rm ff}\) ) lower than 0.3 % in the application of fabrication of multi-junction solar cell with an intrinsic \(E_{\rm ff}\) of 30 % and an open circuit voltage above 1.9 V.     

Optimal N-Point Configurations on the Sphere: “Magic” Numbers and Smale’s 7th Problem

This paper inquires into the concavity of the map \(N\mapsto v_s(N)\) from the integers \(N\ge 2\) into the minimal average standardized Riesz pair-energies \(v_s(N)\) of \(N\) -point configurations on the sphere \(\mathbb {S}^2\) for various \(s\in \mathbb {R}\) . The standardized Riesz pair-energy of a pair of points on \(\mathbb {S}^2\) a chordal distance \(r\) apart is \(V_s(r)= s^{-1}\left( r^{-s}-1 \right) \) , \(s \ne 0\) , which becomes \(V_0(r) = \ln \frac{1}{r}\) in the limit \(s\rightarrow 0\) . Averaging it over the \(\left( \begin{array}{c} N\\ 2\end{array}\right) \) distinct pairs in a configuration and minimizing over all possible \(N\) -point configurations defines \(v_s(N)\) . It is known that \(N\mapsto v_s(N)\) is strictly increasing for each \(s\in \mathbb {R}\) , and for \(s<2\) also bounded above, thus “overall concave.” It is (easily) proved that \(N\mapsto v_{-2}^{}(N)\) is even locally strictly concave, and that so is the map \(2n\mapsto v_s(2n)\) for \(s<-2\) . By analyzing computer-experimental data of putatively minimal average Riesz pair-energies \(v_s^x(N)\) for \(s\in \{-1,0,1,2,3\}\) and \(N\in \{2,\ldots ,200\}\) , it is found that the map \(N\mapsto {v}_{-1}^x(N)\) is locally strictly concave, while \(N\mapsto {v}_s^x(N)\) is not always locally strictly concave for \(s\in \{0,1,2,3\}\) : concavity defects occur whenever \(N\in {\mathcal {C}}^{x}_+(s)\) (an \(s\) -specific empirical set of integers). It is found that the empirical map \(s\mapsto {\mathcal {C}}^{x}_+(s),\ s\in \{-2,-1,0,1,2,3\}\) , is set-theoretically increasing; moreover, the percentage of odd numbers in \({\mathcal {C}}^{x}_+(s),\ s\in \{0,1,2,3\}\) is found to increase with \(s\) . The integers in \({\mathcal {C}}^{x}_+(0)\) are few and far between, forming a curious sequence of numbers, reminiscent of the “magic numbers” in nuclear physics. It is conjectured that these new “magic numbers” are associated with optimally symmetric optimal-log-energy \(N\) -point configurations on \(\mathbb {S}^2\) . A list of interesting open problems is extracted from the empirical findings, and some rigorous first steps toward their solutions are presented. It is emphasized how concavity can assist in the solution to Smale’s \(7\) th Problem, which asks for an efficient algorithm to find near-optimal \(N\) -point configurations on \(\mathbb {S}^2\) and higher-dimensional spheres.     

Split supersymmetry under GUT and current dark matter constraints

There are four types of two-Higgs doublet models under a discrete \(Z_2\) symmetry imposed to avoid tree-level flavor-changing neutral current, i.e. type-I, type-II, type-X, and type-Y models. We investigate the possibility to discriminate the four models in the light of the flavor physics data, including \(B_s\) – \(\bar{B}_s\) mixing, \(B_{s,d} \rightarrow \mu ^+ \mu ^-\) , \(B\rightarrow \tau \nu \) and \(\bar{B} \rightarrow X_s \gamma \) decays, the recent LHC Higgs data, the direct search for charged Higgs at LEP, and the constraints from perturbative unitarity and vacuum stability. After deriving the combined constraints on the Yukawa interaction parameters, we have shown that the correlation between the mass eigenstate rate asymmetry \(A_{\Delta \Gamma }\) of \(B_{s} \rightarrow \mu ^+ \mu ^-\) and the ratio \(R=\mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{exp}/ \mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{SM}\) could be a sensitive probe to discriminate the four models with future precise measurements of the observables in the \(B_{s} \rightarrow \mu ^+ \mu ^-\) decay at LHCb.     

NMR in Magnetic Field of the Earth: Pre-Polarization of Nuclei with Alternating Magnetic Field

The polarization of nuclei in the low static magnetic field \(B_0\) with an alternating magnetic field \(B^{*} (B^{*} \gg B_0)\) at a very low frequency \(f_m\) (but \(f_m\gg 1\) / \({T_1}\) , where \(T_1\) is the spin-lattice relaxation time) has been investigated. The question of the optimization of the energy consumption during the pre-polarization is also considered. The possibilities of the method are illustrated by the observation of nuclear magnetic resonance signals from a few liquids.     

Semiclassical strings in supergravity PFT

In this paper, we introduce the bulk viscosity in the formalism of modified gravity theory in which the gravitational action contains a general function \(f(R,T)\) , where \(R\) and \(T\) denote the curvature scalar and the trace of the energy–momentum tensor, respectively, within the framework of a flat Friedmann–Robertson–Walker model. As an equation of state for a prefect fluid, we take \(p=(\gamma -1)\rho \) , where \(0 \le \gamma \le 2\) and a viscous term as a bulk viscosity due to the isotropic model, of the form \(\zeta =\zeta _{0}+\zeta _{1}H\) , where \(\zeta _{0}\) and \(\zeta _{1}\) are constants, and \(H\) is the Hubble parameter. The exact non-singular solutions to the corresponding field equations are obtained with non-viscous and viscous fluids, respectively, by assuming a simplest particular model of the form of \(f(R,T) = R+2f(T)\) , where \(f(T)=\alpha T\) ( \(\alpha \) is a constant). A big-rip singularity is also observed for \(\gamma <0\) at a finite value of cosmic time under certain constraints. We study all possible scenarios with the possible positive and negative ranges of \(\alpha \) to analyze the expansion history of the universe. It is observed that the universe accelerates or exhibits a transition from a decelerated phase to an accelerated phase under certain constraints of \(\zeta _0\) and \(\zeta _1\) . We compare the viscous models with the non-viscous one through the graph plotted between the scale factor and cosmic time and find that the bulk viscosity plays a major role in the expansion of the universe. A similar graph is plotted for the deceleration parameter with non-viscous and viscous fluids and we find a transition from decelerated to accelerated phase with some form of bulk viscosity.     

The Colored Hofstadter Butterfly for the Honeycomb Lattice

We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number \(\sigma _\mathrm{H}\) is given as the winding number of an eigenvector of a \(2 \times 2\) transfer matrix, as a function of the quasi-momentum \(k\in (0,2\pi )\) . This method is computationally efficient (of order \(\mathcal {O}(n^4)\) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for \(\sigma _\mathrm{H}\) for flux \(p/q\) in the \(r\) -th gap conforms with the Diophantine equation \(r=\sigma _\mathrm{H}\cdot p+ s\cdot q\) , which determines \(\sigma _\mathrm{H}\mod q\) . A window such as \(\sigma _\mathrm{H}\in (-q/2,q/2)\) , or possibly shifted, provides a natural further condition for \(\sigma _\mathrm{H}\) , which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition \(\sigma _\mathrm{H}\in (-q,q)\) .     

Reinforced Brownian Motion: A Prototype

We analyze the Brownian Motion limit of a prototypical unit step reinforced random-walk on the half-line. A reinfoced random walk is one which changes the weight of any edge (or vertex) visited to increase the frequency of return visits. The generating function for the discrete case is first derived for the joint probability distribution of \(S_N\) (the location of the walker at the \(N^{th}\) step) and \(A_N\) , the maximum location the walker achieved in \(N\) steps. Then the bulk of the analysis concerns the statistics of the limiting Brownian walker, and of its “environment”, both parametrized by the amplitude \(\delta \) of the reinforcement. The walker marginal distribution can be interpreted as that of free diffusion with a source serving as a diffusing soft confinement, details depending very much on the value of \(-1< \delta < \infty \) .     

Tracy–Widom at High Temperature

We investigate the marginal distribution of the bottom eigenvalues of the stochastic Airy operator when the inverse temperature \(\beta \) tends to \(0\) . We prove that the minimal eigenvalue, whose fluctuations are governed by the Tracy–Widom \(\beta \) law, converges weakly, when properly centered and scaled, to the Gumbel distribution. More generally we obtain the convergence in law of the marginal distribution of any eigenvalue with given index \(k\) . Those convergences are obtained after a careful analysis of the explosion times process of the Riccati diffusion associated to the stochastic Airy operator. We show that the empirical measure of the explosion times converges weakly to a Poisson point process using estimates proved in Dumaz and Virág (Ann Inst H Poincaré Probab Statist 49(4):915–933, 2013). We further compute the empirical eigenvalue density of the stochastic Airy ensemble on the macroscopic scale when \(\beta \rightarrow 0\) . As an application, we investigate the maximal eigenvalues statistics of \(\beta _N\) -ensembles when the repulsion parameter \(\beta _N\rightarrow 0\) when \(N\rightarrow +\infty \) . We study the double scaling limit \(N\rightarrow +\infty , \beta _N \rightarrow 0\) and argue with heuristic and numerical arguments that the statistics of the marginal distributions can be deduced following the ideas of Edelman and Sutton (J Stat Phys 127(6):1121–1165, 2007) and Ramírez et al. (J Am Math Soc 24:919–944, 2011) from our later study of the stochastic Airy operator.     

Phase Transitions in Layered Systems

We consider the Ising model on \(\mathbb Z\times \mathbb Z\) where on each horizontal line \(\{(x,i), x\in \mathbb Z\}\) , called “layer”, the interaction is given by a ferromagnetic Kac potential with coupling strength \(J_{ \gamma }(x,y)={ \gamma }J({ \gamma }(x-y))\) , where \(J(\cdot )\) is smooth and has compact support; we then add a nearest neighbor ferromagnetic vertical interaction of strength \({ \gamma }^{A}\) , where \(A\ge 2\) is fixed, and prove that for any \(\beta \) larger than the mean field critical value there is a phase transition for all \({ \gamma }\) small enough.     

Efficient Algorithms for the Two-Dimensional Ising Model with a Surface Field

The bond propagation and site propagation algorithms are extended to the two-dimensional (2D) Ising model with a surface field. With these algorithms we can calculate the free energy, internal energy, specific heat, magnetization, correlation functions, surface magnetization, surface susceptibility and surface correlations. The method can handle continuous and discrete bond and surface-field disorder and is especially efficient in the case of bond or site dilution. To test these algorithms, we study the wetting transition of the 2D Ising model, which was solved exactly by Abraham. We can locate the transition point accurately with a relative error of \(10^{-8}\) . We carry out the calculation of the specific heat and surface susceptibility on lattices with sizes up to \(200^2 \times 200\) . The results show that a finite jump develops in the specific heat and surface susceptibility at the transition point as the lattice size increases. For lattice size \(320^2 \times 320\) the parallel correlation length exponent is \(1.86\) , while Abraham’s exact result is \(2.0\) . The perpendicular correlation length exponent for lattice size \(160^2\times 160\) is \(1.05\) , whereas its exact value is \(1.0\) .     

Maximal Distance Travelled by N Vicious Walkers Till Their Survival

We consider N Brownian particles moving on a line starting from initial positions \(\mathbf{{u}}\equiv \{u_1,u_2,\ldots u_N\}\) such that \(0 . Their motion gets stopped at time \(t_s\) when either two of them collide or when the particle closest to the origin hits the origin for the first time. For \(N=2\) , we study the probability distribution function \(p_1(m|\mathbf{{u}})\) and \(p_2(m|\mathbf{{u}})\) of the maximal distance travelled by the \(1^{\text {st}}\) and \(2^{\text {nd}}\) walker till \(t_s\) . For general N particles with identical diffusion constants \(D\) , we show that the probability distribution \(p_N(m|\mathbf{u})\) of the global maximum \(m_N\) , has a power law tail \(p_i(m|\mathbf{{u}}) \sim {N^2B_N\mathcal {F}_{N}(\mathbf{u})}/{m^{\nu _N}}\) with exponent \(\nu _N =N^2+1\) . We obtain explicit expressions of the function \(\mathcal {F}_{N}(\mathbf{u})\) and of the N dependent amplitude \(B_N\) which we also analyze for large N using techniques from random matrix theory. We verify our analytical results through direct numerical simulations.     

A flexible scintillation light apparatus for rare event searches

Compelling experimental evidences of neutrino oscillations and their implication that neutrinos are massive particles have given neutrinoless double beta decay ( \(\beta \beta 0\nu \) ) a central role in astroparticle physics. In fact, the discovery of this elusive decay would be a major breakthrough, unveiling that neutrino and antineutrino are the same particle and that the lepton number is not conserved. It would also impact our efforts to establish the absolute neutrino mass scale and, ultimately, understand elementary particle interaction unification. All current experimental programs to search for \(\beta \beta 0\nu \) are facing with the technical and financial challenge of increasing the experimental mass while maintaining incredibly low levels of spurious background. The new concept described in this paper could be the answer which combines all the features of an ideal experiment: energy resolution, low cost mass scalability, isotope choice flexibility and many powerful handles to make the background negligible. The proposed technology is based on the use of arrays of silicon detectors cooled to 120 K to optimize the collection of the scintillation light emitted by ultra-pure crystals. It is shown that with a 54 kg array of natural CaMoO \(_4\) scintillation detectors of this type it is possible to yield a competitive sensitivity on the half-life of the \(\beta \beta 0\nu \) of \(^{100}\) Mo as high as \(\sim \) \(10^{24}\)  years in only 1 year of data taking. The same array made of \(^{40}\) Ca \(^{\mathrm {nat}}\) MoO \(_4\) scintillation detectors (to get rid of the continuous background coming from the two neutrino double beta decay of \(^{48}\) Ca) will instead be capable of achieving the remarkable sensitivity of \(\sim \) \(10^{25}\)  years on the half-life of \(^{100}\) Mo \(\beta \beta 0\nu \) in only 1 year of measurement.     

Toward a Mathematical Holographic Principle

In work started in [17] and continued in this paper our objective is to study selectors of multivalued functions which have interesting dynamical properties, such as possessing absolutely continuous invariant measures. We specify the graph of a multivalued function by means of lower and upper boundary maps \(\tau _{1}\) and \(\tau _{2}.\) On these boundary maps we define a position dependent random map \(R_{p}=\{\tau _{1},\tau _{2};p,1-p\},\) which, at each time step, moves the point \(x\) to \(\tau _{1}(x)\) with probability \(p(x)\) and to \(\tau _{2}(x)\) with probability \(1-p(x)\) . Under general conditions, for each choice of \(p\) , \(R_{p}\) possesses an absolutely continuous invariant measure with invariant density \(f_{p}.\) Let \(\varvec{\tau }\) be a selector which has invariant density function \(f.\) One of our objectives is to study conditions under which \(p(x)\) exists such that \(R_{p}\) has \(f\) as its invariant density function. When this is the case, the long term statistical dynamical behavior of a selector can be represented by the long term statistical behavior of a random map on the boundaries of \(G.\) We refer to such a result as a mathematical holographic principle. We present examples and study the relationship between the invariant densities attainable by classes of selectors and the random maps based on the boundaries and show that, under certain conditions, the extreme points of the invariant densities for selectors are achieved by bang-bang random maps, that is, random maps for which \(p(x)\in \{0,1\}.\)      

Scaling Limits and Critical Behaviour of the 4-Dimensional n-Component |\varphi |^4 Spin Model

We consider the \(n\) -component \(|\varphi |^4\) spin model on \({\mathbb {Z}}^4\) , for all \(n \ge 1\) , with small coupling constant. We prove that the susceptibility has a logarithmic correction to mean field scaling, with exponent \(\frac{n+2}{n+8}\) for the logarithm. We also analyse the asymptotic behaviour of the pressure as the critical point is approached, and prove that the specific heat has fractional logarithmic scaling for \(n =1,2,3\) ; double logarithmic scaling for \(n=4\) ; and is bounded when \(n>4\) . In addition, for the model defined on the \(4\) -dimensional discrete torus, we prove that the scaling limit as the critical point is approached is a multiple of a Gaussian free field on the continuum torus, whereas, in the subcritical regime, the scaling limit is Gaussian white noise with intensity given by the susceptibility. The proofs are based on a rigorous renormalisation group method in the spirit of Wilson, developed in a companion series of papers to study the 4-dimensional weakly self-avoiding walk, and adapted here to the \(|\varphi |^4\) model.     

Near-Extreme Eigenvalues and the First Gap of Hermitian Random Matrices

We study the phenomenon of “crowding” near the largest eigenvalue \(\lambda _\mathrm{max}\) of random \(N \times N\) matrices belonging to the Gaussian Unitary Ensemble of random matrix theory. We focus on two distinct quantities: (i) the density of states (DOS) near \(\lambda _\mathrm{max}\) , \(\rho _\mathrm{DOS}(r,N)\) , which is the average density of eigenvalues located at a distance \(r\) from \(\lambda _\mathrm{max}\) and (ii) the probability density function of the gap between the first two largest eigenvalues, \(p_\mathrm{GAP}(r,N)\) . In the edge scaling limit where \(r = \mathcal{O}(N^{-1/6})\) , which is described by a double scaling limit of a system of unconventional orthogonal polynomials, we show that \(\rho _\mathrm{DOS}(r,N)\) and \(p_\mathrm{GAP}(r,N)\) are characterized by scaling functions which can be expressed in terms of the solution of a Lax pair associated to the Painlevé XXXIV equation. This provides an alternative and simpler expression for the gap distribution, which was recently studied by Witte et al. in Nonlinearity 26:1799, 2013. Our expressions allow to obtain precise asymptotic behaviors of these scaling functions both for small and large arguments.     

Effect of vacancies on magnetic behaviors of Cu-doped 6H-SiC

Magnetism in Cu-doped, Cu \(\rm _{Si}\) –V \(\rm _{Si}\) codoped, or Cu \(\rm _{Si}\) –V \(\rm _{C}\) codoped 6H-SiC are investigated using the first principle. The total density of states for the ferromagnetic Cu \(\rm _{Si}\) at doping concentration of 0.926 at. \(\%\) shows half-metallic behavior, which leads to the total magnetic moment of 2.84  \(\rm \mu _{B}\) per supercell. The total magnetic moment increases with increasing Cu content. The long-range ferromagnetic interaction between Cu atoms can be attributed to the C-mediated double exchange through the strong \(3d\) ? \(2p\) interaction between Cu and neighboring C ones. It is important to note that both V \(\rm _{Si}\) and V \(\rm _{C}\) play a negative role in ferromagnetic coupling between Cu ions. So, to obtain a larger magnetic moment from Cu-doped 6H–SiC, we should try to avoid the appearance of V \(\rm _{Si}\) and V \(\rm _{C}\) during the process of sample preparation. Our theoretical calculations give a valuable insight on how to get a large magnetic moment from Cu-doped 6H–SiC.     

Is the CNMSSM more credible than the CMSSM?

“Post-sphaleron baryogenesis”, a fresh and profound mechanism of baryogenesis accounts for the matter–antimatter asymmetry of our present universe in a framework of Pati–Salam symmetry. We attempt here to embed this mechanism in a non-SUSY SO(10) grand unified theory by reviving a novel symmetry breaking chain with Pati–Salam symmetry as an intermediate symmetry breaking step and as well to address post-sphaleron baryogenesis and neutron–antineutron oscillation in a rational manner. The Pati–Salam symmetry based on the gauge group \(\mathrm{SU}(2)_L \times \mathrm{SU}(2)_{R} \times \mathrm{SU}(4)_C\) is realized in our model at \(10^{5}\) – \(10^{6}\)  GeV and the mixing time for the neutron–antineutron oscillation process having \(\Delta B=2\) is found to be \(\tau _{n-\bar{n}} \simeq 10^{8}\) – \(10^{10}\)  s with the model parameters, which is within the reach of forthcoming experiments. Other novel features of the model include low scale right-handed \(W^{\pm }_R\) , \(Z_R\) gauge bosons, explanation for neutrino oscillation data via the gauged inverse (or extended) seesaw mechanism and most importantly TeV scale color sextet scalar particles responsible for an observable \(n\) – \(\bar{n}\) oscillation which may be accessible to LHC. We also look after gauge coupling unification and an estimation of the proton lifetime with and without the addition of color sextet scalars.     

Electroweak radiative corrections to $$W^+W^-\gamma $$ production at the ILC

We consider holographic superconductors in a rotating black string spacetime. In view of the mandatory introduction of the \(A_\varphi \) component of the vector potential we are left with three equations to be solved. Their solutions show that the rotation parameter \(a\) influences the critical temperature \(T_\mathrm{c}\) and the conductivity \(\sigma \) in a simple but non-trivial way.     

Ga-vacancy-induced room-temperature ferromagnetic and adjusted-band-gap behaviors in GaN nanoparticles

Room-temperature ferromagnetism has been found in Ga-deficient GaN grown using the direct reaction of Ga \(_{2}\) O \(_{3}\) powder with NH \(_{3}\) gas. The observed magnetism in GaN induced by Ga vacancies is investigated both experimentally and theoretically. First-principles calculations reveal that the spontaneous spin polarization is created by the 3.0  \(\mu _\mathrm{B}\) local moment for GaN and magnetism originates from the polarization of the unpaired 2 \(p\) electrons of N surrounding the Ga vacancy. At the same time, the band gap can be also adjusted by changing the Ga-vacancy concentration.