Transition Radiation (TR) plays an important role in particle identification in high-energy physics and its characteristics provide a feasible method of energy calibration in the energy range up to 10 TeV, which is of interest for dark matter searches in cosmic rays. In a Transition Radiation Detector (TRD), the TR signal is superimposed onto the ionization energy loss signal induced by incident charged particles. In order to make the TR signal stand out from the background of ionization energy loss in a significant way, we optimized both the radiators and the detector. We have designed a new prototype of regular radiator optimized for a maximal TR photon yield, combined with the Side-On TRD which is supposed to improve the detection efficiency of TR. We started a test beam experiment with the Side-On TRD at Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire (CERN), and found that the experimental data is consistent with the simulation results.

Xenon dual-phase time projection chambers designed to search for Weakly Interacting Massive Particles have so far shown a relative energy resolution which degrades with energy above $\sim$200 keV. This has limited their sensitivity in the search for rare events like the neutrinoless double-beta decay of $^{136}$Xe at its $Q$-value, $Q_{\beta\beta}\simeq$ 2.46 MeV. For the XENON1T dual-phase time projection chamber, we demonstrate that the relative energy resolution at 1 $\sigma/\mu$ is as low as (0.79$\pm$0.02) % in its one-ton fiducial mass, and for single-site interactions at $Q_{\beta\beta}$. We achieve this by using a signal correction method to rectify the saturation effects of the signal readout system. The very good energy resolution from keV to MeV energies demonstrated in XENON1T opens up new windows for the xenon dual-phase dark matter detectors to simultaneously search for other rare events.

In the early days of 2017 September, an exceptionally energetic solar active region AR12673 aroused great interest in the solar physics community. It produced four X class flares, more than 20 CMEs and an intense geomagnetic storm, for which the peak value of the Dst index reached up to -142nT at 2017 September 8 02:00 UT. In this work, we check the interplanetary and solar source of this intense geomagnetic storm. We find that this geomagnetic storm was mainly caused by a shock-ICME complex structure, which was formed by a shock driven by the 2017 September 6 CME propagating into a previous ICME which was the interplanetary counterpart of the 2017 September 4 CME. To better understand the role of this structure, we conduct the quantitative analysis about the enhancement of ICME's geoeffectiveness induced by the shock compression. The analysis shows that the shock compression enhanced the intensity of this geomagnetic storm by a factor of two. Without shock compression, there would be only a moderate geomagnetic storm with a peak Dst value of -79 nT. In addition, the analysis of the proton flux signature inside the shock-ICME complex structure shows that this structure also enhanced the solar energetic particles (SEPs) intensity by a factor of ~ 5. These findings illustrate that the shock-ICME complex structure is a very important factor in solar physics study and space weather forecast.

Magnetic flux rope (MFR) is the core structure of the greatest eruptions, i.e., the coronal mass ejections (CMEs), on the Sun, and magnetic clouds are post-eruption MFRs in interplanetary space. There is a strong debate about whether or not a MFR exists prior to a CME and how the MFR forms/grows through magnetic reconnection during the eruption. Here we report a rare event, in which a magnetic cloud was observed sequentially by four spacecraft near Mercury, Venus, Earth and Mars, respectively. With the aids of a uniform-twist flux rope model and a newly developed method that can recover a shock-compressed structure, we find that the axial magnetic flux and helicity of the magnetic cloud decreased when it propagated outward but the twist increased. Our analysis suggests that the `pancaking' effect and `erosion' effect may jointly cause such variations. The significance of the `pancaking' effect is difficult to be estimated, but the signature of the erosion can be found as the imbalance of the azimuthal flux of the cloud. The latter implies that the magnetic cloud was eroded significantly leaving its inner core exposed to the solar wind at far distance. The increase of the twist together with the presence of the erosion effect suggests that the post-eruption MFR may have a high-twist core enveloped by a less-twisted outer shell. These results pose a great challenge to the current understanding on the solar eruptions as well as the formation and instability of MFRs.

