The aim was to assess epidemiological characteristics of the most recent consumption patterns of meat, vegetable, and fruit among representative urban and rural residents aged 60+ years in regional China. In this cross-sectional survey conducted in mid-2018, participants aged 60+ years were randomly chosen from urban and rural communities in Nanjing municipality of China. Meat, vegetable, and fruit intake were assessed with a validated food frequency questionnaire. Multivariate logistic regression models were applied to compute odds ratio (OR) and 95 % confidence interval (CI) to investigate the association of socio-demographic characteristics with a likelihood of meeting intake recommendation. Among the 20 867 participants, 49⋅5 % were men and 45⋅0 % urban elders, and 6⋅5 % aged 80+ years. The mean values of consumption frequency of red meat, white meat, vegetable, and fruit were 2⋅99 ± 2⋅28, 1⋅37 ± 1⋅13, 5⋅24 ± 6⋅43, and 2⋅64 ± 2⋅91 times/week, respectively, among overall participants. Moreover, there were 14⋅9, 23⋅7, and 12⋅1 % of participants meeting intake recommendations of meat, vegetable, and fruit, separately, in this study. After adjustment for potential confounders, age, gender, residence area, and educational attainment each was associated with the likelihood of meeting intake recommendation of meat, vegetable, or fruit. The consumption frequency and proportion of participants meeting intake recommendations of meat, vegetable, or fruit were not high among elders in regional China. Socio-demographic characteristics were associated with intake recommendations of meat, vegetables, and fruit. It has public health implications that participants’ socio-demographic attributes shall be considered for precision intervention on meat, vegetable, and fruit consumption in healthy eating campaigns among elders in China.

Dietary antioxidant indices (DAI) may be potentially associated with relative telomere length (RTL) of leucocytes. This study aimed to investigate the relationship between DAI and RTL. A cross-sectional study involving 1656 participants was conducted. A generalised linear regression model and a restricted cubic spline model were used to assess the correlation of DAI and its components with RTL. Generalised linear regression analysis revealed that DAI (β = 0·005, P = 0·002) and the intake of its constituents vitamin C (β = 0·043, P = 0·027), vitamin E (β = 0·088, P < 0·001), Se (β = 0·075, P = 0·003), and Zn (β = 0·075, P = 0·023) were significantly and positively correlated with RTL. Sex-stratified analysis showed that DAI (β = 0·006, P = 0·005) and its constituents vitamin E (β = 0·083, P = 0·012), Se (β = 0·093, P = 0·006), and Zn (β = 0·092, P = 0·034) were significantly and positively correlated with RTL among females. Meanwhile, among males, only vitamin E intake (β = 0·089, P = 0·013) was significantly and positively associated with RTL. Restricted cubic spline analysis revealed linear positive associations between DAI and its constituents’ (vitamin E, Se and Zn) intake and RTL in the total population. Sex-stratified analysis revealed a linear positive correlation between DAI and its constituents’ (vitamin E, Se and Zn) intake and RTL in females. Our study found a significant positive correlation between DAI and RTL, with sex differences.

This study investigates the effect of vibration on the flow structure transitions in thermal vibrational convection (TVC) systems, which occur when a fluid layer with a temperature gradient is excited by vibration. Direct numerical simulation (DNS) of TVC in a two-dimensional enclosed square box is performed over a range of dimensionless vibration amplitudes $0.001 \le a \le 0.3$ and angular frequencies $10^{2} \le \omega \le 10^{7}$, with a fixed Prandtl number of 4.38. The flow visualisation shows the transition behaviour of flow structure upon the varying frequency, characterising three distinct regimes, which are the periodic-circulation regime, columnar regime and columnar-broken regime. Different statistical properties are distinguished from the temperature and velocity fluctuations at the boundary layer and mid-height. Upon transition into the columnar regime, columnar thermal coherent structures are formed, in contrast to the periodic oscillating circulation. These columns are contributed by the merging of thermal plumes near the boundary layer, and the resultant thermal updrafts remain at almost fixed lateral position, leading to a decrease in fluctuations. We further find that the critical point of this transition can be described nicely by the vibrational Rayleigh number ${{Ra}}_{vib}$. As the frequency continues to increase, entering the so-called columnar-broken regime, the columnar structures are broken, and eventually the flow state becomes a large-scale circulation (LSC), characterised by a sudden increase in fluctuations. Finally, a phase diagram is constructed to summarise the flow structure transition over a wide range of vibration amplitude and frequency parameters.

