The present study aims at studying the flow and field produced by a stellar radiative zone which is initially made to rotate differentially in the presence of a large-scale poloidal magnetic field threading the whole domain. We focus both on the axisymmetric configurations produced by the initial winding-up of the magnetic field lines and on the possible instabilities of those configurations. We aim in particular at assessing the role of the stratification at stabilising the system. We perform 2D and 3D global Boussinesq numerical simulations started from an initial radial or cylindrical differential rotation and a large-scale poloidal magnetic field. Under the conditions of a large rotation frequency compared to the Alfv\'en frequency, a magnetic configuration strongly dominated by its toroidal component is built. The parameters of the simulations are chosen to respect the ordering of time scales of a typical stellar radiative zone. In this framework, the axisymmetric evolution is studied by varying the relative effects of the thermal diffusion, the Br\"unt-V\"ais\"al\"a frequency, the rotation and the initial poloidal field strength. We find that the axisymmetric state only depends on $t_{es}/t_{Ap}$, the ratio between the Eddington-Sweet circulation time scale and the Alfv\'en time scale. A scale analysis of the Boussinesq equations allows us to recover this result. In the cylindrical case, a magneto-rotational instability develops when the thermal diffusivity is sufficiently high to enable the favored wavenumbers to be insensitive to the effects of the stable stratification. In the radial case, the magneto-rotational instability is driven by the latitudinal shear created by the back-reaction of the Lorentz force on the flow. Increasing the level of stratification then leaves the growth rate of the instability mainly unaffected while its horizontal length scale grows.

An extremely weak circularly polarized signature was recently discovered in spectral lines of the chemically peculiar Am star Sirius A. A weak surface magnetic field was proposed to account for the observed polarized signal, but the shape of the phase-averaged signature, dominated by a prominent positive lobe, is not expected in the standard theory of the Zeeman effect. We aim at verifying the presence of weak circularly polarized signatures in two other bright Am stars, beta UMa and theta Leo, and investigating the physical origin of Sirius-like polarized signals further. We present here a set of deep spectropolarimetric observations of beta UMa and theta Leo, observed with the NARVAL spectropolarimeter. We analyzed all spectra with the Least Squares Deconvolution multiline procedure. To improve the signal-to-noise ratio and detect extremely weak signatures in Stokes V profiles, we co-added all available spectra of each star (around 150 observations each time). Finally, we ran several tests to evaluate whether the detected signatures are consistent with the behavior expected from the Zeeman effect. The line profiles of the two stars display circularly polarized signatures similar in shape and amplitude to the observations previously gathered for Sirius A. Our series of tests brings further evidence of a magnetic origin of the recorded signal. These new detections suggest that very weak magnetic fields may well be present in the photospheres of a significant fraction of intermediate-mass stars. The strongly asymmetric Zeeman signatures measured so far in Am stars (featuring a dominant single-sign lobe) are not expected in the standard theory of the Zeeman effect and may be linked to sharp vertical gradients in photospheric velocities and magnetic field strengths.

Recent spectropolarimetric surveys of main-sequence intermediate-mass stars have exhibited a dichotomy in the distribution of the observed magnetic field between the kG dipoles of Ap/Bp stars and the sub-Gauss magnetism of Vega and Sirius. We would like to test whether this dichotomy is linked to the stability versus instability of large-scale magnetic configurations in differentially rotating radiative zones. We computed the axisymmetric magnetic field obtained from the evolution of a dipolar field threading a differentially rotating shell. A full parameter study including various density profiles and initial and boundary conditions was performed with a 2D numerical code. We then focused on the ratio between the toroidal and poloidal components of the magnetic field and discuss the stability of the configurations dominated by the toroidal component using local stability criteria and insights from recent 3D numerical simulations. The numerical results and a simple model show that the ratio between the toroidal and the poloidal magnetic fields is highest after an Alfv\'en crossing time of the initial poloidal field. For high density contrasts, this ratio converges towards an asymptotic value that can thus be extrapolated to realistic stellar cases. We then consider the stability of the magnetic configurations to non-axisymmetric perturbations and find that configurations dominated by the toroidal component are likely to be unstable if the shear strength is significantly higher than the poloidal Alfv\'en frequency. An expression for the critical poloidal field below which magnetic fields are likely to be unstable is found and is compared to the lower bound of Ap/Bp magnetic fields.