The spectrum of the MeV tail detected in the black-hole candidate Cygnus X-1 remains controversial as it appeared much harder when observed with the INTEGRAL Imager IBIS than with the INTEGRAL spectrometer SPI or CGRO. We present an independent analysis of the spectra of Cygnus X-1 observed by IBIS in the hard and soft states. We developed a new analysis software for the PICsIT detector layer and for the Compton mode data of the IBIS instrument and calibrated the idiosyncrasies of the PICsIT front-end electronics. The spectra of Cygnus X-1 obtained for the hard and soft states with the INTEGRAL imager IBIS are compatible with those obtained with the INTEGRAL spectrometer SPI, with CGRO, and with the models that attribute the MeV hard tail either to hybrid thermal/non-thermal Comptonisation or to synchrotron emission.

As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.

Closure quantities measured by very long baseline interferometry (VLBI) observations are independent of instrumental and propagation instabilities and antenna gain factors, but are sensitive to source structure. A new method is proposed to calculate a structure index based on the median values of closure quantities rather than the brightness distribution of a source. The results are comparable to structure indices based on imaging observations at other epochs and demonstrate the flexibility of deriving structure indices from exactly the same observations as used for geodetic analysis and without imaging analysis. A three-component model for the structure of source 3C371 is developed by model-fitting closure phases. It provides a real case of tracing how the structure effect identified by closure phases in the same observations as the delay observables affects the geodetic analysis, and investigating which geodetic parameters are corrupted to what extent by the structure effect. Using the resulting structure correction based on the three-component model of source 3C371, two solutions, with and without correcting the structure effect, are made. With corrections, the overall rms of this source is reduced by 1 ps, and the impacts of the structure effect introduced by this single source are up to 1.4 mm on station positions and up to 4.4 microarcseconds on Earth orientation parameters. This study is considered as a starting point for handling the source structure effect on geodetic VLBI from geodetic sessions themselves.

Gamma-ray polarimetry is a new powerful tool to study the processes responsible for the emission from astrophysical sources and the environments in which this emission takes place. Few successful polarimetric measurements have however been performed thus far in the gamma-ray energy band due to the difficulties involved. POLAR is a dedicated polarimeter designed to perform high precision measurements of the polarization of the emission from gamma-ray burst in the 50-500 keV energy range. This new polarimeter is expected to detect approximately 50 gamma-ray bursts per year while performing high precision polarization measurements on approximately 10 bursts per year. The instrument was launched into lower earth orbit as part of the second Chinese space lab, the Tiangong-2, on September 15th 2016 and has been taking data successfully since being switched on one week after. The instrument uses a segmented scintillator array consisting of 1600 plastic scintillator bars, read out by 25 flat-panel multi-anode photomultipliers, to measure the Compton scattering angles of incoming photons. The small segmentation and relatively large uniform effective area allow the instrument to measure the polarization of a large number of transient events, such as gamma-ray bursts, with an unprecedented precision during its two year life-time. The final flight model underwent detailed calibration prior to launch as well as intensive space qualification tests, a summary of which will be presented in this paper. The instrument design will be discussed first followed by an overview of the on-ground tests, finally the in-orbit behavior as measured during the first weeks of the mission will be presented.