A high-power all polarization-maintaining (PM) chirped pulse amplification (CPA) system operating in the 2.0 μm range is experimentally demonstrated. Large mode area (LMA) thulium-doped fiber (TDF) with a core/cladding diameter of 25/400 μm is employed to construct the main amplifier. Through dedicated coiling and cooling of the LMA-TDF to manage the loss of the higher order mode and thermal effect, a maximum average power of 314 W with a slope efficiency of 52% and polarization extinction ratio of 20 dB is realized. The pulse duration is compressed to 283 fs with a grating pair, corresponding to a calculated peak power of 10.8 MW, considering the compression efficiency of 88% and the estimated Strehl ratio of 89%. Moreover, through characterizing the noise properties of the laser, an integrated relative intensity noise of 0.11% at 100 Hz−1 MHz is obtained at the maximum output power, whereas the laser timing jitter is degraded by the final amplifier from 318 to 410 fs at an integration frequency of 5 kHz to 1 MHz, owing to the self-phase modulation effect-induced spectrum broadening. The root-mean-square of long-term power fluctuation is tested to be 0.6%, verifying the good stability of the laser operation. To the best of our knowledge, this is the highest average power of an ultrafast laser realized from an all-PM-fiber TDF-CPA system ever reported.

We numerically investigated the global linear instability and bifurcations in electro-thermo-convection (ETC) of a dielectric liquid confined in a two-dimensional (2-D) concentric annulus subjected to a strong unipolar injection. Seven kinds of solutions exist in this ETC system due to the complex bifurcations, i.e. saddle-node, subcritical and supercritical Hopf bifurcations. These bifurcation routes constitute at most four solution branches. Global linear instability analysis and energy analysis were conducted to explain the instability mechanism and transition of different solutions and to predict the local instability regions. The linearized lattice Boltzmann method (LLBM) for global linear instability analysis, first proposed by Pérez et al. (Theor. Comput. Fluid Dyn., vol. 31, 2017, pp. 643–664) to analyse incompressible flows, was extended here to solve the whole set of coupled linear equations, including the linear Navier–Stokes equations, the linear energy equation, Poisson's equation and the linear charge conservation equation. A multiscale analysis was also performed to recover the macroscopic linearized Navier–Stokes equations from the four different discrete lattice Boltzmann equations (LBEs). The LLBM was validated by calculating the linear critical value of 2-D natural convection; it has an error of 1.39% compared with the spectral method. Instability with global travelling wave behaviour is a unique behaviour in the annulus configuration electrothermohydrodynamic system, which may be caused by the baroclinity. Finally, the chaotic behaviour was quantitatively analysed through calculation of the fractal dimension and Lyapunov exponent.

We present the results of two population surveys conducted 10 years apart (December 2010–February 2011 and December 2020–January 2021) of the Critically Endangered white-headed langur Trachypithecus leucocephalus in the Chongzuo White-Headed Langur National Nature Reserve, Guangxi Province, China. In the first survey, we recorded 818 individuals in 105 groups and 16 solitary adult males. In the second survey, we recorded 1,183 individuals in 128 groups and one solitary adult male. As a result of government policies, poaching for food and traditional medicine is no longer a primary threat to these langurs. However, severe forest loss and fragmentation caused by human activities could limit any future increase of this langur population.

