The third-order law links energy transfer rates in the inertial range of magneto- hydrodynamic (MHD) turbulence with third-order structure functions. Anisotropy, a typical property in the solar wind, challenges the applicability of the third-order law with the isotropic assumption. To shed light on the energy transfer process in the presence of anisotropy, we conducted direct numerical simulations of forced MHD turbulence with normal and hyper-viscosity under various strengths of the external magnetic field ($B_0$), and calculated three forms of third-order structure function with or without averaging of the azimuthal or polar angles with respect to $B_0$ direction. Correspondingly, three estimated energy transfer rates were obtained. The result shows that the peak of normalized third-order structure function occurs at larger scales closer to the $B_0$ direction, and the maximum of longitudinal transfer rates shifts away from the $B_0$ direction at larger $B_0$. Compared with normal viscous cases, hyper-viscous cases can attain better separated inertial range, thus facilitating the estimation of the energy cascade rates. We find that the widespread use of the isotropic form of the third-order law in estimating the energy transfer rates is questionable in some cases, especially when the anisotropy arising from the mean magnetic field is inevitable. In contrast, the direction-averaged third-order structure function properly accounts for the effect of anisotropy and predicts the energy transfer rates and inertial range accurately, even at very high $B_0$. With limited statistics, the third-order structure function shows a stronger dependence on averaging of azimuthal angles than the time, especially for high $B_0$ cases. These findings provide insights into the anisotropic effect on the estimation of energy transfer rates.

We study the following rainbow version of subgraph containment problems in a family of (hyper)graphs, which generalizes the classical subgraph containment problems in a single host graph. For a collection $\mathit {\mathbf {G}}=\{G_1, G_2,\ldots , G_{m}\}$ of not necessarily distinct k-graphs on the same vertex set $[n]$, a (sub)graph H on $[n]$ is rainbow if there exists an injection $\varphi : E(H)\rightarrow [m]$, such that $e\in E(G_{\varphi (e)})$ for each $e\in E(H)$. Note that if $|E(H)|=m$, then $\varphi $ is a bijection, and thus H contains exactly one edge from each $G_i$.

Our main results focus on rainbow clique-factors in (hyper)graph systems with minimum d-degree conditions. Specifically, we establish the following:

1. (1) A rainbow analogue of an asymptotical version of the Hajnal–Szemerédi theorem, namely, if $t\mid n$ and $\delta (G_i)\geq (1-\frac {1}{t}+\varepsilon )n$ for each $i\in [\frac {n}{t}\binom {t}{2}]$, then $\mathit {\mathbf {G}}$ contains a rainbow $K_t$-factor;

2. (2) Essentially, a minimum d-degree condition forcing a perfect matching in a k-graph also forces rainbow perfect matchings in k-graph systems for $d\in [k-1]$.

The degree assumptions in both results are asymptotically best possible (although the minimum d-degree condition forcing a perfect matching in a k-graph is in general unknown). For (1), we also discuss two directed versions and a multipartite version. Finally, to establish these results, we in fact provide a general framework to attack this type of problem, which reduces it to subproblems with finitely many colors.

The plant Camellia fascicularis, belonging to family Theaceae, has high ornamental and medicinal value, and rare gene resources for genetic improvement of Camellia crops, but is currently threatened with extinction because of the unique and extremely small wild populations. Molecular markers have clarified the wild plant species’ genetic diversity structure, new gene resources and relationship with crops. This will be beneficial for conservation of these valuable crop-related wild species and crop improvement. In this study, we identified 95,979 microsatellite loci from 155,011 transcriptome unigenes, and developed 14 polymorphic expressed sequence tag-derived simple sequence repeat (EST-SSR) microsatellite markers for C. fascicularis. The number of alleles (Na) per locus was 2–8 with a mean of 4.86. The genetic diversity of 40 individuals from four natural populations of C. fascicularis was analysed using these polymorphic markers. The number of alleles (Na) for EST-SSR ranged from 2 to 5, with the expected heterozygosities (He) and observed heterozygosities (Ho) in all loci ranging from 0.183 to 0.683, and from 0.201 to 0.700, respectively, implying a rich genetic variation present in wild C. fascicularis populations. Moreover, the phylogenetic analysis among four populations, using the 14 EST-SSR markers developed in this study, grouped 40 individuals into three groups, which coincide with their geographic distribution. These results showed that 14 EST-SSR markers are available for the analysis of genetic variation in C. fascicularis populations and genetic improvement of new Camellias cultivars by interspecific hybridization, and are beneficial for conservation of the endangered species.