Modified gravity theories have the potential of explaining the recent acceleration of the Universe without resorting to the mysterious concept of dark energy. In particular, it has been pointed out that matter-geometry coupling may be responsible for the recent cosmological dynamics of the Universe, and matter itself may play a more fundamental role in the description of the gravitational processes that usually assumed. We study the quantum cosmology of the $f(R,T)$ gravity theory, in which the effective Lagrangian of the gravitational field is given by an arbitrary function of the Ricci scalar, and the trace of the matter energy-momentum tensor, respectively. For the background geometry we adopt the Friedmann--Robertson--Walker metric, and we assume that matter content of the Universe consists of a perfect fluid. We obtain the general form of the gravitational Hamiltonian, of the quantum potential, and of the canonical momenta, respectively. This allows us to formulate the full Wheeler-de Witt equation describing the quantum properties of this modified gravity model. As a specific application we consider in detail the quantum cosmology of the $f(R,T)=F^0(R)+\theta RT$ model, in which $F^0(R)$ is an arbitrary function of the Ricci scalar, and $\theta $ is a function of the scale factor only. The Hamiltonian form of the equations of motion, and the Wheeler-de Witt equations are obtained, and a time parameter for the corresponding dynamical system is identified, which allows to formulate the Schr\"{o}dinger--Wheeler--de Witt equation for the quantum-mechanical description of the model under consideration. A perturbative approach for the study of this equation is developed, and the energy levels of the Universe are obtained by using a twofold degenerate perturbation approach. A second quantization approach for the description of quantum time is also proposed, and briefly discussed.

The CONT14 campaign with state-of-the-art VLBI data has observed the source 0642+449 with about one thousand observables each day during a continuous observing period of fifteen days, providing tens of thousands of closure delays---the sum of the delays around a closed loop of baselines. The closure delay is independent of the instrumental and propagation delays and provides valuable additional information about the source structure. We demonstrate the use of this new "observable" for the determination of the structure in the radio source 0642+449. This source, as one of the defining sources in the second realization of the International Celestial Reference Frame (ICRF2), is found to have two point-like components with a relative position offset of -426 microarcseconds in right ascension and -66 microarcseconds in declination. The two components are almost equally bright with a flux-density ratio of 0.92. The standard deviation of closure delays for source 0642+449 was reduced from 139 ps to 90 ps by using this two-component model. Closure delays larger than one nanosecond are found to be related to the source structure, demonstrating that structure effects for a source with this simple structure could be up to tens of nanoseconds. The method described in this paper does not rely on a priori source structure information, such as knowledge of source structure determined from direct (Fourier) imaging of the same observations or observations at other epochs. We anticipate our study to be a starting point for more effective determination of the structure effect in VLBI observations.

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

In recent years, more and more gamma-ray bursts with late rebrightenings in multi-band afterglows unveil the late-time activities of the central engines. GRB 100814A is a special one among the well-sampled events, with complex temporal and spectral evolution. The single power-law shallow decay index of the optical light curve observed by GROND between 640 s and 10 ks is $\alpha_{\rm opt} = 0.57 \pm 0.02$, which apparently conflicts with the simple external shock model expectation. Especially, there is a remarkable rebrightening in the optical to near infrared bands at late time, challenging the external shock model with synchrotron emission coming from the interaction of the blast wave with the surrounding interstellar medium. In this paper, we invoke a magnetar with spin evolution to explain the complex multi-band afterglow emission of GRB 100814A. The initial shallow decay phase in optical bands and the plateau in X-ray can be explained as due to energy injection from a spin-down magnetar. At late time, with the falling of materials from the fall-back disk onto the central object of the burster, angular momentum of the accreted materials is transferred to the magnetar, which leads to a spin-up process. As a result, the magnetic dipole radiation luminosity will increase, resulting in the significant rebrightening of the optical afterglow. It is shown that the observed multi-band afterglow emission can be well reproduced by the model.

