Splitting Parameters Research Articles

Investigation of PCVSi− defect in 4H–SiC as a candidate for a qubit

Abstract Exploration of spin defects in semiconductors for possible qubits encourages the development of the quantum field. Silicon carbide (SiC) is a suitable platform to carry spin defects, due to its excellent electrical, mechanical and optical properties, together with its convenience for crystallographic growth and doping processes. In this study, a negatively charged phosphorus-vacancy (PCVSi −) defect, consisting of a silicon vacancy and nearby substitution of a phosphorus atom to a carbon atom in 4H–SiC, is investigated by first-principles calculations. This defect is demonstrated to possess a high spin (S = 1) with relatively low formation energy. Computed zero-phonon line energy and zero-field splitting parameters of this defect are close to those of neutral divacancy (VCVSi 0), negatively charged nitrogen-vacancy center (NCVSi −) and some other color centers, which indicate a similarity of both optical and spin properties among them. Moreover, the electron spin coherence time of this defect turns out to be 1.15–1.40 ms. Such a long coherence time provides the defect with reliability for quantum information processing. Our results show that the PCVSi − defect can be a promising candidate for a qubit.

Transition-Metal-Related Quantum Emitters in Wurtzite AlN and GaN.

Transition-metal centers exhibit a paramagnetic ground state in wide-bandgap semiconductors and are promising for nanophotonics and quantum information processing. Specifically, there is a growing interest in discovering prominent paramagnetic spin defects that can be manipulated using optical methods. Here, we investigate the electronic structure and magneto-optical properties of Cr and Mn substitutional centers in wurtzite AlN and GaN. We use state-of-the-art hybrid density functional theory calculations to determine level structure, stability, optical signatures, and magnetic properties of these centers. The excitation energies are calculated using the constrained occupation approach and rigorously verified with the complete active space configuration interaction approach. Our simulations of the photoluminescence spectra indicate that in AlN and in GaN are responsible for the observed narrow quantum emission near 1.2 eV. We compute the zero-field splitting (ZFS) parameters and outline an optical spin polarization protocol for and . Our results demonstrate that these centers are promising candidates for spin qubits.

Local Structure and Optical Studies of Mn2+ Doped L-histidine-4- nitrophenolate 4-nitrophenol Single Crystal

Background: The zero-field splitting parameters of Mn2+ doped L-histidine-4- nitrophenolate 4-nitrophenol single crystals are evaluated Methods: The superposition model and perturbation theory are used to obtain zero-field splitting parameters for Mn2+ ion-doped LHPP single crystals. The optical spectra of the system are computed using the crystal field parameters from the superposition model as input into the crystal field analysis program. Results: The evaluated zero field splitting parameters are in good match with the experimental values when local distortion is taken into account. The experimental finding that the Mn2+ ion enters the L-histidine-4-nitrophenolate 4-nitrophenol lattice at the interstitial position is supported by the theoretical result. Conclusion: It is found that the calculated and experimental band positions agree fairly well. Thus, the theoretical study supports the experimental observation.

High resolution analysis of the CD[formula omitted] deuterated methane: Extended investigation of the pentad region

A highly accurate rotational–vibrational analysis of Fourier transform infrared spectra of the 12CD4 molecule is presented. The high resolution infrared spectra were measured with a IFS125 HR Fourier transform interferometer from Bruker at an optical resolution of 0.003 cm−1 and analyzed in the 1750–2400 cm−1 region. Here the 2ν2, ν2+ν4, 2ν4, ν1 and ν3 bands (altogether, nine sub-bands of different symmetry) of the pentad are located. The number of 1213/1993/1576/77/1582 transitions with the Jmax = 23/23/23/14/32 were assigned to the 2ν2, ν2+ν4, 2ν4, ν1 and ν3 bands of 12CD4. The obtained experimental data were used for the determination of the upper ro-vibrational energy values. To provide more correct values of the upper energies, more than 7800 highly accurate “hot” transitions from the dyad region were additionally processed. In general, 4088 upper ro-vibrational energies of the pentad (for comparison, 2525 upper ro-vibrational energies with the value of Jmax=20 are known in the modern literature up to now) were determined, which were used then in the weighted fit procedure with a goal to determine the spectroscopic parameters (band centers, rotational, centrifugal distortion, tetrahedral splitting and resonance interaction parameters) of the effective Hamiltonian. The obtained drms value is 5.5×10−4 cm−1 which is almost one hundred times better than the reproduction of the same set of experimental data by the parameters known in the earlier literature.

