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.     

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.     

Pesin’s Entropy Formula for Systems Between C^1 and C^{1+\alpha }

In this article we give a new observation of Pesin’s entropy formula, motivated from Mañé’s proof of (Ergod Theory Dyn Sys 1:95–102, 1981). Let \(M\) be a compact Riemann manifold and \(f:\,M\rightarrow M\) be a \(C^1\) diffeomorphism on \(M\) . If \(\mu \) is an \(f\) -invariant probability measure which is absolutely continuous relative to Lebesgue measure and nonuniformly-H \(\ddot{\text {o}}\) lder-continuous(see Definition 1.1), then we have Pesin’s entropy formula, i.e., the metric entropy \(h_\mu (f)\) satisfies $$\begin{aligned} h_{\mu }(f)=\int \sum _{\lambda _i(x)> 0}\lambda _i(x)d\mu , \end{aligned}$$ where \(\lambda _1(x)\ge \lambda _2(x)\ge \cdots \ge \lambda _{dim\,M}(x)\) are the Lyapunov exponents at \(x\) with respect to \(\mu .\) Nonuniformly-H \(\ddot{\text {o}}\) lder-continuous is a new notion from probabilistic perspective weaker than \(C^{1+\alpha }.\)      

Comprehensive Bayesian analysis of rare (semi)leptonic and radiative \varvec{B} decays

The available data on \(|\Delta B| = |\Delta S| = 1\) decays are in good agreement with the Standard Model when permitting subleading power corrections of about \(15\,\%\) at large hadronic recoil. Constraining new-physics effects in \(\mathcal {C}_{7}^{\mathrm {}}\) , \(\mathcal {C}_{9}^{\mathrm {}}\) , \(\mathcal {C}_{10}^{\mathrm {}}\) , the data still demand the same size of power corrections as in the Standard Model. In the presence of chirality-flipped operators, all but one of the power corrections reduce substantially. The Bayes factors are in favor of the Standard Model. Using new lattice inputs for \(B\rightarrow K^*\) form factors and under our minimal prior assumption for the power corrections, the favor shifts toward models with chirality-flipped operators. We use the data to further constrain the hadronic form factors in \(B\rightarrow K\) and \(B\rightarrow K^*\) transitions.     

Modeling of HOT (111) HgCdTe MWIR detector for fast response operation

The paper reports on photoelectrical performance of the mid-wave infrared (MWIR) (111) HgCdTe high operating temperature detector for the fast response conditions. Detector structure was simulated with software APSYS by Crosslight Inc. The detailed analysis of the time response as a function of device architecture and applied voltage was performed pointing out optimal working conditions. The time response of the MWIR HgCdTe detector with 50 % cut-off wavelength of \(\lambda _{c} \approx 5.3\, \upmu \hbox {m}\) at \(T = 200\)  K was estimated at the level of \(\tau _{s} \approx \) 2,500 ps for \(V = 100\)  mV and series resistance \(R_{Series} = 510\,\Omega \) . The series resistance’s reduction enables to reach \(\tau _{s}\approx 60\!-\!500\)  ps.     

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)\) .     

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.     

Microcanonical Analysis of the Curie–Weiss Anisotropic Quantum Heisenberg Model in a Magnetic Field

The anisotropic quantum Heisenberg model with Curie-Weiss-type interactions is studied analytically in several variants of the microcanonical ensemble. (Non)equivalence of microcanonical and canonical ensembles is investigated by studying the concavity properties of entropies. The microcanonical entropy \(s(e,\varvec{m})\) is obtained as a function of the energy \(e\) and the magnetization vector \({\varvec{m}}\) in the thermodynamic limit. Since, for this model, \(e\) is uniquely determined by \({\varvec{m}}\) , the same information can be encoded either in \(s(\varvec{m})\) or \(s(e,m_1,m_2)\) . Although these two entropies correspond to the same physical setting of fixed \(e\) and \({\varvec{m}}\) , their concavity properties differ. The entropy \(s_{{\varvec{h}}}(u)\) , describing the model at fixed total energy \(u\) and in a homogeneous external magnetic field \({\varvec{h}}\) of arbitrary direction, is obtained by reduction from the nonconcave entropy \(s(e,m_1,m_2)\) . In doing so, concavity, and therefore equivalence of ensembles, is restored. \(s_{{\varvec{h}}}(u)\) has nonanalyticities on surfaces of co-dimension 1 in the \((u,\varvec{h})\) -space. Projecting these surfaces into lower-dimensional phase diagrams, we observe that the resulting phase transition lines are situated in the positive-temperature region for some parameter values, and in the negative-temperature region for others. In the canonical setting of a system coupled to a heat bath of positive temperatures, the nonanalyticities in the microcanonical negative-temperature region cannot be observed, and this leads to a situation of effective nonequivalence even when formal equivalence holds.     

Line intensity measurements of methane’s \nu _3-band using a cw-OPO

We report on absolute line strength measurements of P(1), R(0) and R(1) singlet lines in the \(3.3\,\upmu\hbox {m}\,\nu _{3}\) (C–H stretching) band of methane \(^{12}\hbox {CH}_4\) at reference temperature \(T=296\)  K. Line strength measurements are performed at low pressure \((P \le 1\hbox { Torr})\) using direct absorption spectroscopy technique based on a widely tunable continuous-wave singly resonant optical parametric oscillator. The \(1\sigma \) overall accuracy in line strength determinations ranges between 7 and 8 % mostly limited by pressure and frequency measurements. A comparison with previous reported values is made. Our results show good agreement with the HITRAN 2012 database.     

