Magnetic field models for HD 116458 and HD 126515

We have modeled the magnetic fields of the slowly rotating stars HD 116458 and HD 126515 using the “magnetic charge” technique. HD 116458 has a small angle between its rotation axis and dipole axis (β = 12°), whereas this angle is large for HD 126515 (β = 86°). Both stars can be described with a decentered-dipole model, with the respective displacements being r = 0.07 and r = 0.24 in units of the stellar radius. The decentered-dipole model is able to satisfactorily explain the phase relations for the effective field, B_e(P), and the mean surface field, B_s(P), for both stars, along with the fact that the B_e(P) phase relation for HD 126515 is anharmonic. We discuss the role of systematic measurement errors possibly resulting from instrumental or methodical effects in one or both of the phase relations. The displacement of the dipole probably reflects real asymmetry of the stellar field structure, and is not due to measurement errors. Using both phase relations, B_e(P) and B_s(P), in the modeling considerably reduces the influence of the nonuniform distribution of chemical elements on the stellar surface.     

Effect of 2-D Random Field Discretization on Failure Probability and Failure Mechanism in Probabilistic Slope Stability

This paper investigates, using the random field theory and Monte Carlo simulation, the effects of random field discretization on failure probability, p _f, and failure mechanism of cohesive soil slope stability. The spatial sizes of the discretized elements in random field Δx, Δy in horizontal and vertical directions, respectively, are assigned a series of combinational values in order to model the discretization accuracy. The p _f of deterministic critical slip surface (DCSS) and that of the slope system both are analyzed. The numerical simulation results have demonstrated that both the ratios of Δy/λ _y (λ _y  = scale of fluctuation in vertical direction) and Δx/λ _x (λ _x  = scale of fluctuation in horizontal direction) contribute in a similar manner to the accuracy of p _f of DCSS. The effect of random field discretization on the p _f can be negligible if both the ratios of Δx/λ _x and Δy/λ _y are no greater than 0.1. The normalized discrepancy tends to increase at a linear rate with Δy/λ _y when Δx/λ _x is larger than 0.1, and vice versa for p _f of DCSS. The random field discretization tends to have more considerable influence on the p _f of DCSS than on that of the slope system. The variation of p _f versus λ _x and λ _y may exhibit opposite trends for the cases where the limit state functions of slope failure are defined on DCSS and on the slope system as well. Finally, the p _f of slope system converges in a more rapid manner to that of DCSS than the failure mechanism does to DCSS as the spatial variability of soil property grows from significant to negligible.     

Factors of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O fractionation in corundum estimated from the calculation of isotopic shifts on vibration frequencies

The β-factors of corundum were estimated on the basis of DFT calculations of vibrational frequency changes due to ¹⁶O–¹⁸O isotope substitution in a harmonic approximation using an all-electron Gaussian-type basis set and the B3LYP hybrid functional (the CRYSTAL program). Calculations were performed accounting for eight phonon wave vectors within the first Brillouin zone. The results are approximated by the relation 1000ln β _crn = 9.19874x–0.12326x ² + 0.00213x ³ (x = 10⁶/T(K)², 400 < T(K) < 1300), which can be used in isotope geochemical studies in combination with the known temperature effects on the β-factors of other phases.     

Atmospheric elemental abundances for the components of the multiple system ADS 11061. 41 Draconis

We obtained speckle interferometric and spectroscopic observations of the system 41 Dra during its periastron passage in 2001. The components’ lines are resolved in the spectral interval 3700–9200 Å. The observed wavelength dependence of the brightness difference between the components is used to estimate the B-V indices separately for each of the components: B-V = 0.511 for component a and B-V = 0.502 for component b. We derived improved effective temperatures of the components from their B-V values and hydrogen-line profiles. The observations can be described with the parameters for the components T _eff ^a = 6370 K, log g^a = 4.05 and T _eff ^b = 6410 K, log g^b = 4.20. The iron, carbon, nitrogen, and oxygen abundances in the atmospheres of the components are log N(Fe)^a = 7.55, log N(Fe)^b = 7.60, log N(C)^a = 8.52, log N(C)^b = 8.58, log N(N)^a = 8.05, log N(N)^b = 7.99, log N(O)^a = 8.73, log N(O)^b = 8.76.     

Global solar magnetology and reference points of the solar cycle

