Contrafreeloading (CFL) refers to animals’ tendency to prefer obtaining food through effort rather than accessing food that is freely available. Researchers have proposed various hypotheses to explain this intriguing phenomenon, but few studies have provided a comprehensive analysis of the factors influencing this behaviour. In this study, we observed the choice of alternative food containers in budgerigars (Melopsittacus undulatus) to investigate their CFL tendencies and the effects of pre-training, food deprivation, and effort required on the CFL tasks. The results showed that budgerigars did not exhibit significant difference in their first choices or the time interacting with less challenging versus more challenging food containers. Moreover, when evaluating each budgerigar’s CFL level, only half of them were identified as strong contrafreeloaders. Thus, we suggest that budgerigars exhibit an intermediate CFL level that lies somewhere between a strong tendency and the absence of such behaviour. Furthermore, we also found that food-deprived budgerigars tended to select less challenging food containers, and pre-trained budgerigars were more likely to choose highly challenging food containers than moderately challenging food containers, which means that the requirement of only a reasonable effort (access to food from moderately challenging food containers in this study) and the experience of pre-training act to enhance their CFL levels, whereas the requirement of greater effort and the experience of food deprivation act to decrease their CFL levels. Studying animal CFL can help understand why animals choose to expend effort to obtain food rather than accessing it for free, and it also has implications for setting feeding environments to enhance the animal welfare of captive and domesticated animals.

Recent experiments and simulations have sparked growing interest in the study of Rayleigh–Bénard convection in very slender cells. One pivotal inquiry arising from this interest is the elucidation of the flow structure within these very slender cells. Here we employ tomographic particle image velocimetry, for the first time, to capture experimentally the full-field three-dimensional and three-component velocity field in a very slender cylindrical cell with aspect ratio $\Gamma =1/10$. The experiments cover a Rayleigh number range $5.0 \times 10^8 \leqslant Ra \leqslant 5.0 \times 10^9$ and Prandtl number 5.7. Our experiments reveal that the flow structure in the $\Gamma =1/10$ cell is neither in the multiple-roll form nor in the simple helical form; instead, the ascending and descending flows can intersect and cross each other, resulting in the crossing events. These crossing events separate the flow into segments; within each segment, the ascending and descending flows ascend or descend side by side vertically or in the twisting manner, and the twisting is not unidirectional, while the segments near the boundary can also be in the form of a donut like structure. By applying the mode decomposition analyses to the measured three-dimensional velocity fields, we identified the crossing events as well as the twisting events for each instantaneous flow field. Statistical analysis of the modes reveals that as $Ra$ increases, the average length of the segments becomes smaller, and the average number of segments increases from 2.5 to 3.9 in the $Ra$ range of our experiments.

Although it is well established that gestational diabetes mellitus (GDM) is associated with fetal overgrowth in singleton pregnancies, little is known about its role in twins. We aimed to explore the relationship between GDM and the longitudinal fetal growth in twin pregnancies. This was a retrospective matched cohort study of GDM and non-GDM twin pregnancies delivered ≥36 weeks without other complications. All the women performed ≥3 ultrasounds after 22 weeks. Linear mixed models (LMMs) were used to explore the relationships between longitudinal fetal growth trajectories and GDM. Group-based trajectory modeling (GBTM) and generalized estimating equation (GEE) were applied to identify the latent growth patterns and investigate their relationships with GDM. In total, 215 GDM and 645 non-GDM twins were included, the majority of the patients did not require medication therapy (n = 202, GDMA1). LMM revealed that, compared with non-GDM, GDM was associated with an average increase in fetal weight of 4.36 g (95% CI [1.25, 7.48]) per week. GBTM and GEE further revealed that GDM increased the odds of fetal weight trajectory to nearly 40% of the total fetal weight trajectory, classified into the high-speed group (aOR = 1.39, 95% CI [1.03, 1.88]), associating with a 49.44 g (95% CI [11.41, 87.48]) increase in birth weight. Subgroup analysis revealed that all these differences were only significant among the GDMA1 pregnancies (p < .05). GDM (GDMA1) is significantly associated with an increase in fetal weight during gestation in twin pregnancies. However, this acceleration is mild, and its significance requires further exploration.