Straightplasma channels are widely used to guide relativistic intense laser pulses over several Rayleigh lengths for laser wakefield acceleration. Recently, a curved plasma channel with gradually varied curvature was suggested to guide a fresh intense laser pulse and merge it into a straight channel for staged wakefield acceleration [Phys. Rev. Lett. 120, 154801 (2018)]. In this work, we report the generation of such a curved plasma channel from a discharged capillary. Both longitudinal and transverse density distributions of the plasma inside the channel were diagnosed by analyzing the discharging spectroscopy. Effects of the gas-filling mode, back pressure and discharging voltage on the plasma density distribution inside the specially designed capillary are studied. Experiments show that a longitudinally uniform and transversely parabolic plasma channel with a maximum channel depth of 47.5 μm and length of 3 cm can be produced, which is temporally stable enough for laser guiding. Using such a plasma channel, a laser pulse with duration of 30 fs has been successfully guided along the channel with the propagation direction bent by 10.4°.

Since 2016, China has been conducting military flybys around Taiwan, while the US has approved arms sales to Taiwan on several occasions and sent warplanes and battleships through the Taiwan Strait. How does Taiwanese public opinion respond to the Chinese and US military presence in the Strait? Is the public likely to become less supportive of de jure independence for Taiwan on account of China's military deterrence or more supportive owing to a perceived likelihood of US military assistance? In this report, we provide answers to these questions based on evidence from a survey experiment conducted in Taiwan in October–November 2020. We find that Taiwanese are less sensitive to the Chinese military presence in the Taiwan Strait but have become more supportive of de jure independence after seeing the US aircraft in the area. Our findings contribute to studies of cross-Strait relations and US foreign policy on the Taiwan Strait.

An all-fiber high-power linearly polarized chirped pulse amplification (CPA) system is experimentally demonstrated. Through stretching the pulse duration to a full width of approximately 2 ns with two cascaded chirped fiber Bragg gratings (CFBGs), a maximum average output power of 612 W is achieved from a high-gain Yb-doped fiber that has a core diameter of 20 μm with a slope efficiency of approximately 68% at the repetition rate of 80 MHz. At the maximum output power, the polarization degree is 92.5% and the M2 factor of the output beam quality is approximately 1.29; the slight performance degradations are attributed to the thermal effects in the main amplifier. By optimizing the B-integral of the amplifier and finely adjusting the higher-order dispersion of one of the CFBGs, the pulse width is compressed to 863 fs at the highest power with a compression efficiency of 72%, corresponding to a maximum compressed average power of 440.6 W, single pulse energy of 5.5 μJ and peak power of about 4.67 MW. To the best of our knowledge, this is the highest average power of a femtosecond laser directly generated from an all-fiber linearly polarized CPA system.

In this paper, we prove that the ratio of the modulus of the iterates of two points in an escaping Fatou component could be bounded even if the orbit of the component contains a sequence of annuli whose moduli tend to infinity, and this cannot happen when the maximal modulus of the meromorphic function is uniformly large enough. In this way we extend certain related results for entire functions to meromorphic functions with infinitely many poles.

The new mineral, gysinite-(La), with the ideal formula PbLa(CO3)2(OH)⋅H2O, has been discovered in lujavrite from the Saima alkaline complex, Liaoning Province, China. It commonly occurs as subhedral to anhedral, granular and platy crystals of 5 to 50 μm in size, in interstices or enclosed in microcline, aegirine and nepheline. Associated minerals include nepheline, aegirine, microcline, natrolite, eudialyte, lamprophyllite, bastnäsite-(Ce), parasite-(Ce), ancylite-(La), ancylite-(Ce), bobtraillite, britholite-(Ce), thorite, calcite and galena. The crystallisation of gysinite-(La) may be related to the post-magmatic carbonation event. Gysinite-(La) crystals are generally transparent, colourless, or pale yellow, with a vitreous lustre and white streak. It is brittle with an uneven fracture, and the estimated Mohs hardness is 3½ to 4. The calculated density is 5.007 g/cm3. Optically, gysinite-(La) is biaxial (–), α= 1.832(2), β= 1.849(4), γ = 1.862(5) in white light and 2Vmeas = 81.6°. The empirical formula of gysinite-(La) is (La0.93Pb0.61Nd0.23Pr0.14Sr0.04Gd0.02Sm0.01Eu0.01Ca0.01)Σ2(CO3)2(OH)1.34⋅0.66H2O, which is calculated on the basis of general formula (REExM2+2–x)(CO3)2(OH)x⋅(2–x)H2O. The strongest eight lines of its powder X-ray diffraction pattern [d, Å (I, %) (hkl)] are: 5.596 (21) (011), 4.349 (100) (110), 3.732 (68) (111), 2.984 (61) (121), 2.667 (21) (031), 2.363 (48) (131), 2.090 (29) (221) and 2.028 (21) (212). Gysinite-(La) is orthorhombic, in the space group Pmcn, and unit-cell parameters refined from single-crystal X-ray diffraction data are: a = 5.0655(2) Å, b = 8.5990(3) Å, c = 7.3901(4) Å, V = 321.90(2) Å3 and Z = 2. It is a new member of the ancylite group and isostructural with gysinite-(Nd), but with La and Pb dominant in the metal cation sites in the structure.