Shame leads to devaluation of the social self, and thus to a desire to improve self-esteem. Money, which is related to the notion of one’s ability, may help people demonstrate competence and gain self-esteem and respect from others. Based on the perspectives of feelings-as-information and threatened ego, we tested the hypothesis that a sense of shame heightens the desire for money, prompting self-interested behaviors as reflected by monetary donations and social value orientation. The results showed that subjects in the shame condition donated less money (Experiment 1) and exhibited more self-interested choices in the modified decomposed game (Experiment 2). The desire for money as reflected in overestimated coin sizes mediated the effect of shame on self-interested behavior. Our findings suggest that shame elicits the desire to acquire money to amend the threatened social self and improve self-esteem; however, it may induce a self-interested inclination that could harm social relationships.

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.

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.

Based on the measurements conducted over the landfast sea ice in Prydz Bay, East Antarctica during the sea-ice growth season in 2016, various parameterization schemes in the high-resolution thermodynamic snow/ice model HIGHTSI are evaluated. The parameterization scheme of turbulent fluxes produces the largest errors compared with the parameterization schemes for other surface heat fluxes. However, the sea-ice thickness simulation is most sensitive to the differences in upward longwave radiation at the surface. In addition, the sea-ice thickness simulation during the growth season is highly sensitive to the oceanic heat flux, and a new oceanic heat flux parameterization scheme based on the bulk method is proposed. The new parameterization scheme is tested in a second year, and it significantly improves the model performance relative to the standard configuration when compared against observations. Finally, the seasonal variation in the heat budget and its influence on the sea-ice thickness variation are analyzed. The net shortwave radiation, sensible heat flux and conductive heat flux (the net longwave radiation and latent heat flux) are found to be the surface heat sources (heat sinks) during the growth season. The larger conductive heat flux and the smaller oceanic heat flux can intensify the growth of sea ice.

As the corona virus disease 2019 (COVID-19) pandemic continues around the world, understanding the transmission characteristics of COVID-19 is vital for prevention and control. We conducted the first study aiming to estimate and compare the relative risk of secondary attack rates (SARs) of COVID-19 in different contact environments. Until 26 July 2021, epidemiological studies and cluster epidemic reports of COVID-19 were retrieved from SCI, Embase, PubMed, CNKI, Wanfang and CBM in English and Chinese, respectively. Relative risks (RRs) were estimated in pairwise comparisons of SARs between different contact environments using the frequentist NMA framework, and the ranking of risks in these environments was calculated using the surface under the cumulative ranking curve (SUCRA). Subgroup analysis was performed by regions. Thirty-two studies with 68 260 participants were identified. Compared with meal or gathering, transportation (RR 10.55, 95% confidence interval (CI) 1.43–77.85), medical care (RR 11.68, 95% CI 1.58–86.61) and work or study places (RR 10.15, 95% CI 1.40–73.38) had lower risk ratios for SARs. Overall, the SUCRA rankings from the highest to the lowest were household (95.3%), meal or gathering (81.4%), public places (58.9%), daily conversation (50.1%), transportation (30.8%), medical care (18.2%) and work or study places (15.3%). Household SARs were significantly higher than other environments in the subgroup of mainland China and sensitive analysis without small sample studies (<100). In light of the risks, stratified personal protection and public health measures need to be in place accordingly, so as close contacts categorising and management.

Arctic landfast sea ice (LFSI) represents an important quasi-stationary coastal zone. Its evolution is determined by the regional climate and bathymetry. This study investigated the seasonal cycle and interannual variations of LFSI along the northwest coast of Kotelny Island. Initial freezing, rapid ice formation, stable and decay stages were identified in the seasonal cycle based on application of the visual inspection approach (VIA) to MODIS/Envisat imagery and results from a thermodynamic snow/ice model. The modeled annual maximum ice thickness in 1995–2014 was 2.02 ± 0.12 m showing a trend of −0.13 m decade−1. Shortened ice season length (−22 d decade−1) from model results associated with substantial spring (2.3°C decade−1) and fall (1.9°C decade−1) warming. LFSI break-up resulted from combined fracturing and melting, and the local spatiotemporal patterns of break-up were associated with the irregular bathymetry. Melting dominated the LFSI break-up in the nearshore sheltered area, and the ice thickness decreased to an average of 0.50 m before the LFSI disappeared. For the LFSI adjacent to drift ice, fracturing was the dominant process and the average ice thickness was 1.56 m at the occurrence of the fracturing. The LFSI stages detected by VIA were supported by the model results.