The solar cycle can be described as a complex interaction of large-scale/global and local magnetic fields. In general, this approach agrees with the traditional dynamo scheme, although there are numerous discrepancies in the details. Integrated magnetic indices introduced earlier are studied over long time intervals, and the epochs of the main reference points of the solar cycles are refined. A hypothesis proposed earlier concerning global magnetometry and the natural scale of the cycles is verified. Variations of the heliospheric magnetic field are determined by both the integrated photospheric i(B _r )_ph and source surface i(B _r )_ss indices, however, their roles are different. Local fields contribute significantly to the photospheric index determining the total increase in the heliospheric magnetic field. The i(B _r )_ss index (especially the partial index ZO, which is related to the quasi-dipolar field) determines narrow extrema. These integrated indices supply us with a “passport” for reference points, making it possible to identify them precisely. A prominent dip in the integrated indices is clearly visible at the cycle maximum, resulting in the typical double-peak form (the Gnevyshev dip), with the succeeding maximum always being higher than the preceding maximum. At the source surface, this secondary maximum significantly exceeds the primary maximum. Using these index data, we can estimate the progression expected for the 23rd cycle and predict the dates of the ends of the 23rd and 24th cycles (the middle of 2007 and December 2018, respectively).     

Determination of atrazine isotherms on agricultural soils

For determination of atrazine isotherms in agricultural soils of Fars Province, composite soil samples from 0 to 5 cm depth with textures of silty clay loam, clay loam and loam were collected. In order to form the atrazine isotherms, 10, 50 and 100 µg atrazine g^?1 soil was added to the soil samples. Soluble atrazine in water:soil ratios of 10:1, 50:1 and 200:1 was measured after 3-h shaking. Finally, for each cases of applied atrazine, water extractable atrazine was determined and quantified using gas chromatography instrument. The results indicated that there was a linear relationship between the logarithms of water extractable atrazine and added atrazine for different water:soil ratios. A general equation of WEA = K(WS) ^α (AA) ^β is obtained experimentally between water extractable atrazine, µg g^?1(WEA), and added atrazine, µg g^?1 (AA), where K, α and β are absorption constants; WS is the water:soil ratio, g g^?1. For the loam, silty clay loam and clay loam soil textures, the α were 0.49, 0.23 and 0.13, respectively, the β were 0.55, 0.806 and 0.21, respectively, and the K were 1.44, 0.78 and 25.38, respectively.     

Atmospheric chemical abundances of the components of the quadruple system ADS 11061. 40 Draconis

As part of our study of the components of the hierarchic quadruple system ADS 11061, we acquired spectroscopic observations of the binary 40 Dra. Echelle spectra showing the separation of the components’ lines were obtained in the spectral range 3700–9200 Å. Effective temperatures and surface gravities were derived for the components from BV photometry and the hydrogen-line profiles. The components of the 40 Dra system have parameters close to T _eff ^a = 6420 K, log g _a = 4.17, T _eff ^b = 6300 K, and log g _b = 4.20. We find the microturbulence velocity in the component atmospheres to be V _t = 2.6 km/s. The abundances of iron, carbon, nitrogen, and oxygen in the atmospheres of both components are estimated to be log N(Fe)^a = 7.50, log N(Fe)^b = 7.46, log N(C)^a = 8.39, log N(C)^b = 8.45, log N(N)^a = 8.12, log(N)^b = 8.15, log N(O)^a = 8.77, log N(O)^b = 8.74.     

Seyfert nuclei as sources of ultrahigh-energy cosmic rays

Particles can be accelerated to ultrahigh energies E≈10²¹ eV in moderate Seyfert nuclei. This acceleration occurs in shock fronts in relativistic jets. The maximum energy and chemical composition of the accelerated particles depend on the magnetic field in the jet, which is not well known; fields in the range ～5–1000 G are considered in the model. The highest energies of E≈10²¹ eV are acquired by Fe nuclei when the field in the jet is B≈16 G. When B～(5–40) G, nuclei with Z<10 are accelerated to E≤10²⁰ eV, while nuclei with Z≥10 acquire energies E≥2×10²⁰ eV. Only particles with Z≥23 acquire energies E≤10²⁰ eV when B～1000 G. Protons are accelerated to E<4×10¹⁹ eV, and do not fall into the range of energies of interest for any magnetic field B. The particles lose a negligible amount of their energy in interactions with infrared photons in the accretion disk; losses in the thick gas-dust torus are also negligible if the luminosity of the galaxy is L≤10⁴⁶ erg/s and the angle between the normal to the galactic plane and the line of sight is sufficiently small, i.e., if the axial ratio of the galactic disk is comparatively high. The particles do not lose energy to curvature radiation if their deviations from the jet axis do not exceed 0.03–0.04 pc at distances from the center of R≈40–50 pc. Synchrotron losses are small, since the magnetic field frozen in the galactic wind at R≤40–50 pc is directed (as in the jet) primarily in the direction of motion. If the model considered is valid, the detected cosmic-ray protons could be either fragments of Seyfert nuclei or be accelerated in other sources. The jet magnetic fields can be estimated both from direct astronomical observations and from the energy spectrum and chemical composition of cosmic rays.     