We report that vertical vibration with small amplitude and high frequency can tame convective heat transport in Rayleigh–Bénard convection in a turbulent regime. When vertical vibration is applied, a dynamically averaged ‘anti-gravity’ results that stabilizes the thermal boundary layer and inhibits the eruption of thermal plumes. This eventually leads to the attenuation of the intensity of large-scale mean flow and a significant suppression of turbulent heat transport. Accounting for both the thermally led buoyancy and the vibration-induced anti-gravitational effects, we propose an effective Rayleigh number that helps to extend the Grossmann–Lohse theory to thermal vibrational turbulence. The prediction of the reduction on both the Nusselt and Reynolds numbers obtained by the extended model is found to agree well with the numerical data. In addition, vibrational influences on the mean flow structure and the temporal evolution of Nusselt and Reynolds numbers are investigated. The non-uniform characteristic of vibration-induced ‘anti-gravity’ is discussed. The present findings provide a powerful basis for studying thermal vibrational turbulence and put forward a novel strategy for actively controlling thermal turbulence.

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

We study the two-dimensional Rayleigh–Bénard instability subject to the combined effects of a solid–liquid phase boundary and shear using linear stability theory and energy analysis. We consider two thermal states of the solid (isothermal and conducting), and two types of shear that can arise in different contexts. When the melting temperature is equal to that held at the top boundary, three instability modes can arise with the increase of Reynolds number Re that characterizes the shear intensity: a boundary mode, a mixed boundary–bulk mode and a bulk flow mode. When the melting temperature lies between the top and the bottom boundaries, the introduction of Couette flow, independent of its intensity, always leads to the mixed mode, whereas the instability with Poiseuille flow is dominated by the bulk flow mode once Re exceeds a critical value, below which the mixed mode dominates. The energy analysis suggests that there exist two mechanisms by which the shear flow affects the system: one is by inhibiting the upward heat flux and another is by absorbing energy from the perturbed hydrodynamic field. These two mechanisms can play totally different roles in different cases. Results in the high-Re regime indicate that, when Re exceeds its classical threshold, i.e. Re = 5772.2 for Poiseuille flow, the Tollmien–Schlichting instability will be dominant in the present system.