A fever clinic within a hospital plays a vital role in pandemic control because it serves as an outpost for pandemic discovery, monitoring and handling. As the outbreak of coronavirus disease 2019 (COVID-19) in Wuhan was gradually brought under control, the fever clinic in the West Campus of Wuhan Union Hospital introduced a new model for construction and management of temporary mobile isolation wards. A traditional battlefield hospital model was combined with pandemic control regulations, to build a complex of mobile isolation wards that used adaptive design and construction for medical operational, medical waste management and water drainage systems. The mobile isolation wards allowed for the sharing of medical resources with the fever clinic. This increased the capacity and efficiency of receiving, screening, triaging and isolation and observation of patients with fever. The innovative mobile isolation wards also controlled new sudden outbreaks of COVID-19. We document the adaptive design and construction model of the novel complex of mobile isolation wards and explain its characteristics, functions and use.

We report on the generation of a mid-infrared (mid-IR) frequency comb with a maximum average output power of 250 mW and tunability in the 2.7–4.0 μm region. The approach is based on a single-stage difference frequency generation (DFG) starting from a compact Yb-doped fiber laser system. The repetition rate of the near-infrared (NIR) comb is locked at 75 MHz. The phase noise of the repetition rate in the offset-free mid-IR comb system is measured and analyzed. Except for the intrinsic of NIR comb, environmental noise at low frequency and quantum noise at high frequency from the amplifier chain and nonlinear spectral broadening are the main noise sources of broadening the linewidth of comb teeth, which limits the precision of mid-IR dual-comb spectroscopy.

In this paper, we propose and study an almost periodic reaction–diffusion epidemic model in which disease latency, spatial heterogeneity and general seasonal fluctuations are incorporated. The model is given by a spatially nonlocal reaction–diffusion system with a fixed time delay. We first characterise the upper Lyapunov exponent λ* for a class of almost periodic reaction–diffusion equations with a fixed time delay and provide a numerical method to compute it. On this basis, the global threshold dynamics of this model is established in terms of λ* It is shown that the disease-free almost periodic solution is globally attractive if λ* < 0, while the disease is persistent if λ* > 0. By virtue of numerical simulations, we investigate the effects of diffusion rate, incubation period and spatial heterogeneity on disease transmission.

We aimed to describe the clinical features in coronavirus disease 2019 (COVID-19) cases. We studied 134 critically ill COVID-19 cases from 30 December 2019 to 20 February 2020 in an intensive care unit (ICU) at Wuhan Jinyintan Hospital. Demographics, underlying diseases, therapy strategies and test results were collected and analysed from patients on admission, admission to the ICU and 48 h before death. The non-survivors were older (65.46 (s.d. 9.74) vs. 46.45 (s.d. 11.09)) and were more likely to have underlying diseases. The blood group distribution of the COVID-19 cases differed from that of the Han population in Wuhan, with type A being 43.85%; type B, 26.92%; type AB, 10% and type O, 19.23%. Non-survivors tend to develop more severe lymphopaenia, with higher C-reactive protein, interleukin-6, procalcitonin, D-dimer levels and gradually increased with time. The clinical manifestations were non-specific. Compared with survivors, non-survivors more likely to have organ function injury, and to receive mechanical ventilation, either invasively or noninvasively. Multiple organ failure and secondary bacterial infection in the later period is worthy of attention.