VSP wavefield separation using the gray-scale Hough transform: synthetic data

The principal objective in Vertical Seismic Profile (VSP) processing is the separation of the downgoing and upgoing wavefields. Several methods have been suggested in this area. This paper presents a new approach based on the gray-scale Hough transform (GSHT) which is an extension of the conventional Hough transform used to detect straight lines and other curves. The technique, we suggest here, directly maps the gray-scale VSP image, including the downgoing and upgoing linear events, in image coordinate space (x,t,g) to the gray Hough parameter counting space (θ,ρ,g), where θ and ρ are the polar parameters and g is the gray-scale value. In this new space, the downgoing events appear in the negative angles θ quadrant and the upgoing in the positive quadrant, owning to their opposite apparent velocities. The inverse GSHT algorithm, we developed in this study, is performed for extracting separately these two wavefields by considering the straight lines that satisfy the corresponding filtering conditions. The experimental results on synthetic VSP datasets are convincing. The wave separation is well performed, even in the presence of loud noise levels, with signal to noise ratio improvement and amplitude preservation, in contrast to median filtering.     

Masses of the visual components and black holes in X-ray novae: Effects of proximity of the components

It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V _rot sin i and underestimation of the component mass ratio q = M _x /M _v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M _x and visual components M _v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M _x /M _v , reducing M _v , while M _x changes only slightly. The resulting distribution of M _v attains its maximum near M _v ? 0.35M _⊙, in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M _x also differs strongly from the mass distribution for massive stars in the Galaxy.     

Comparison of the magnetic properties of leading and following spots and the overlying ultraviolet emission

SDO/HMI and SDO/AIA data for the 24th solar-activity cycle are analyzed using a quicker and more accurate method for resolving π ambiguities in the transverse component of the photospheric magnetic field, yielding new results and confirming some earlier results on the magnetic properties of leading and following magnetically connected spots and single spots. The minimum inclination of the field lines to the positive normal to the solar surface α _min within umbrae is smaller in leading than in following spots in 78% of the spot pairs considered; the same trend is found for the mean angle 〈α〉 in 83% of the spot pairs. Positive correlations between the α _min values and the 〈α〉 values in leading and following spots are also found. On average, in umbrae, the mean values of 〈B〉, the umbra area S, and the angles α _min and 〈α〉 decrease with growth in the maximum magnetic field B _max in both leading and following spots. The presence of a positive correlation between B _max and S is confirmed, and a positive correlation between 〈B〉 and S in leading and following spots has been found. Themagnetic properties of the umbrae of magnetically connected pairs of spots are compared with the contrast of the He II 304 emission above the umbrae, C ₃₀₄. Spots satisfying certain conditions display a positive correlation between C _304?L and 〈α _L 〉 for the leading (L) spots, and between C _304?L /C _304?F and l _L /l _F , where l _L (l _F ) are the lengths of the field lines connecting leading (L) or following (F) spots from the corresponding spot umbrae to the apex of the field line.     

Water determination in iron oxyhydroxides and iron ores by Karl Fischer titration

Protohematite (Fe_2?x/3(OH) _x O_3?x 1 ≤ x < 0.5) and hydrohematite (Fe_2?x/3(OH) _x O_3?x 0.5 ≤ x < 0) are iron-defective phases containing hydroxyl groups in their structures. These species were described in prior studies mainly with the aid of X-ray diffraction and Infrared spectroscopy. The existence of these phases in soils might have influence in redox processes, and they were considered as a possible water reservoir in Martian soils. In this study, we have used for the first time the Karl Fischer titration method to determine the amount of water released after heating several synthetic samples of goethite, hematite and natural iron ores at 105, 400, 600 and 900 °C. It was found that heating at 105 °C did not remove all moisture from the samples, and higher temperatures were necessary to completely remove all the absorbed water. The water contents determined at 400, 600 and 900 °C were found to be the same within the experimental errors, suggesting the inexistence of both protohematite and hydrohematite in the investigated samples. Therefore, the above-mentioned effects of these phases in soils might have to be reevaluated.     