Fluorsigaiite, ideally Ca2Sr3(PO4)3F, is a new Sr analogue of fluorphosphohedyphane and a new member of the apatite supergroup. It was discovered in lujavrite from the Saima alkaline complex, Liaoning Province, China. Fluorsigaiite commonly occurs as individual prismatic, columnar and platy crystals of 10 to 50 μm in size, associated with microcline, nepheline, aegirine, natrolite, eudialyte, fluorapatite, a fluorstrophite-like mineral, stronadelphite and calcite. Occasionally, crystals of fluorsigaiite form prismatic aggregates in the interstices of lujavrite. Fluorsigaiite is translucent to transparent, colourless to yellowish white with a vitreous lustre and without fluorescence. The estimated Mohs hardness is 5, and the tenacity is brittle with uneven fractures. The calculated density is 3.842 g/cm3. Optically, fluorsigaiite is uniaxial (–) with ω = 1.64(1) and ɛ = 1.63(1) in white light and without dispersion. The mean chemical composition (in wt.%) of fluorsigaiite is Na2O 0.75, CaO 15.17, SrO 44.44, La2O3 3.64, Ce2O3 2.22, Pr2O3 0.19, Nd2O3 0.13, Sm2O3 0.05, Gd2O3 0.23, P2O5 31.87, F 1.91, H2O 0.46, sum 100.26, giving the empirical formula (Sr2.82Ca1.79Na0.16La0.15Ce0.09Pr0.01Nd0.01Gd0.01)Σ5.04P2.97O12[F0.66(OH)0.34]Σ1, which is calculated on the basis of 13 total anions and F+(OH) = 1. The strongest eight lines of its powder X-ray diffraction pattern [d, Å (I, %) (hkl)] are: 3.563 (15) (002), 3.275 (15) (102), 3.144 (19) (120), 2.876 (100) (121), 2.861 (96) (112), 2.772 (27) (300), 1.991 (17) (222) and 1.895 (23) (213). Fluorsigaiite is hexagonal, in the space group P63/m and unit-cell parameters refined from single-crystal X-ray diffraction data are: a = 9.6101(2) Å, c = 7.1311(1) Å, V = 570.35(3) Å3 and Z = 2. It is isostructural with hedyphane-group minerals, and contains different prevailing (species-defining) Ca and Sr cations at the Ca1 and Ca2 sites, respectively. Fluorsigaiite was probably formed from Sr-rich fluids at late-magmatic or hydrothermal stage of the Saima lujavrite.

In this paper, we report that reversals of the large-scale circulation in two-dimensional Rayleigh–Bénard (RB) convection can be suppressed by imposing sinusoidally distributed heating to the bottom plate. We examine how the frequency of flow reversals depends on the dimensionless wavenumber $k$ of the spatial temperature modulation with various modulation amplitude $A$. For sufficiently large $k$, the flow reversal frequency is close to that in the standard RB convection under uniform heating. However, when $k$ decreases, the frequency of flow reversal gradually becomes lower and can even be largely reduced. Furthermore, we examine the growth of the corner roll and the global flow structure based on Fourier mode decomposition, and reveal that the size of the corner roll diminishes as the wavenumber decreases. The reason is that the regions occupied by the cold phase can absorb heat from the hot plumes and thus lower their temperature, which reduces the corner roll's kinetic energy input provided by the buoyancy force, and weakens the feeding process of the corner rolls. This results in the locking of the corner roll into a smaller region near the corner, making it harder for a reversal to occur. Using the concept of horizontal convection caused by non-uniform heating, we find a relevant parameter $k/A$ to describe briefly how the reversal frequency depends on wavenumber and modulation amplitude. The present work provides a new way to control the flow reversals in RB convection through modifying temperature boundary conditions.

Modal global linear stability analysis of thermal convection is performed with the linearized lattice Boltzmann method (LLBM). The onset of Rayleigh–Bénard convection in rectangular cavities with conducting and adiabatic sidewalls and the instability of two-dimensional (2-D) and three-dimensional (3-D) natural convection in cavities are studied. The method of linearizing the local equilibrium probability distribution function that was first proposed by Pérez et al. (Theor. Comp. Fluid Dyn., vol. 31, 2017, pp. 643–664) is extended to solve the coupled linear Navier–Stokes equations together with the linear energy equation in this work. A multiscale analysis is also performed to recover the macroscopic linear Navier–Stokes equations from the discrete lattice Boltzmann equations for both the single and multiple relaxation time models. The present LLBM is implemented in the framework of the Palabos library. It is validated by calculating the linear critical value of 2-D natural convection that the LLBM with the multiple relaxation time model has an error less than 1 % compared with the spectral method. The instability mechanism of the flow is explained by kinetic energy transfer analysis. It is shown that the buoyancy mechanism and inertial mechanism tend to stabilize the Hopf bifurcation of the 2-D natural convection at Pr < 0.08 and Pr > 1, respectively. For 3-D natural convection, subcritical bifurcation of the Hopf type is found for low-Prandtl-number fluids (Pr < 0.1).