Unraveling the Stability and Magnetic Properties of Bis-Hydrated Mn(II) Complexes via Tailored Ligand Design.

Exploring the electronic structure and dynamic behavior of Mn(II) complexes reveals fascinating magnetic properties and prospective biomedical applications. In this study, we investigate the solvent phase dynamics of heptacoordinated Mn(II) complexes through ab initio molecular dynamics simulations and density functional theory (DFT) calculations with effectively varying temperatures. We observed that the complex with high stability ([Mn(pmpa)(H2O)2]) remains relatively rigid as the temperature increases to 90 °C, with only a minor change in its radial distribution functions (RDFs), compared to the RDF peaks at 25 °C. To elucidate the impact of halogens on the magnetic anisotropy of seven-coordinated Mn(II) complexes, we performed both DFT and multireference calculations. This shows that the zero-field splitting (ZFS) parameter D follows the order D(I)> D(Br)> D(Cl). We observed a significant increase in the D-value following the substitution of soft Se-donors in the equatorial position and heavier halogens in the axial position. The D-value of halogen derivatives of Se-analogues varies in the order of D(Cl) < D(I) < D(Br), deviating from the regular spectrochemical series with the discrepancy between the covalency of the Mn(II)-Se bond and the ligand field strength. We anticipate that this study will enhance our understanding of the solvent phase dynamics and structural aspects of ZFS in various Mn(II) complexes with different electronic environments.

Anisotropy and XKS splitting from geodynamic models of double subduction: testing the limits of interpretation

SUMMARY The analysis of the splitting signature of XKS phases is crucial for constraining seismic anisotropy patterns, especially in complex subduction settings such as outward-dipping double subduction. A natural example of this is found in the Central Mediterranean, where the Apennine and the Dinaride slabs subduct in opposite directions, with the Adriatic plate separating them. To assess the capability of XKS-splitting analysis in revealing anisotropic seismic properties, such as fast polarization directions and shear wave anisotropy (in per cent), we use three-dimensional numerical geodynamic models combined with texture evolution simulations. In these models, two identical outward-dipping oceanic plates are separated by a continental plate. Using the full elastic tensors – directly derived from the texture evolution simulations – we compute anisotropic seismic properties and synthetic teleseismic waveforms. From these waveforms synthetic observables are determined, including apparent splitting parameters (fast polarization directions and delay times) and splitting intensities. Based on these observables, we (1) derive models for a single anisotropic layer (one-layer model), (2) identify regions with significant depth-dependent anisotropic seismic properties, and (3) perform inversions at selected locations in terms of two anisotropic layers (two-layer model). We consider two geodynamic models: one with a strong (M1) and one with a weak (M2) continental plate. Model M1 exhibits significant retreat of the subducting plates with no horizontal stretching of the continental plate, whereas Model M2 shows less retreat, substantial horizontal stretching, and detachment of the subducting plates. These different subduction styles result in distinct flow and deformation patterns in the upper mantle, which are reflected in the anisotropic seismic properties. In Model M1, the fast polarization directions below the continental plate are predominantly trench-parallel, whereas in Model M2, they are mostly trench-normal. In most regions of both models, the one-layer models are sufficient to resolve the anisotropic seismic properties, as these properties are nearly constant with depth. However, for both models, we identify some isolated regions – primarily near the tips of the subducting plates and beneath the continental plate – where fast polarization directions exhibit significant variations with depth. Inverting the apparent splitting parameters in these regions yields multiple two-layer models at each location that excellently fit the observables. However, their anisotropic seismic properties can vary significantly, and not all these two-layer models adequately approximate the true depth variations. This ambiguity can be partially reduced by selecting two-layer models in which the summed shear wave anisotropy closely matches that of one of the one-layer models, as these models better capture the true variations.

Structure-optical properties of CdO glass matrix with low levels of cobalt impurities