Electron cyclotron resonance ion thrusters (ECRITs) have the potential to be used for space gravitational wave detection due to their wide thrust range. However, an unclear understanding of dynamic processes of ECRITs with strongly coupled multi-operating parameters limits further improvements on thrust noise and response velocity by feedback control systems. An integrative mathematical model considering the non-Maxwell electron energy distribution function for ECRITs is validated by experiments and used to study the influence of operating parameters on the dynamic processes of thrusters, which provides a new simplified grid model. Simulation results show the response processes with microwave (MW) power can be divided into two stages. The characteristic times of the first and second stages are respectively several microseconds and 10 ms, which are respectively dominated by plasma motion and the volume effect. The overshoot of screen grid (SG) current decreases, while its response time remains unchanged when the response time of MW power is prolonged. The response time of SG current with a step increase of flow rate is approximately 10 ms, consistent with the volume effect. The SG current decreases with rise of flow rate for high flow rate operations due to the small increment of ion density limited by low electron temperature, the decrease of ion Bohm velocity and reduction of sheath extraction area. The influence of grid voltage on the dynamic process of the SG current depends on variation ranges of extraction capabilities. When variations of sheath extraction area are limited, the response time is 5 μs, consistent with plasma response time. It is prolonged to 0.5 ms if sheath extraction area variations are large because they cause obvious variations of plasma parameters in the discharge chamber. These dynamic results can not only facilitate designing feedback controllers of micro-propulsion systems for high-precision space missions, but also provide guidance for ion sources to generate highly stable or rapid-response ion beam.

Rhopalosiphum padi is an important grain pest, causing severe losses during crop production. As a systemic insecticide, flonicamid can control piercing-sucking pests efficiently. In our study, the lethal effects of flonicamid on the biological traits of R. padi were investigated via a life table approach. Flonicamid is highly efficiently toxic to R. padi, with an LC50 of 9.068 mg L−1. The adult longevity and fecundity of the R. padi F0 generation were markedly reduced under the LC25 and LC50 concentrations of flonicamid exposure. In addition, negative transgenerational effects on R. padi were observed under exposure to lethal concentrations of flonicamid, with noticeable decreases in the reproductive period, adult longevity, total longevity, and total fecundity of the F1 generation under the LC25 concentration of flonicamid. Furthermore, the third nymph stage (N3), preadult stage, duration of the adult pre-reproductive period, duration of the total pre-reproductive period, reproductive period, adult longevity, total longevity, and total fecundity of the F1 generation were significantly lower under treatment with the LC50 concentration of flonicamid. The life table parameters were subsequently analysed, revealing that the intrinsic rate of increase (rm) and the net reproductive rate (R0) were significantly lower but that the finite rate of increase (λ) and the mean generation time (T) were not significantly different under the LC25 and LC50 concentrations of flonicamid. These data are beneficial for grain aphid control and are critical for exploring the role of flonicamid in the integrated management of this key pest.

Observations of the intracluster medium (ICM) in the outskirts of galaxy clusters reveal shocks associated with gas accretion from the cosmic web. Previous work based on non-radiative cosmological hydrodynamical simulations have defined the shock radius, $r_{\text{shock}}$, using the ICM entropy, $K \propto T/{n_\mathrm{e}}^{2/3}$, where T and $n_{\text{e}}$ are the ICM temperature and electron density, respectively; the $r_{\text{shock}}$ is identified with either the radius at which K is a maximum or at which its logarithmic slope is a minimum. We investigate the relationship between $r_{\text{shock}}$, which is driven by gravitational hydrodynamics and shocks, and the splashback radius, $r_{\text{splash}}$, which is driven by the gravitational dynamics of cluster stars and dark matter and is measured from their mass profile. Using 324 clusters from The Three Hundred project of cosmological galaxy formation simulations, we quantify statistically how $r_{\text{shock}}$ relates to $r_{\text{splash}}$. Depending on our definition, we find that the median $r_{\text{shock}} \simeq 1.38 r_{\text{splash}} (2.58 R_{200})$ when K reaches its maximum and $r_{\text{shock}} \simeq 1.91 r_{\text{splash}} (3.54 R_{200})$ when its logarithmic slope is a minimum; the best-fit linear relation increases as $r_{\text{shock}} \propto 0.65 r_{\text{splash}}$. We find that $r_{\text{shock}}/R_{200}$ and $r_{\text{splash}}/R_{200}$ anti-correlate with virial mass, $M_{200}$, and recent mass accretion history, and $r_{\text{shock}}/r_{\text{splash}}$ tends to be larger for clusters with higher recent accretion rates. We discuss prospects for measuring $r_{\text{shock}}$ observationally and how the relationship between $r_{\text{shock}}$ and $r_{\text{splash}}$ can be used to improve constraints from radio, X-ray, and thermal Sunyaev-Zeldovich surveys that target the interface between the cosmic web and clusters.