QCD analysis of forward neutron production in DIS

Observing light-by-light scattering at the large hadron collider (LHC) has received quite some attention and it is believed to be a clean and sensitive channel to possible new physics. In this paper, we study the diphoton production at the LHC via the process \({{pp}}\rightarrow {{p}}\gamma \gamma {{p}}\rightarrow {{p}}\gamma \gamma {{p}}\) through graviton exchange in the large extra dimension (LED) model. Typically, when we do the background analysis, we also study the double Pomeron exchange of \(\gamma \gamma \) production. We compare its production in the quark–quark collision mode to the gluon–gluon collision mode and find that contributions from the gluon–gluon collision mode are comparable to the quark–quark one. Our result shows, for extra dimension \(\delta =4\) , with an integrated luminosity \(\mathcal{L} = 200\,\mathrm{fb}^{-1}\) at the 14 TeV LHC, that diphoton production through graviton exchange can probe the LED effects up to the scale \({M}_{S}=5.06 (4.51, 5.11)\,\mathrm{TeV}\) for the forward detector acceptance \(\xi _1 (\xi _2, \xi _3)\) , respectively, where \(0.0015<\xi _1<0.5\) , \(0.1<\xi _2<0.5\) , and \(0.0015<\xi _3<0.15\) .     

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.     

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.     

Terahertz generation by nonlinear mixing of laser and its second harmonic in a rippled density plasma

Terahertz radiation generation by second-order nonlinear mixing of laser $ (\omega_{1} ,\,\vec{k}_{1} ) $ and its frequency shifted second harmonic $ \omega_{2} = 2\omega_{1} - \omega ,\,\,\vec{k}_{2} \, $ $ (\omega \ll \omega_{1} ) $ in a plasma, in the presence of an obliquely inclined density ripple of wave number $ \vec{q} $ , are investigated. The lasers exert ponderomotive force on electrons and drive density perturbations at $ (2\omega_{1} ,\,2\vec{k}_{1} - \vec{q}) $ and $ (\omega_{1} - \omega_{2} ,\,\vec{k}_{1} - \vec{k}_{2} - \vec{q}) $ . These perturbations beat with the electron oscillatory velocities due to the lasers to produce a nonlinear current at $ \omega ,\,\vec{k} = 2\vec{k}_{1} - \vec{k}_{2} - \vec{q} $ , resonantly driving the terahertz radiation when $ \vec{q} $ satisfies the phase matching condition. The radiated THz intensity depends on the relative polarization of the lasers and scales as the square of intensity of the fundamental laser and linearly with the square root of the intensity of the second harmonic. The THz emission is maximized when the polarization of the lasers is aligned. These results are consistent with the recent experimental results.     

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.     

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.     

Antiproton low-energy collisions with Ps-atoms and true muonium atoms (μ + μ −)

Three-charge-particle collisions with participation of ultra-slow antiprotons ( \(\overline {\rm {p}}\) ) is the subject of this work. Specifically we compute the total cross sections and corresponding thermal rates of the following three-body reactions: \(\overline {\rm p}+(e^+e^-) \rightarrow \overline {\rm {H}} + e^-\) and \(\overline {\rm p}+(\mu ^+\mu ^-) \rightarrow \overline {\rm {H}}_{\mu } + \mu ^-\) , where \(e^-(\mu ^-)\) is an electron (muon) and \(e^+(\mu ^+)\) is a positron (antimuon) respectively, \(\overline {\rm {H}}=(\overline {\rm p}e^+)\) is an antihydrogen atom and \(\overline {\rm {H}}_{\mu }=(\overline {\rm p}\mu ^+)\) is a muonic antihydrogen atom, i.e. a bound state of \(\overline {\rm {p}}\) and μ ⁺. A set of two-coupled few-body Faddeev-Hahn-type (FH-type) equations is numerically solved in the framework of a modified close-coupling expansion approach.     

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.     

Complex Vector Formalism of Harmonic Oscillator in Geometric Algebra: Particle Mass,Spin and Dynamics in Complex Vector Space

Elementary particles are considered as local oscillators under the influence of zeropoint fields. Such oscillatory behavior of the particles leads to the deviations in their path of motion. The oscillations of the particle in general may be considered as complex rotations in complex vector space. The local particle harmonic oscillator is analyzed in the complex vector formalism considering the algebra of complex vectors. The particle spin is viewed as zeropoint angular momentum represented by a bivector. It has been shown that the particle spin plays an important role in the kinematical intrinsic or local motion of the particle. From the complex vector formalism of harmonic oscillator, for the first time, a relation between mass $m$ and bivector spin $S$ has been derived in the form $\varvec{\sigma }_3 mc^2{\mathcal {J}}_{\pm } =\lambda \Omega _{\mathbf{s}} \cdot \mathrm{{S}} {\mathcal {J}}_{\pm }$ . Where, $\Omega _{s}$ is the angular velocity bivector of complex rotations, $c$ is the velocity of light. The unit vector $\varvec{\sigma }_3$ acts as an operator on the idempotents ${\mathcal {J}}_{+}$ and ${\mathcal {J}}_{-}$ to give the eigen values $\lambda =\pm 1.$ The constant $\lambda $ represents two fold nature of the equation corresponding to particle and antiparticle states. Further the above relation shows that the mass of the particle may be interpreted as a local spatial complex rotation in the rest frame. This gives an insight into the nature of fundamental particles. When a particle is observed from an arbitrary frame of reference, it has been shown that the spatial complex rotation dictates the relativistic particle motion. The mathematical structure of complex vectors in space and spacetime is developed.