The impacted role of cadmium oxide (CdO) on the structure-optical properties of current borate glasses with low cobalt oxide impurities was studied. The study was performed on a glass system of the formula [x CdO – (80-x) B2O3 – 0.4 CoO – 19.6 Na2O]; x = 0, 6, 12, 18 and 24 mol%, and produced by melt quenching method. Structure and optical absorption spectroscopy were done to explore the electronic transitions of Co cations inside the glass matrix. The obtained data indicated that structure and optical properties were significantly changed following CdO additives. Bands of optical absorption within both visible and near-infrared ranges showed signatures from Co2+ cations in the tetrahedral and octahedral sites, while Co3+ cations in the octahedral sites within the glass network. Optical absorption also changed significantly for the CdO-doped glasses which exhibited red shifts in visible absorption bands and blue shifts in the NIR absorption bands. The obtained absorption bands in both visible and NIR were employed to obtain the crystal field splitting parameter (Δt) showing augmented records from 3364 cm−1 to 3377 cm−1 with more CdO addition. The Tauc gap energy values were then attained showing decreased values from 3.50 to 3.27 eV when further CdO is added. Further, the impacted role of CdO addition on the nonlinearity optical properties of the present glasses was finally estimated. Potential increases in the linear and nonlinear optical are discussed based on polarizability and basicity increments inside the glass matrix.

Magnetic studies and Hirshfeld surface analysis of acetato bridged 2D Mn(II) coordination polymer with 4-aminobenzohydrazide ligand

A new acetato bridged 2D coordination polymer with general formula of [Mn(µ-4-ambh)(µ-OAc)(OAc)]n (1) was prepared by the reaction of Mn(OAc)2·4H2O with 4-aminobenzohydrazide (4-ambh) in equimolar ratio. The reaction was done in methanol under reflux condition and the single crystals of 1 were obtained by evaporation of the solvent. The structure of this compound was determined by X-ray analysis, and its spectroscopic properties and thermal stability were also studied. These studies indicated that in this compound, the manganese ions are connected by a mixture of acetate and 4-ambh bridges to form a 2D polymeric network at the bc crystallographic plane. The intermolecular and intramolecular interactions were investigated by Hirshfeld surface analysis. Magnetic studies revealed a dominant antiferromagnetic intermetallic coupling of J = −0.78 cm−1 magnitude. Magnetization measurements suggest a diamagnetic ground state and a small gap of 1.59 cm−1 with the first excited state; while the axial zero-field splitting parameter, D, was calculated at −0.21 cm−1.

Teleseismic measurements of Upper Mantle Shear-Wave Anisotropy in Southern Mexico

The Mexican subduction system is an ideal region to study 3-D mantle deformation patterns in response to changes in slab geometry and the presence of tears. Shear-wave splitting measurements were made using SKS, SKKS, and PKS waves in southern Mexico, where the Cocos slab subducts beneath the North American and western Caribbean plates. For most of southern Mexico, the results are consistent with predominantly trench-normal fast polarization directions that can be interpreted as a consequence of sub-slab entrained flow and 2-D corner flow in the mantle wedge in the presence of A-type olivine fabric (or similar). This pattern of trench-perpendicular fast axes extends northward to the region southeast of the Trans-Mexican Volcanic Belt. Beneath its eastern end, fast axes rotate ∼20° clockwise and are likely controlled by the absolute motion of the North American plate. In southeastern Mexico, along the coast and above the mantle wedge tip, the fast axes are trench-normal and the delay times are the shortest. They were interpreted to result from a possibly serpentinized mantle wedge tip. In the same region above the mantle wedge core, the splitting parameters appear to result from different flow patterns in the mantle wedge and the sub-slab mantle.

Impact of SnO2 doping on spectroscopic characteristics of the SnO2-Na2O-CoO-B2O3 borate glass matrix

The role of tin oxide (SnO2) on the structure-optical features of borate glasses containing low cobalt oxide impurities was studied. The study was carried out on a glass system [x SnO2 – (80-x) B2O3 – 19.5 Na2O – 0.5 CoO]; x = 0.0, 0.5, 1.0, 2.0 and 3.0 mol%. Glass synthesis was carried out by the melt-quench technique. Structure and optical absorption spectroscopy were carried out to explore the optical transitions of Co cations inside the present glass matrix. The obtained results indicate that, with the further addition of SnO2 content, the molar volume decreased from 31.675 cm3/mol to 29.147 cm3/mol. FTIR studies indicated the basic structural units of trigonal BO3 units and BO4 tetrahedra. Herein also, additional contents of tin oxide increased the nonbridging oxygen contents. The reason behind such increment is attributed to the existence of Sn4+ ions in octahedral coordination acting as network modifiers. Optical absorption spectra indicated the existence of Co ions in both Co2+ and Co3+ oxidation states. Moreover, the existence of significant SnO2 concurred with the significant blue shift in the visible absorption bands and red shift in near infrared absorption bands. The positions absorption bands of Co ions were analyzed in the context of ligand field parameters determinations. The obtained data of the crystal field splitting parameter (10Dqt) showed decreased values from 3451 cm−1 to 3364 cm−1 with more SnO2 addition. While Racah parameter B values were also, attained showing increased values from 967 cm−1 to 985 cm−1 when additional SnO2 is added. Further, the impacted role of SnO2 addition on the structural and optical properties of the present glasses was finally estimated.