The High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station is planned for operation starting around 2020 for about 10 years. It is designed as a next generation space facility focused on indirect dark matter search, precise cosmic ray spectrum and composition measurements up to the knee energy, and high energy gamma-ray monitoring and survey. The calorimeter plays an essential role in the main scientific objectives of HERD. A 3-D cubic calorimeter filled with high granularity crystals as active material is a very promising choice for the calorimeter. HERD is mainly composed of a 3-D calorimeter (CALO) surrounded by silicon trackers (TK) from all five sides except the bottom. CALO is made of 9261 cubes of LYSO crystals, corresponding to about 55 radiation lengths and 3 nuclear interaction lengths, respectively. Here the simulation results of the performance of CALO with GEANT4 and FLUKA are presented: 1) the total absorption CALO and its absorption depth for precise energy measurements (energy resolution: 1\% for electrons and gamma-rays beyond 100 GeV, 20\% for protons from 100 GeV to 1 PeV); 2) its granularity for particle identification (electron/proton separation power better than $10^{-5}$); 3) the homogenous geometry for detecting particles arriving from every unblocked direction for large effective geometrical factor ($>$3 ${\rm m}^{2}{\rm sr}$ for electron and diffuse gamma-rays, $>$2 $ {\rm m}^{2}{\rm sr}$ for cosmic ray nuclei); 4) expected observational results such as gamma-ray line spectrum from dark matter annihilation and spectrum measurement of various cosmic ray chemical components.

Well-sampled optical and radio light curves of BL Lacertae in B, V, R, I bands and 4.8, 8.0, 14.5 GHz from 1968 to 2014 were presented in this paper. A possible $1.26 \pm 0.05$ yr period in optical bands and a $7.50 \pm 0.15$ yr period in radio bands were detected based on discrete correlation function, structure function as well as Jurkevich method. Correlations among different bands were also analyzed and no reliable time delay was found between optical bands. Very weak correlations were detected between V band and radio bands. However, in radio bands the variation at low frequency lagged that at high frequency obviously. The spectrum of BL Lacertae turned mildly bluer when the object turned brighter, and stronger bluer-when-brighter trends were found for short flares. A scenario including a precessing helical jet and periodic shocks was put forward to interpret the variation characteristics of BL Lacertae.

There exist cislunar and trans-lunar libration points near the Moon, which are referred as the LL1 and LL2 points respectively and can generate the different types of low-energy trajectories transferring from Earth to Moon. The time-dependent analytic model including the gravitational forces from the Sun, Earth and Moon is employed to investigate the energy-minimal and practical transfer trajectories. However, different from the circular restricted three-body problem, the equivalent gravitational equilibria are defined according to the geometry of instantaneous Hill's boundary due to the gravitational perturbation from the Sun. The relationship between the altitudes of periapsis and eccentricities is achieved from the Poincar\'e mapping for all the lunar captured trajectories, which presents the statistical feature of the fuel cost and captured orbital elements rather than generating a specified Moon-captured segment. The minimum energy required by the captured trajectory on a lunar circular orbit is deduced in the spatial bi-circular model. It is presented that the asymptotical behaviors of invariant manifolds approaching to/from the libration points or halo orbits are destroyed by the solar perturbation. In fact, the energy-minimal cislunar transfer trajectory is acquired by transiting LL1 point, while the energy-minimal trans-lunar transfer trajectory is obtained by transiting LL2 point. Finally, the transfer opportunities for the practical trajectories escaped from the Earth and captured by the Moon are yielded by transiting halo orbits near LL1 and LL2 points, which can be used to generate the whole trajectories.

We show that the photospheres of "failed" Gamma-Ray Bursts (GRBs), whose bulk Lorentz factors are much lower than 100, can be outside of internal shocks. The resulting radiation from the photospheres is thermal and bright in UV/Soft X-ray band. The photospheric emission lasts for about one thousand seconds with luminosity about several times 10^46 erg/s. These events can be observed by current and future satellites. It is also shown that the afterglows of failed GRBs are peculiar at the early stage, which makes it possible to distinguish failed GRBs from ordinary GRBs and beaming-induced orphan afterglows.