A novel ionic liquid/α-ZrP (C16MIM/α-ZrP) lamellar nanocomposite was fabricated via the electrostatic self-assembly deposition technique by using exfoliated α-ZrP nanosheets and guest molecules (1-hexadecyl-3-methylimidazolium bromide) as building blocks under mild conditions. C16MIM/α-ZrP nanocomposite was characterized by various analytical techniques such as X-ray diffraction (XRD), transmission electron microscopy (TEM), scanning electron microscope (SEM), Fourier transform infrared spectroscopy, and synchronous thermal analyzer. The net interlayer spacing of α-ZrP determined by XRD confirmed that the C16MIM cations formed a monolayer arrangement between the α-ZrP nanosheets. The morphology and microstructure of C16MIM/α-ZrP composite were observed using SEM and TEM. The C16MIM/α-ZrP modified glass carbon electrode exhibited excellent electrocatalytic activity toward the oxidation of nitrite in weak base media. The results obtained with differential pulse voltammetry demonstrated that the C16MIM/α-ZrP hybrid detected nitrite linearly in the concentration range from 7.3 μM to 1.25 mM with the detection limit of 1.26 μM (S/N = 3). Additionally, the prepared sensor showed outstanding reproducibility, high stability, and anti-interference capability.

Snow and ice were monitored by thermistor-string-based Snow and Ice Mass Balance Array (SIMBA) in Lake Orajärvi in northern Finland. An existing automatic SIMBA-algorithm was further developed to derive air/snow, snow/ice and ice/water interfaces based on the SIMBA environment temperature (ET) profiles. The identified interfaces agreed with in situ observations made in 2011/12 winter season. The method was capable to identify upward-moving snow/ice interface that was also visible from SIMBA heating temperature (HT) profiles, which responds to differences in the thermal diffusivities of air, snow, ice and water. The SIMBA data obtained in winters 2017/18 and 2018/19 were used to investigate snow and ice mass balance. An upward-moving snow/ice interface was detected as a result of meteoric ice (snow ice and superimposed ice) formation. Snow contributed to granular lake ice formation up to 40–55% of the total ice thickness on the seasonal mean. Heavy snowfalls and low air temperature in early winter are favourable for granular ice formation. The seasonal mean snow depth on nearby land was 2.7–2.9 times of that on the lake. The estimation of freeboard from snow and ice mass-balance measurement is sensitive to the snow density. Accurate ice freeboard calculation is still a challenge.

A Fast Ice Prediction System (FIPS) was constructed and is the first regional land-fast sea-ice forecasting system for the Antarctic. FIPS had two components: (1) near-real-time information on the ice-covered area from MODIS and SAR imagery that revealed, tidal cracks, ridged and rafted ice regions; (2) a high-resolution 1-D thermodynamic snow and ice model (HIGHTSI) that was extended to perform a 2-D simulation on snow and ice evolution using atmospheric forcing from ECMWF: either using ERA-Interim reanalysis (in hindcast mode) or HERS operational 10-day predictions (in forecast mode). A hindcast experiment for the 2015 season was in good agreement with field observations, with a mean bias of 0.14 ± 0.07 m and a correlation coefficient of 0.98 for modeled ice thickness. The errors are largely caused by a cold bias in the atmospheric forcing. The thick snow cover during the 2015 season led to modeled formation of extensive snow ice and superimposed ice. The first FIPS operational service was performed during the 2017/18 season. The system predicted a realistic ice thickness and onset of snow surface melt as well as the area of internal ice melt. The model results on the snow and ice properties were considered by the captain of R/V Xuelong when optimizing a low-risk route for on-ice transportation through fast ice to the coastal Zhongshan Station.

We examine the relative effect of warming events (storms) and snow cover on thermodynamic growth of Arctic sea ice in winter. We use a 1-D snow and ice thermodynamic model to perform sensitivity experiments. Observations from the winter period of the Norwegian young sea ICE (N-ICE2015) campaign north of Svalbard are used to initiate and force the model. The N-ICE2015 winter was characterized by frequent storm events that brought pulses of heat and moisture, and a thick snow cover atop the sea ice (0.3–0.5 m). By the end of the winter, sea-ice bottom growth was negligible. We show that the thermodynamic effect of storms to the winter sea-ice growth is controlled by the amount of snow on sea ice. For 1.3 m initial ice thickness, the decrease in ice growth caused by the warming events ranged from −1.4% (for 0.5 m of snow) to −7.5% (for snow-free conditions). The decrease in sea-ice growth caused by the thick snow (0.5 m) was more important, ranging from −17% (with storms) to −23% (without storms). The results showcase the critical role of snow on winter Arctic sea-ice growth.

We describe an ultra-wide-bandwidth, low-frequency receiver recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band ( ${\sim}60\%$ ), the system temperature is approximately 22 K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver, including its astronomical objectives, as well as the feed, receiver, digitiser, and signal processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified, including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration, and timing stability.