Lateral and vertical variations of seismic anisotropy in the lithosphere–asthenosphere system underneath Central Europe from long-term splitting measurements

SUMMARY Backazimuthal variations in the shear wave splitting of core-refracted shear waves (SKS, SKKS and PKS phases, jointly referred to as XKS) at the Black Forest Observatory (BFO, Southwest Germany) indicate small-scale lateral and partly vertical variations of the seismic anisotropy. However, existing anisotropy studies and models for the nearby Upper Rhine Graben (URG) area in the northern Alpine foreland are mostly based on short-term recordings and by this suffer from a limited backazimuthal coverage and averaging over a wide or the whole backazimuth range. To identify and delimit laterally confined anisotropy regimes in this region, we carry out XKS splitting measurements at six neighbouring (semi-)permanent broad-band seismological recording stations (interstation distance 10–80 km). We manually analyse long-term (partly &amp;gt; 20 yr) recordings to achieve a sufficient backazimuthal coverage to resolve complex anisotropy. The splitting parameters (fast polarization direction $\phi $, delay time $\delta t$) are determined in a single- and multi-event analysis. We test structural anisotropy models with one layer with horizontal or tilted symmetry axis and with two layers with horizontal symmetry axes (transverse isotropy). To account for lateral variations around a single recording site, modelling is compared for the whole and for limited backazimuth ranges. Based on this, we provide a 3-D block model with spatial variation of anisotropic properties. Based on delay times &amp;gt; 0.3 s and missing discrepancies between SKS and SKKS phases, which do not support lower mantle anisotropy, the found anisotropy is placed in the lithosphere and asthenosphere. The spatial distribution as well as the lateral and backazimuthal variations of the splitting parameters confirm lateral and partly vertical variations in anisotropy. On the east side of the URG, we suggest two anisotropic layers in the Moldanubian Zone (south) and one anisotropic layer in the Saxothuringian Zone (north). In the Moldanubian Zone, a change of the fast polarization directions is observed between the east and the west side of the URG, indicating different textures. At the boundary between the two terranes, an inclined anisotropy is modelled which may be related with deformation during Variscan subduction. Regarding the observation of numerous null measurements and inconsistent splitting parameters, especially (southwest of BFO) in the southern URG, different hypothesis are tested: scattering of the seismic wavefield due to small-scale lateral heterogeneities, a vertical a-axis due to a vertical mantle flow related to the Kaiserstuhl Volcanic Complex, as well as a different preferred orientation of the olivine crystals (not A-type, but C-type) due to specific ambient conditions (high temperature, water content).

Modeling of Cr3+ doped Cassiterite (SnO2) Single Crystals

Using the superposition model, the crystal field and zero-field splitting parameters of Cr3+ doped cassiterite (tin oxide), SnO2 single crystals are computed. For calculations, the appropriate locations for Cr3+ ions in SnO2 with distortion are taken into account. The experimental values and the zero-field splitting parameters in theory with local distortion agree fairly well. Using the Crystal Field Analysis Program and crystal field parameters, the optical energy bands for Cr3+ in SnO2 are calculated. The findings indicate that in SnO2 single crystals, one of the Sn4+ ions is replaced by Cr3+ ions.

High-field/high-frequency electron spin resonances of Fe-doped β−Ga2O3 by terahertz generalized ellipsometry: Monoclinic symmetry effects