Gamma ray bursts (GRBs) have great advantages for their huge burst energies, luminosities and high redshifts in probing the Universe. A few interesting luminosity correlations of GRBs have been used to test cosmology models. Especially, for a subsample of long GRBs with known redshifts and a plateau phase in the afterglow, a correlation between the end time of the plateau phase (in the GRB rest frame) and the corresponding X-ray luminosity has been found. In this paper, we re-analyze the subsample and found that a significantly tighter correlation exists when we add a third parameter, i.e. the isotropic $\gamma$-ray energy release, into the consideration. Additionally, both long and intermediate duration GRBs are consistent with the same three-parameter correlation equation. It is argued that the new three-parameter correlation is consistent with the hypothesis that the subsample of GRBs with a plateau phase in the afterglow be associated with the birth of rapidly rotating magnetars, and that the plateau be due to the continuous energy-injection from the magnetar. It is suggested that the newly born millisecond magnetars associated with GRBs might provide a good standard candle in the Universe.

We propose an off-axis relativistic jet model for the Type Ic supernova SN 2007gr. Most of the energy ($\sim2\times10^{51}$ erg) in the explosion is contained in non-relativistic ejecta which produces the supernova. The optical emission is coming from the decay process of $\rm ^{56}Ni$ synthesized in the bulk SN ejecta. Only very little energy ($\sim10^{48}$ erg) is contained in the relativistic jet with initial velocity about 0.94 times the speed of light. The radio and X-ray emission comes from this relativistic jet. With some typical parameters of a Wolf-Rayet star (progenitor of Type Ic SN), i.e., the mass loss rate $\dot{M}=1.0 \times10^{-5} M_{\odot} \rm yr^{-1}$ and the wind velocity $v_{\rm w}=1.5\times10^{3} \rm km s^{-1}$ together with an observing angle of $\theta_{\rm obs} = 63.3^{\circ}$, we can obtain the multiband light curves that fit the observations well. All the observed data are consistent with our model. Thus we conclude that SN 2007gr contains a weak relativistic jet and we are observing the jet from off-axis.

Early rebrightenings at a post-burst time of 10^2 - 10^4 s have been observed in the afterglows of some gamma-ray bursts (GRBs). Unlike X-ray flares, these rebrightenings usually last for a relatively long period. The continuous energy injection mechanism usually can only produce a plateau in the afterglow light curve, but not a rebrightening. Also, a sudden energy injection can give birth to a rebrightening, but the rebrightening is a bit too rapid. Here we argued that the early rebrightenings can be produced by the ring-shaped jet model. In this scenario, the GRB outflow is not a full cone, but a centrally hollowed ring. Assuming that the line of sight is on the central symmetry axis of the hollow cone, we calculate the overall dynamical evolution of the outflows and educe the multiband afterglow light curves. It is found that the early rebrightenings observed in the afterglows of a few GRBs, such as GRBs 051016B, 060109, 070103 and 070208 etc, can be well explained in this framework. It is suggested that these long-lasted early rebrightenings in GRB afterglows should be resulted from ring-shaped jets.

A long plateau phase and an amazing brightness have been observed in the Xray afterglow of GRB 060729. This peculiar light curve is likely due to long-term energy injection in external shock. Here we present a detailed numerical study on the energy injection process of magnetic dipole radiation from a strongly magnetized millisecond pulsar and model the multi-band afterglow observations. It is found that this model can successfully explain the long plateaus in the observed X-ray and optical afterglow light curves. The sharp break following the plateaus should be due to the rapid decline of the emission power of the central pulsar. At an even late time (~5*10^6s), an obvious jet break appears, which implies a relatively large half opening angle of theta~0.3 for the GRB ejecta. Due to the energy injection, the Lorentz factor of the outflow is still larger than two 10^7s post the GRB trigger, making the X-ray afterglow of this burst detectable by Chandra even 642 days after the burst.

The Galactic hard X-ray transient SWIFT J195509+261406 was first observed as a gamma-ray burst GRB 070610. Within 3 days after the burst, more than forty optical flares had been observed. Here we propose that this peculiar event should be associated with a white dwarf. The hard X-ray burst itself may be triggered by the collision between two planets orbiting the white dwarf. Some cracked frag-ments produced in the collision then fall onto the surface of the white dwarf in several days, giving birth to the observed optical flares via cyclotron radiation. Our model can satisfactorily explain the basic features of the observations.