Diagnostic classification models (DCMs) have seen wide applications in educational and psychological measurement, especially in formative assessment. DCMs in the presence of testlets have been studied in recent literature. A key ingredient in the statistical modeling and analysis of testlet-based DCMs is the superposition of two latent structures, the attribute profile and the testlet effect. This paper extends the standard testlet DINA (T-DINA) model to accommodate the potential correlation between the two latent structures. Model identifiability is studied and a set of sufficient conditions are proposed. As a byproduct, the identifiability of the standard T-DINA is also established. The proposed model is applied to a dataset from the 2015 Programme for International Student Assessment. Comparisons are made with DINA and T-DINA, showing that there is substantial improvement in terms of the goodness of fit. Simulations are conducted to assess the performance of the new method under various settings.

Time limits are imposed on many computer-based assessments, and it is common to observe examinees who run out of time, resulting in missingness due to not-reached items. The present study proposes an approach to account for the missing mechanisms of not-reached items via response time censoring. The censoring mechanism is directly incorporated into the observed likelihood of item responses and response times. A marginal maximum likelihood estimator is proposed, and its asymptotic properties are established. The proposed method was evaluated and compared to several alternative approaches that ignore the censoring through simulation studies. An empirical study based on the PISA 2018 Science Test was further conducted.

Developing large-eddy simulation (LES) wall models for separated flows is challenging. We propose to leverage the significance of separated flow data, for which existing theories are not applicable, and the existing knowledge of wall-bounded flows (such as the law of the wall) along with embedded learning to address this issue. The proposed so-called features-embedded-learning (FEL) wall model comprises two submodels: one for predicting the wall shear stress and another for calculating the eddy viscosity at the first off-wall grid nodes. We train the former using the wall-resolved LES (WRLES) data of the periodic hill flow and the law of the wall. For the latter, we propose a modified mixing length model, with the model coefficient trained using the ensemble Kalman method. The proposed FEL model is assessed using the separated flows with different flow configurations, grid resolutions and Reynolds numbers. Overall good a posteriori performance is observed for predicting the statistics of the recirculation bubble, wall stresses and turbulence characteristics. The statistics of the modelled subgrid-scale (SGS) stresses at the first off-wall grids are compared with those calculated using the WRLES data. The comparison shows that the amplitude and distribution of the SGS stresses and energy transfer obtained using the proposed model agree better with the reference data when compared with the conventional SGS model.