We demonstrate detection and measurement of electron paramagnetic spin resonances (EPR) of iron defects in β−Ga2O3 utilizing generalized ellipsometry at frequencies between 110 and 170 GHz. The experiments are performed on an Fe-doped single crystal in a free-beam configuration in reflection at 45∘ and magnetic fields between 3 and 7 T. In contrast with low-field, low-frequency EPR measurements, we observe all five transitions of the s=5/2 high-spin state Fe3+ simultaneously. We confirm that ferric Fe3+ is predominantly found at octahedrally coordinated Ga sites. We obtain the full set of fourth-order monoclinic zero-field splitting parameters for both octahedrally and tetrahedrally coordinated sites by employing measurements at multiple sample azimuth rotations. The capability of high-field EPR allows us to demonstrate that simplified second-order orthorhombic spin Hamiltonians are insufficient, and fourth-order terms as well as consideration of the monoclinic symmetry are needed. These findings are supported by computational approaches based on density-functional theory for second-order and on ligand-field theory for fourth-order parameters of the spin Hamiltonian. Terahertz ellipsometry is a way to measure spin resonances in a cavity-free setup. Its possibility of varying the probe frequency arbitrarily without otherwise changing the experimental setup offers unique means of truly disentangling different components of highly anisotropic spin Hamiltonians. Published by the American Physical Society 2024

Optimal Combination of the Splitting–Linearizing Method to SSOR and SAOR for Solving the System of Nonlinear Equations

The symmetric successive overrelaxation (SSOR) and symmetric accelerated overrelaxation (SAOR) are conventional iterative methods for solving linear equations. In this paper, novel approaches are presented by combining a splitting–linearizing method with SSOR and SAOR for solving a system of nonlinear equations. The nonlinear terms are decomposed at two sides through a splitting parameter, which are linearized around the values at the previous step, obtaining a linear equation system at each iteration step. The optimal values of parameters are determined to minimize the reciprocal of the maximal projection, which are sought in preferred ranges using the golden section search algorithm. Numerical tests assess the performance of the developed methods, namely, the optimal splitting symmetric successive over-relaxation (OSSSOR), and the optimal splitting symmetric accelerated over-relaxation (OSSAOR). The chief advantages of the proposed methods are that they do not need to compute the inverse matrix at each iteration step, and the computed orders of convergence by OSSSOR and OSSAOR are between 1.5 and 5.61; they, without needing the inner iterations loop, converge very fast with saving CPU time to find the true solution with a high accuracy.

Time Varying Crustal Anisotropy at Whakaari/White Island Volcano

AbstractWhakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear‐wave splitting analyses to investigate variations in anisotropy between 2018 and 2020, during quiescence, unrest, and the eruption. We examined spatial and temporal variations in 3,499 shear‐wave splitting and 2,656 Vp/Vs ratio measurements. Comparing shear‐wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

An iterative weighted least-square fitting method for crustal anisotropy using receiver functions

SUMMARY The harmonic variation of the P-to-S converted phases (i.e. Pms) observed from receiver functions (RFs) includes information on crustal azimuthal anisotropy. However, this harmonic analysis is easily influenced by low-quality RF traces, and the measurements may be misleading. Here, we propose an improved method, named the iterative weighted least-square method (IWLS), to extract the splitting parameters of the crust and simultaneously retrieve the two- and four-lobed components of backazimuthal variation. The quality and weights of different RF traces are estimated properly in the IWLS method. The weight function is related to the sharpness of the Pms phase and the smearing of other signals. We conduct many synthetic tests, and the IWLS method provides stable measurements for poor backazimuthal coverage, strong noise, weak P-wave azimuthal anisotropy and multiple anisotropic layers. We apply the IWLS method to observational data from two temporary stations on the southeastern Tibetan Plateau and North China Craton, respectively. The measurements are comparable to previous results and provide insight into crustal deformation.

Axial Coordinated Manganese(III) Porphyrin/Tetraazaporphyrin - 4-(10-phenylanthracen-9-yl)Pyridine Dyads: Self-Assembly, Structure and Spectral Properties in Ground and Excited States.