The conventional magnetic induction equation that governs hydromagnetic dynamo action is transformed into an equivalent integral equation system. An advantage of this approach is that the computational domain is restricted to the region occupied by the electrically conducting fluid and to its boundary. This integral equation approach is first employed to simulate kinematic dynamos excited by Beltrami-like flows in a finite cylinder. The impact of externally added layers around the cylinder on the onset of dynamo actions is investigated. Then it is applied to simulate dynamo experiments within cylindrical geometry including the von Karman sodium (VKS) experiment and the Riga dynamo experiment. A modified version of this approach is utilized to investigate magnetic induction effects under the influence of externally applied magnetic fields which is also important to measure the proximity of a given dynamo facility to the self-excitation threshold.

The asymmetric shape of reversals of the Earth's magnetic field indicates a possible connection with relaxation oscillations as they were early discussed by van der Pol. A simple mean-field dynamo model with a spherically symmetric $\alpha$ coefficient is analysed with view on this similarity, and a comparison of the time series and the phase space trajectories with those of paleomagnetic measurements is carried out. For highly supercritical dynamos a very good agreement with the data is achieved. Deviations of numerical reversal sequences from Poisson statistics are analysed and compared with paleomagnetic data. The role of the inner core is discussed in a spectral theoretical context and arguments and numerical evidence is compiled that the growth of the inner core might be important for the long term changes of the reversal rate and the occurrence of superchrons.

The intention of the ''von Karman sodium'' (VKS) experiment is to study the hydromagnetic dynamo effect in a highly turbulent and unconstrained flow. Much effort has been devoted to the optimization of the mean flow and the lateral boundary conditions in order to minimize the critical magnetic Reynolds number and hence the necessary motor power. The main focus of this paper lies on the role of ''lid layers'', i.e. layers of liquid sodium between the impellers and the end walls of the cylinder. First, we study an analytical test flow to show that lid layers can have an ambivalent effect on the efficiency of the dynamo. The critical magnetic Reynolds number shows a flat minimum for a small lid layer thickness, but increases for thicker layers. For the actual VKS geometry it is shown that static lid layers yield a moderate increase of the critical magnetic Reynolds number by approximately 12 per cent. A more dramatic increase by 100 until 150 per cent can occur when some rotational flow is taken into account in those layers. Possible solutions of this problem are discussed for the real dynamo facility.

In a recent paper (Phys. Rev. Lett. 94 (2005), 184506; physics/0411050) it was shown that a simple mean-field dynamo model with a spherically symmetric helical turbulence parameter alpha can exhibit a number of features which are typical for Earth's magnetic field reversals. In particular, the model produces asymmetric reversals, a positive correlation of field strength and interval length, and a bimodal field distribution. All these features are attributable to the magnetic field dynamics in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. The negative slope of the growth rate curve between the nearby local maximum and the exceptional point makes the system unstable and drives it to the exceptional point and beyond into the oscillatory branch where the sign change happens. A weakness of this reversal model is the apparent necessity to fine-tune the magnetic Reynolds number and/or the radial profile of alpha. In the present paper, it is shown that this fine-tuning is not necessary in the case of higher supercriticality of the dynamo. Numerical examples and physical arguments are compiled to show that, with increasing magnetic Reynolds number, there is strong tendency for the exceptional point and the associated local maximum to move close to the zero growth rate line. Although exemplified again by the spherically symmetric alpha^2 dynamo model, the main idea of this ''self-tuning'' mechanism of saturated dynamos into a reversal-prone state seems well transferable to other dynamos. As a consequence, reversing dynamos might be much more typical and may occur much more frequently in nature than what could be expected from a purely kinematic perspective.