In this work, the shape of a bluff body is optimized to mitigate velocity fluctuations of turbulent wake flows based on large-eddy simulations (LES). The Reynolds-averaged Navier–Stokes method fails to capture velocity fluctuations, while direct numerical simulations are computationally prohibitive. This necessitates using the LES method for shape optimization given its scale-resolving capability and relatively affordable computational cost. However, using LES for optimization faces challenges in sensitivity estimation as the chaotic nature of turbulent flows can lead to the blowup of the conventional adjoint-based gradient. Here, we propose using the regularized ensemble Kalman method for the LES-based optimization. The method is a statistical optimization approach that uses the sample covariance between geometric parameters and LES predictions to estimate the model gradient, circumventing the blowup issue of the adjoint method for chaotic systems. Moreover, the method allows for the imposition of smoothness constraints with one additional regularization step. The ensemble-based gradient is first evaluated for the Lorenz system, demonstrating its accuracy in the gradient calculation of the chaotic problem. Further, with the proposed method, the cylinder is optimized to be an asymmetric oval, which significantly reduces turbulent kinetic energy and meander amplitudes in the wake flows. The spectral analysis methods are used to characterize the flow field around the optimized shape, identifying large-scale flow structures responsible for the reduction in velocity fluctuations. Furthermore, it is found that the velocity difference in the shear layer is decreased with the shape change, which alleviates the Kelvin–Helmholtz instability and the wake meandering.

We demonstrated a method to improve the output performance of a Ti:sapphire laser in the long-wavelength low-gain region with an efficient stimulated Raman scattering process. By shifting the wavelength of the high-gain-band Ti:sapphire laser to the long-wavelength low-gain region, high-performance Stokes operation was achieved in the original long-wavelength low-gain region of the Ti:sapphire laser. With the fundamental wavelength tuning from 870 to 930 nm, first-order Stokes output exceeding 2.5 W was obtained at 930–1000 nm, which was significantly higher than that directly generated by the Ti:sapphire laser, accompanied by better beam quality, shorter pulse duration and narrower linewidth. Under the pump power of 42.1 W, a maximum first-order Stokes power of 3.24 W was obtained at 960 nm, with a conversion efficiency of 7.7%. Furthermore, self-mode-locked modulations of first- and second-order Stokes generation were observed in Ti:sapphire intracavity solid Raman lasers for the first time.

The incidence of obesity-related glomerulopathy (ORG) is rising worldwide with very limited treatment methods. Paralleled with the gut–kidney axis theory, the beneficial effects of butyrate, one of the short-chain fatty acids (SCFA) produced by gut microbiota, on metabolism and certain kidney diseases have gained growing attention. However, the effects of butyrate on ORG and its underlying mechanism are largely unexplored. In this study, a mice model of ORG was established with a high-fat diet feeding for 16 weeks, and sodium butyrate treatment was initiated at the 8th week. Podocyte injury, oxidative stress and mitochondria function were evaluated in mice kidney and validated in vitro in palmitic acid-treated-mouse podocyte cell lines. Further, the molecular mechanisms of butyrate on podocytes were explored. Compared with controls, sodium butyrate treatment alleviated kidney injuries and renal oxidative stress in high-fat diet-fed mice. In mouse podocyte cell lines, butyrate ameliorated palmitic acid-induced podocyte damage and helped maintain the structure and function of the mitochondria. Moreover, the effects of butyrate on podocytes were mediated via the GPR43-Sirt3 signal pathway, as evidenced by the diminished effects of butyrate with the intervention of GPR43 or Sirt3 inhibitors. In summary, we conclude that butyrate has therapeutic potential for the treatment of ORG. It attenuates high-fat diet-induced ORG and podocyte injuries through the activation of the GPR43-Sirt3 signalling pathway.

This study examined the sour grapes/sweet lemons rationalization through 2 conditions: ‘attainable’ (sweet lemons) and ‘unattainable’ (sour grapes), reflecting China’s 2019-nCoV vaccination strategy. The aim was to find ways to change people’s beliefs and preferences regarding vaccines by easing their safety concerns and encouraging more willingness to get vaccinated. An online survey was conducted from January 22 to 27, 2021, with 3,123 residents across 30 provinces and municipalities in the Chinese mainland. The direction of belief and preference changed in line with the sour grapes/sweet lemons rationalization. Using hypothetical and real contrasts, we compared those for whom the vaccine was relatively unattainable (‘sour grapes’ condition) with those who could get the vaccine easily (‘sweet lemons’). Whether the vaccine was attainable was determined in the early stage of the vaccine roll-out by membership in a select group of workers that was supposed to be vaccinated to the greatest extent possible, or, by being in the second stage when the vaccine was available to all. The attainable conditions demonstrated higher evaluation in vaccine safety, higher willingness to be vaccinated, and lower willingness to wait and see. Hence, we propose that the manipulation of vaccine attainability, which formed the basis of the application of sour grapes/sweet lemons rationalization, can be utilized as a means to manipulate the choice architecture to nudge individuals to ease vaccine safety concerns, reducing wait-and-see tendencies, and enhancing vaccination willingness. This approach can expedite universal vaccination and its associated benefits in future scenarios resembling the 2019-nCoV vaccine rollout.