Self-assembly of new donor-acceptor systems based on (5,10,15,20-tetraphenylporphinato)manganese(III)/(5,10,15,20-tetra-4-tert-butylphenylporphinato)manganese(III)/(octakis(4-tert-butylphenyl)tetraazaporphinato)manganese(III) acetate ((AcO)MnTPP/(AcO)MnTBPP/(AcO)MnTAP) and 4-(10-phenylanthracen-9-yl)pyridine (PyAn) was studied using fluorescence spectroscopy and mass spectrometry. It was found that the coordination complexes of 1 : 1 composition (dyads) are formed in toluene. The spectral properties, the chemical structures and redox behavior of the dyads were described using 1H NMR, IR, ESR spectroscopy and cyclic voltammetry, respectively. The dynamic processes and the characteristics in the excited state of the dyads were obtained using the femtosecond transient absorption spectroscopy method. Density functional theory (DFT), time-dependent DFT methods were used to elucidate the dyad electronic structures and to establish the differences in their frontier molecular orbitals. The analysis of the lambda parameter and the distance of hole-pair interaction was indicated more favorable charge transfer between the macrocycle and the axial PyAn fragment in (AcO)(PyAn)MnTAP. The calculated values of the zero-field splitting parameters D and E/D, together with the g tensors of the lowest spin-orbit state for (AcO)MnTPP and (AcO)(PyAn)MnTPP were obtained using the combination of DFT and Multireference Perturbation Theory (CASSCF/NEVPT2) simulations. The data obtained develop the fundamental basis in the field of photovoltaics and show the prospects for the study of molecular systems of this class.

Discriminating Lithospheric and Asthenospheric Anisotropy Beneath Northern Oman: Sharp Contrast Observed at the Semail Gap Fault Zone

AbstractTo gain a deeper understanding of the extensive and varied lithospheric deformations beneath northern Oman, we examine seismic anisotropy in this region using splitting analysis of teleseismic shear wave data. Our study utilizes data from a dense network consisting of 13 permanent and 45 temporary seismic stations, which were operational for approximately 2.5 years starting from October 2013. By examining the azimuthal distribution of shear wave splitting (SWS) parameters, we were able to divide the study area into three sub‐regions. The stations located to the west of the Hawasina window exhibit relatively azimuthally invariant SWS parameters suggesting a single anisotropic layer. On the other hand, most of the stations located in the central and eastern regions display variations versus back‐azimuth, indicating the potential presence of depth‐dependent anisotropy. The General NW‐SE trend of the Fast Polarization Directions (FPDs) of the one‐layer anisotropy in the west and FPDs of the upper layers in the east is concordant with the strike of the structures resulting from the collision between the continental and oceanic plates. A clear contrast in SWS parameters is observed in the Semail Gap Fault Zone (SGFZ), suggesting that the SGFZ can be a lithospheric‐scale structure that hampers the intrusion of mafic magma from the southeast. Furthermore, the FPDs of the lower layer in the east exhibit an NE‐SW trend, which may be indicative of the large‐scale mantle flow resulting from the present‐day plate motion.

Understanding the multiplet fine structures of MgO: Cr3+ with combination method of first−principles calculations and ligand field theory

The interpretative and predictive capabilities on the energy fine structures of transition metal ions in solids are explored by integrating the first -principles calculations with ligand-field (LF) theory. The charge-compensated site of MgO: Cr3+ is chosen as a case and the aim is to determine the energy-level fine structure and EPR g factor g// and g⊥. The approach involves initially computing all the required parameters for LF theory through first-principles calculations, and these parameters are subsequently fed into the LF model Hamiltonian, which fully incorporates interactions between different spectrum terms through the LF operator, spin-orbit coupling (SOC) operator and other relevant factors. The calculated zero-field splitting parameter D and g factors of the ground term A24(t23) are reasonably good agreement with the experimental observations. However, the energy levels of excited states are slightly underestimated compared to experimental values due to the underestimation of Racah parameters in the first-principles calculations, except the T24[t22(T13)e] term which is insensitive to Racah parameters for it is equal to 10Dq relative to ground term A24(t23) with possible minor corrections due to perturbations.

Estimation of Frequencies and Splitting of &lt;sub&gt;2&lt;/sub&gt;S&lt;sub&gt;1&lt;/sub&gt; and &lt;sub&gt;3&lt;/sub&gt;S&lt;sub&gt;1&lt;/sub&gt; Modes from the IGETS Superconducting Gravimeter Network Data after the 2013 Sea of Okhotsk Earthquake

For the first time, the 2S1 and 3S1 modes of the Earth’s free oscillations following the largest deep-focus earthquake in the Sea of Okhotsk on May 24, 2013, have been analyzed using the maximum likelihood method. Data from 14 superconducting gravimeters of the IGETS network deployed in the central part of Europe were used for the analysis. Estimates of the time of mode excitation after the earthquake origin time, degenerate frequencies, and mode splitting parameters are obtained.