Electron paramagnetic resonance studies on clinochlore from Longitudinal Valley area,northeastern Taiwan

A natural sample of clinochlore from the Longitudinal Valley area of northeastern Taiwan has been characterized by using the powder X-ray diffraction (XRD), differential thermal analysis and electron paramagnetic resonance (EPR) spectroscopic techniques. The lattice parameters of the monoclinic (IIb) clinochlore with the composition (Mg_2.988 Al_1.196 Fe_1.6845 Mn_0.026)_5.8945 (Si_2.559 Al_1.441)₄ O₁₀ (OH)₈ have been calculated from the powder XRD data and are found to be a = 5.347 Å, b = 9.223 Å, c = 14.250 Å, β = 97.2° and Z = 2. The thermal behaviour of the sample showed the typical behaviour of clinochlore with a hydroxyl content of 12.5 wt%. The EPR spectrum at room temperature exhibits two resonance signals centred at g ≈ 2.0 and g ≈ 8.0. The signal at g ≈ 2.0 shows a six-line hyperfine structure which is a characteristic of Mn²⁺ ions in octahedral symmetry. The resonance signal at g ≈ 8.0 is a characteristic of Fe³⁺ ions. The EPR spectra have also been recorded at different temperatures (123–295 K). The population of spin levels (N) has been calculated for g ≈ 2.0 and g ≈ 8.0 resonance signals. It is observed that N increases with decreasing temperature. From EPR spectra, the spin-Hamiltonian parameters have been evaluated. The zero-field splitting parameter (D) is found to be temperature dependent. The peak-to-peak width of the g ≈ 8.0 resonance signal is found to increase with decrease in temperature.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

Separation of chemical elements in the atmospheres of CP stars under the action of light induced drift

A mechanism for the separation of chemical elements and isotopes in the atmospheres of chemically peculiar (CP) stars due to light-induced drift (LID) of ions is discussed. The efficiency of separation due to LID is proportional to the relative difference of the transport frequencies for collisions of ions of heavy elements located in the excited state (collision frequency ν _e ) and ground state (collision frequency ν _g ) with neutral buffer particles (hydrogen and helium), (ν _e ? ν _g )/ν _g . The known interaction potentials are used to numerically compute the relative difference (ν _e ^H ? ν _g ^H )/νg H for collisions between the ions Be⁺, Mg⁺, Ca⁺, Sr⁺, Cd⁺, Ba⁺, Al⁺, and C⁺ and hydrogen atoms. These computations show that, at the temperatures characteristic of the atmospheres of CP stars, T = 7000?20 000 K, values of |ν _e ^H ?ν _g ^H |/ν _g ^H ≈ 0.1?0.4 are obtained. With such relative differences in the transport collision frequencies, the LID rate of ions in the atmospheres of coolCP stars (T < 10000 K) can reach ~0.1 cm/s,which exceeds the drift rate due to light pressure by an order of magnitude. This means that, under these conditions, the separation of chemical elements under the action of LID of ions could be an order of magnitude more efficient than separation due to light pressure. Roughly the same manifestations of LID and light pressure are also expected in the atmospheres of hotter stars (20 000 > T > 10 000 K). LID of heavy ions is manifest only weakly in very hot stars (T > 20 000 K).     

Generation of solar magnetic fields by the αδΘ dynamo

We consider a solar dynamo mechanism that generates large-scale magnetic fields due to the combined action of cyclonic flows (the α effect), differential rotation (the Θ effect), and the non-uniformity of large-scale magnetic fields (the Θ × J effect). Our results are based on numerical model which takes into account currently available data on the differential rotation of the convection zone and the intensity of convective flows in the solar interior. A reasonable choice of parameters characterizing the intensity of magnetic-field generation by the α and Θ × J mechanisms can account for an oscillatory dynamo regime with properties similar to the 22-year magnetic-activity cycle of the Sun. We analyze the nonlinear saturation of the generation effects in the large-scale magnetic field, due to either magnetic stresses or the conservation of magnetic helicity. Allowance for the helicity of the small-scale magnetic fields is of crucial importance in limiting the energy of the generated large-scale magnetic field.     

The magnetic-field structure in a stationary accretion disk