The analytical solutions for the general-four-wave-mixing hyperpolarizabilities $\chi^{(3)}(-(w_1+w_2+w_3);w_1,w_2,w_3)$ on infinite chains under both Su-Shrieffer-Heeger and Takayama-Lin-Liu-Maki models of trans-polyacetylene are obtained through the scheme of dipole-dipole correlation. Analytical expressions of DC Kerr effect $\chi^{(3)}(-w;0,0,w)$, DC-induced second harmonic generation $\chi^{(3)}(-2w;0,w,w)$, optical Kerr effect $\chi^{(3)}(-w;w,-w,w)$ and DC-electric-field-induced optical rectification $\chi^{(3)}(0;w,-w,0)$ are derived. By including or excluding ${\bf \nabla_k}$ terms in the calculations, comparisons show that the intraband contributions dominate the hyperpolarizabilities if they are included. $\nabla_k$ term or intraband transition leads to the break of the overall permutation symmetry in $\chi^{(3)}$ even for the low frequency and non-resonant regions. Hence it breaks the Kleinman symmetry that is directly based on the overall permutation symmetry. Our calculations provide a clear understanding of the Kleinman symmetry breaks that are widely observed in many experiments. We also suggest a feasible experiment on $\chi^{(3)}$ to test the validity of overall permutation symmetry and our theoretical prediction. Finally, our calculations show the following trends for the various third-order nonlinear optical processes in the low frequency and non-resonant region: $\chi^{(3)}(-3w;w,w,w)> \chi^{(3)}(-2w;0,w,w)> \chi^{(3)}(-w;w,-w,w)>\chi^{(3)}(-w; 0,0,w)>= \chi^{(3)}(0;w,-w,0)$, and in the resonant region: $\chi^{(3)}(-w;0,0,w)> \chi^{(3)}(-w;w,-w,w)> \chi^{(3)}(-2w;0,w,w)>\chi^{(3)}(0;w,-w,0)>\chi^{(3)}(-3w;w,w,w)$. (w=\omega)

With only a few exceptions, the numerical simulation of cosmic and laboratory hydromagnetic dynamos has been carried out in the framework of the differential equation method. However, the integral equation method is known to provide robust and accurate tools for the numerical solution of many problems in other fields of physics. The paper is intended to facilitate the use of integral equation solvers in dynamo theory. In concrete, the integral equation method is employed to solve the eigenvalue problem for a hydromagnetic dynamo model with a spherically symmetric, isotropic helical turbulence parameter alpha. Three examples of the function alpha(r) with steady and oscillatory solutions are considered. A convergence rate proportional to the inverse squared of the number of grid points is achieved. Based on this method, a convergence accelerating strategy is developed and the convergence rate is improved remarkably. Typically, quite accurate results can be obtained with a few tens of grid points. In order to demonstrate its suitability for the treatment of dynamos in other than spherical domains, the method is also applied to alpha^2 dynamos in rectangular boxes. The magnetic fields and the electric potentials for the first eigenvalues are visualized.

The homogeneous dynamo effect is at the root of cosmic magnetic field generation. With only a very few exceptions, the numerical treatment of homogeneous dynamos is carried out in the framework of the differential equation approach. The present paper tries to facilitate the use of integral equations in dynamo research. Apart from the pedagogical value to illustrate dynamo action within the well-known picture of Biot-Savart's law, the integral equation approach has a number of practical advantages. The first advantage is its proven numerical robustness and stability. The second and perhaps most important advantage is its applicability to dynamos in arbitrary geometries. The third advantage is its intimate connection to inverse problems relevant not only for dynamos but also for technical applications of magnetohydrodynamics. The paper provides the first general formulation and application of the integral equation approach to time-dependent kinematic dynamos in finite domains. For the spherically symmetric alpha^2 dynamo model it is shown how the general formulation is reduced to a coupled system of two radial integral equations for the defining scalars of the poloidal and the toroidal field components. The integral equation formulation for spherical dynamos with general velocity fields is also derived. Two numerical examples, the alpha^2 dynamo model with radially varying alpha, and the Bullard-Gellman model illustrate the equivalence of the approach with the usual differential equation method.