The AIMTB rapid test assay is an emerging test, which adopted a fluorescence immunochromatographic assay to measure interferon-γ (IFN-γ) production following stimulation of effector memory T cells in whole blood by mycobacterial proteins. The aim of this article was to explore the ability of AIMTB rapid test assay in detecting Mycobacterium tuberculosis (MTB) infection compared with the widely applied QuantiFERON-TB Gold Plus (QFT-Plus) test among rural doctors in China. In total, 511 participants were included in the survey. The concordance between the QFT-Plus test and the AIMTB rapid test assay was 94.47% with a Cohen’s kappa coefficient (κ) of 0.84 (95% CI, 0.79–0.90). Improved concordance between the two tests was observed in males and in participants with 26 or more years of service as rural doctors. The quantitative values of the QFT-Plus test was higher in individuals with a result of QFT-Plus-/AIMTB+ as compared to those with a result of QFT-Plus-/AIMTB- (p < 0.001). Overall, our study found that there was an excellent consistency between the AIMTB rapid test assay and the QFT-Plus test in a Chinese population. As the AIMTB rapid test assay is fast and easy to operate, it has the potential to improve latent tuberculosis infection testing and treatment at the community level in resource-limited settings.

Depression is highly prevalent in haemodialysis patients, and diet might play an important role. Therefore, we conducted this cross-sectional study to determine the association between dietary fatty acids (FA) consumption and the prevalence of depression in maintenance haemodialysis (MHD) patients. Dietary intake was assessed using a validated FFQ between December 2021 and January 2022. The daily intake of dietary FA was categorised into three groups, and the lowest tertile was used as the reference category. Depression was assessed using the Patient Health Questionnaire-9. Logistic regression and restricted cubic spline (RCS) models were applied to assess the relationship between dietary FA intake and the prevalence of depression. As a result, after adjustment for potential confounders, a higher intake of total FA [odds ratio (OR)T3 vs. T1 = 1·59, 95 % confidence interval (CI) = 1·04, 2·46] and saturated fatty acids (SFA) (ORT3 vs. T1 = 1·83, 95 % CI = 1·19, 2·84) was associated with a higher prevalence of depressive symptoms. Significant positive linear trends were also observed (P < 0·05) except for SFA intake. Similarly, the prevalence of depression in MHD patients increased by 20% (OR = 1.20, 95% CI = 1.01–1.43) for each standard deviation increment in SFA intake. RCS analysis indicated an inverse U-shaped correlation between SFA and depression (Pnonlinear > 0·05). Additionally, the sensitivity analysis produced similar results. Furthermore, no statistically significant association was observed in the subgroup analysis with significant interaction. In conclusion, higher total dietary FA and SFA were positively associated with depressive symptoms among MHD patients. These findings inform future research exploring potential mechanism underlying the association between dietary FA and depressive symptoms in MHD patients.

We determine the order of the k-core in a large class of dense graph sequences. Let $G_n$ be a sequence of undirected, n-vertex graphs with edge weights $\{a^n_{i,j}\}_{i,j \in [n]}$ that converges to a graphon $W\colon[0,1]^2 \to [0,+\infty)$ in the cut metric. Keeping an edge (i,j) of $G_n$ with probability ${a^n_{i,j}}/{n}$ independently, we obtain a sequence of random graphs $G_n({1}/{n})$. Using a branching process and the theory of dense graph limits, under mild assumptions we obtain the order of the k-core of random graphs $G_n({1}/{n})$. Our result can also be used to obtain the threshold of appearance of a k-core of order n.