The magnetic-field structure in regions of stationary, planar accretion disks around active galactic nuclei where general-relativistic effects can be neglected (from 10 to 200 gravitational radii) is considered. It is assumed that the magnetic field in the outer edges of the disk, which forms in the magnetosphere of the central black hole during the creation of the relativisitic jets, corresponds to the field of a magnetic dipole perpendicular to the plane of the disk. In this case, the azimuthal field component B_φ in the disk arises due to the presence of the radial field B_ρ and the azimuthal velocity component U_φ. The value of the magnetic field at the inner radius of the disk is taken to correspond to the solution of the induction equation in a diffusion approximation. Numerical solutions of the induction equation are given for a number of cases.     

Geodynamics and underlying bedrock of the magnetically active crust layer of the Lut block,Eastern Iran

The Curie point depth map of Eastern Iran was constituted from spectral analysis of the aeromagnetic data. The reduction to pole (RTP) was applied to the magnetic anomaly data. The Curie point depth values from 165 overlapping blocks, 100 × 100 km in size, have been estimated. The Curie point depth method provides a relationship between the 2-D FFT power spectrum of the magnetic anomalies and the depth of magnetic sources by transforming the spatial data into the frequency domain. The centroid and top depth of the magnetic sources (respectively Z₀ and Z_t) is calculated from radially averaged log power spectrum for each block. Finally, the Curie point depth of Eastern Iran is obtained by Z_b = 2Z₀–Z_t. The highest value of 24 km is located in eastern and western boundaries of the Lut block, and the lowest value of 12 km is located at north of study area. The shallow depths in the Curie-point depth map are well correlated with the young volcanic areas and geothermal potential fields. Geothermal gradient ranging from 24 to 45°C/km. The deduced thermal structure in eastern Iran has a relationship with orogenic collapse associated with delamination of thickened lithospheric root between the Lut and Afghan continental blocks.     

Structural systematics of Ge substitution in primitive pyroxenes

Pyroxenes of general stoichiometry Mg(Ge _x Si_1?x )O₃ were encountered in attempts to synthesise Ge-substituted talcs at 0.2 GPa, 650–700 °C. Orthopyroxenes (Pbca) of compositions x = 0.21, 0.30, and 0.34 were identified, and also a P2₁/c clinopyroxene of composition x = 0.63, and C2/c clinopyroxenes of compositions x = 0.91 and 1. End-member clinoenstatite MgSiO₃-P2₁/c synthesised at 16 GPa, 1300 °C and transformed from C2/c was also included in the study. Crystal structure refinements using single-crystal XRD data showed that unit-cell parameters vary linearly with Si–Ge for the Pbca and P2₁/c pyroxenes, both of which have two symmetrically non-equivalent tetrahedral chains. Refinement of Si–Ge occupancies at tetrahedral sites showed that the two chains of all primitive pyroxenes have very different compositions, with X_Ge(TB) ? X_Ge(TA). This difference arises from the greater flexibility of the B-chain to rotate in response to tetrahedral expansion due to increasing Ge content. The TA-M2 shared polyhedral edge imposes significant constraints on the flexibility of the A-chain, which can accommodate much less Ge than the B-chain. Linear trends of cell parameters, site occupancies, and structural parameters for the primitive pyroxenes, when extrapolated to published data for MgGeO₃–Pbca, extend across the entire Si–Ge join.     

Pre-flare changes in the turbulence regime for the photospheric magnetic field in a solar active region

Observations of the total magnetic field in the active region NOAA 6757 have been used to study the turbulence regime from 2.5 h before the onset of a 2B/X1.5 flare until two minutes after its maximum. The curvature of the exponent ζ(q) for the structure functions of the B _z field increases monotonically before the flare (i.e., the multifractal character of the B _z field becomes more complex) but straightens at the flare maximum and coincides with a linear Kolmogorov dependence (implying a monofractal structure for the B _z field). The observed deviations of ζ(q) from a Kolmogorov line can be used for short-term forecasting of strong flares. Analysis of the power spectra of the B _z field and the dissipation of magnetic-energy fluctuations shows that the beginning of the flare is associated with the onset of a new turbulence regime, which is closer to a classical Kolmogorov regime. The scaling parameter (cancellation index) of the current helicity of the magnetic field, k _h , remains at a high level right up until the last recording of the field just before the flare but decreases considerably at the flare maximum. The variations detected in the statistical characteristics of the turbulence can be explained by the formation and amplification of small-scale flux tubes with strong fields before the flare. The dissipation of magnetic energy before the flare is primarily due to reconnection at tangential discontinuities of the field, while the dissipation after the flare maximum is due to the anomalous plasma resistance. Thus, the flare represents an avalanche dissipation of tangential discontinuities.