Chest X-ray (CXR) is the first tool globally employed to detect cardiopulmonary pathologies. These acquisitions are highly affected by scattered photons due to the large field-of-view required. Scatter in CXRs introduces background in the images, which reduces their contrast. We developed three deep learning-based models to estimate and correct the scatter contribution to CXRs. We used a Monte Carlo (MC) ray-tracing model to simulate CXRs from human models obtained from CT scans using different configurations (depending on the availability of dual-energy acquisitions). The simulated CXRs contained the separated contribution of direct and scattered X-rays in the detector. These simulated datasets were then used as the reference for the supervised training of several NN. Three NN models (single and dual energy) were trained with the MultiResUNet architecture. The performance of the NN models was evaluated on CXRs obtained with the MC code from chest CTs of patients affected by COVID-19. The results showed that the NN models were able to estimate and correct the scatter contribution to CXRs with an error <5%, being robust to variations in the simulation setup and improving the contrast in soft tissue. The single-energy model was tested with real CXRs, providing robust estimations of the scatter-corrected CXRs.

(1) Background: The imaging energy range of a typical Compton camera is limited due to the fact that scattered gamma photons are seldom fully absorbed when the incident energies are above 3 MeV. Further improving the upper energy limit of gamma-ray imaging has important application significance in active interrogation of special nuclear materials and chemical warfare agents, as well as range verification of proton therapy; (2) Methods: To realize gamma-ray imag-ing in a wide energy range of 0.3~7 MeV, a principle prototype, named a portable three-layer Compton camera, was developed using the scintillation detector that consist of an silicon photo-multiplier array coupled with a Gd3Al2Ga3O12:Ce pixelated scintillator array. Implemented in a list-mode maximum likelihood expectation maximization algorithm, a far-field energy-domain imaging method based on two-interaction and three-interaction events was applied to estimate the initial energy and spatial distribution of gamma-ray sources. The simulation model of detec-tors was established based on Monte Carlo simulation toolkit Geant4. The reconstructed images of a 133Ba, a 137Cs and a 60Co point-like sources were successfully obtained with our prototype in laboratory tests and compared with simulation studies; (3) Results: The proportion of effective imaging events accounted for about 2%, which made our prototype realize reconstructing the distribution of a 0.05 μSv/h 137Cs source in 10 seconds. The angular resolution for resolving two 137Cs point-like sources was 15°. Additional simulated imaging of the 6.13 MeV gamma-rays from 14.1 MeV neutron scattering with water preliminarily demonstrated the imaging capability for high incident energy; (4) Conclusions: We concluded that the prototype had good imaging per-formance in a wide energy range (0.3~7 MeV), which showed potential in several MeV gam-ma-ray imaging applications.

With data for SARS-CoV-2 and with many countries entering the second wave of contagion required the improvement of the forecasting model, structuring its model to forecast the peak of the first and second contagion wave in Brazil. The Model Estimation for SARS-CoV-2 Peak Contagion (MESPC) was structured, capable of estimating the peak of contagion for SARS-CoV-2 in the first and second waves, as the main objective of this work. Using the MESPC model, it was possible to estimate, with a certain reliability degree, the peak of contagion for the first and second waves in Brazil, with one day difference from the real to the forecast. It is possible to use MESPC to forecast the peak of contagion for several regions, provided that the necessary structure and calibration are respected.

The tobacco epidemic is one of the leading public health threats the world has ever faced and public health policy that seeks to limit the problem may not only have to target the price of tobacco but also the initiation stage in a smoker’s life – the adolescent stage. This research contributes to the health economics literature by using a Bayesian hierarchical logistic model, estimated using Hamiltonian Monte Carlo (HMC) methods to empirically identify what drives the intentions to quit smoking among adolescent smokers in Zambia. Results suggest that among the junior secondary school-going adolescent smokers in Zambia, about 63% have plans to quit smoking. We find socio-demographic characteristics and several tobacco-smoking-related factors as the main drivers of adolescent smokers’ plans to quit smoking. Most importantly, we provide insights that could be useful to help adolescent smokers realize their quitting plans.

During the recent decades, global warming by greenhouse gas evolution has attracted worldwide attention and evermore strict regulations thereon have become institutionalized as international policies. In the process, more environment-friendly power generation technologies have been developed utilizing fossil fuels with a view to timely commercialization. As one of such “clean coal” technology, Integrated Gasification Fuel Cell system is a promising power generation means and COS hydrolysis reactor is installed downstream of coal syngas to remove acidic gas constituents such as H2S and COS. The most significant design parameters affecting performance of the COS hydrolysis reactor were selected to be GHSV, (catalytic) reaction temperature, and length ratio and numerical modeling was performed considering heat and fluid flow transfer as well as chemical reaction kinetics. Effect of the selected design parameters on the variation of conversion rate and reactant gas mixture concentration were comprehensively investigated to predict performance of COS hydrolysis reactor. Stochastic modeling of reactor performance was finally performed using Monte Carlo simulation and linear regression fitting.

Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators involve a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function.

Layered halide perovskites have emerged as a promising contender in solid-state lighting thanks to their tunable bandgap, high luminescence quantum yield, and high color purity. However, the fabrication of perovskite light-emitting devices in laboratories usually experiences low de-vice-to-device reproducibility since perovskite crystallization is highly sensitive to ambient condi-tions, such as humidity and temperature. Although device processing inside gloveboxes is primarily used to reduce the influence of oxygen and moisture, several extraneous variables, including gas-eous pollution traces and thermal fluctuations in the inert atmosphere or contaminations from re-sidual solvents confined in the enclosed environment, can destabilize the crystallization process and alter the properties of the emissive layers. Here, we examine different experimental set-ups for spin-coating deposition of layered perovskites in inert atmospheres and their influences on the formation of polycrystalline thin films. Our results demonstrate that fluctuations in the glovebox properties (concentrations of residual O2 and H2O or solvent traces), even in very short timescales, can negatively impact the consistency of the perovskite film formation, while thermal variation plays a relatively minor role in this phenomenon. Furthermore, careful storage of chemical species inside the workstation is critical for reproducing high-quality perovskite layers. Consequently, when applying our most controlled environment for perovskite deposition, the photoluminescence life-time of perovskite thin films shows a standard deviation of only 3%, whereas the reference set-up yields a 15% standard deviation. Regarding complete perovskite light-emitting diodes, the uncer-tainties in statistical luminance and EQE data are significantly reduced from 230% and 140% to 38% and 42%, respectively.

We present evidence of the anchoring bias and disposition effect (DE) in investor trading decisions across different asset markets (stable and volatile assets) and market scenarios (stable and volatile markets), as well as how these biases affect investors' portfolio performance. The study employs a quasi-experimental design, allowing subjects to engage in interactive trading with four securities, two with potential negative returns and two with potential positive returns, within a simulated asset market. Our findings reveal the presence of the disposition effect and anchoring bias among individual investors in India. Second, market scenarios or volatility can influence the investors’ behavioral biases. The disposition effect was more prominent in volatile markets and the total market, while the anchoring bias was significant in stable and total markets. Additionally, investors exhibiting the disposition effect exhibit lower portfolio performance, while those demonstrating anchoring bias perform relatively better. The study's findings can help individual investors understand their behavioral biases, avoid them in future trading decisions, and improve portfolio performance. Financial institutions and regulatory agencies could identify investment products and supportive investor interactions towards better market trade decisions.

This paper proposes a parameter coordination optimization design of Power System Stabilizer (PSS) based on an improved Quasi-Affine Transformation Evolutionary (QUATRE) algorithm to suppress low-frequency oscillation and improve the dynamic stability of power system. To begin, the Simulated Annealing (SA) algorithm randomly updates the globally optimal solution of each QUATRE iteration and matches the inferior solution with a certain probability to escape the local extreme point. This new algorithm is first applied to power system. Second, Since damping ratio is one of the criteria to measure the dynamic stability of power system, this paper sets the objective function according to the principle of maximization of damping coefficient of electromechanical mode, and uses SA-QUATRE to search a group of global optimal PSS parameter combinations to improve the stability margin of the system as much as possible. Finally, the method's rationality and validity were validated by applying it to the simulation examples of IEEE 4-machine 2-area system with different operation states. The comparison with the traditional optimization algorithm shows that the proposed method has more advantages for multi-machine PSS parameter coordination optimization, and can restrain the low-frequency oscillation of power system more effectively and enhance the system stability.

Gradient based methods are popularly used in training neural networks and can be broadly categorized into first and second order methods. Second order methods have shown to have better convergence compared to first order methods, especially in solving highly nonlinear problems. The BFGS quasi-Newton method is the most commonly studied second order method for neural network training. Recent methods have shown to speed up the convergence of the BFGS method using the Nesterov’s acclerated gradient and momentum terms. The SR1 quasi-Newton method though less commonly used in training neural networks, are known to have interesting properties and provide good Hessian approximations when used with a trust-region approach. Thus, this paper aims to investigate accelerating the Symmetric Rank-1 (SR1) quasi-Newton method with the Nesterov’s gradient for training neural networks and briefly discuss its convergence. The performance of the proposed method is evaluated on a function approximation and image classification problem.

Splines and quasi-interpolation operators are important both in approximation theory and applications. In this paper, we construct a family of quasi-interpolation operators for the bivariate quintic spline spaces S₅³ (∆_mn⁽²⁾). Moreover, the properties of the proposed quasi-interpolation operators are studied, as well as its applications for solving two-dimensional Burgers’ equation and image reconstruction. Some numerical examples show that these methods, which are easy to implement, provide accurate results.

Seismic travel time tomography using surface waves is an effective tool for three-dimensional crustal imaging. Historically, these surface waves are the result of active seismic sources or earthquakes. More recently, however, also surface waves retrieved through the application of seismic interferometry are exploited. Conventionally, two-step inversion algorithms are employed to solve the tomographic inverse problem. That is, a first inversion results in frequency-dependent, two-dimensional maps of phase velocity, which then serve as input for a series of independent, one-dimensional frequency-to-depth inversions. As such, a two-dimensional grid of localized depth-dependent velocity profiles are obtained. Stitching these separate profiles together subsequently yields a three-dimensional velocity model. Relatively recently, a one-step three-dimensional non-linear tomographic algorithm has been proposed. The algorithm is rooted in a Bayesian framework using Markov chains with reversible jumps, and is referred to as transdimensional tomography. Specifically, the three-dimensional velocity field is parameterized by means of a polyhedral Voronoi tessellation. In this study, we investigate the potential of this algorithm for the purpose of recovering the three-dimensional surface-wave-velocity structure from ambient noise recorded on and around the Reykjanes Peninsula, southwest Iceland. To that end, we design a number of synthetic tests that take into account the station configuration of the Reykjanes seismic network. We find that the algorithm is able to recover the 3D velocity structure at various scales in areas where station density is high. In addition, we find that the standard deviation on the recovered velocities is low in those regions. At the same time, the velocity structure is less well recovered in parts of the peninsula sampled by fewer stations. This implies that the algorithm successfully adapts model resolution to the density of rays. Also, it adapts model resolution to the amount of noise on the travel times. Because the algorithm is computationally demanding, we modify the algorithm such that computational costs are reduced while sufficiently preserving non-linearity. We conclude that the algorithm can now be applied adequately to travel times extracted from (time-averaged) station-station cross correlations by the Reykjanes seismic network.

An important aspect of PET imaging in clinical application is the localization and detection of tumors and lesions by administering a predetermined amount of radiotracer. This allows for example, a detailed view of what is going inside the patient body in cellular level. The quality of PET image is strongly dependent on the amount of administrated radiotracer and the patient’s body parameters. As the amount of injection radiotracer increases, the quality of resulting image increases and the lesion detection efficiency increases. The PET examiner society recognizes that any dose of radiotracer is associated with some possible radiation risks. It can harmful to the patient if essential PET imaging session is not made due to fear of radiation risk. In order to ensure the highest quality diagnosis and the smallest radiation risk, the patient should receive the smallest amount of radiotracer that provide image with sufficient quality. Our study is focused on proposing an efficient PET simulation tool that predict the smallest possible amount of administrated radiotracer to provides the appropriate diagnostic information based on significant patient’s body parameters (weight, age) at fixed scanning to improve the clinical diagnostic process in term of tumor detecting and localization. We have built a model of particular PET scanner and model of patient based on real MRI image and digital anthropomorphic phantom of our region of interest (brain). We have performed Monte Carlo simulation for whole PET procedure with special parameter In evaluating stage, a dataset of 60 patient is used and 11 independent dose prediction simulation for each patient are performed. We conclude that our simulator perform a desirable and efficient prediction of injection radiotracer amount that optimizes the current clinical amount up to 28%. In addition, we found that the total injected radiotracer dosage for adult patients are mostly affected when considering patient weight rather than patient age.

The aim of this work is to study the zonally asymmetric stratopause that occurred in the Arctic winter of 2019/2020, when the polar vortex was particularly strong and there was no sudden stratospheric warming. Aura Microwave Limb Sounder temperature data were used to analyze the evolution of the stratopause with a particular focus on its zonally asymmetric wave 1 pattern. There was a rapid descent of the stratopause height below 50 km in the anticyclone region in mid-December 2019. The descended stratopause persisted until mid-January 2020 and was accompanied by a slow descent of the higher stratopause in the vortex region. The results show that the stratopause in this event was inclined and lowered from the mesosphere in the polar vortex to the stratosphere in the anticyclone. It was found that the vertical amplification of wave 1 between 50 km and 60 km closely coincides in time with the rapid stratopause descent in the anticyclone. Overall, the behavior contrasts with the situation during sudden stratospheric warmings when the stratopause reforms at higher altitudes following wave amplification events. We link the mechanism responsible for coupling between the vertical wave 1 amplification and this form of zonally asymmetric stratopause descent to the unusual disruption of the quasi-biennial oscillation that occurred in late 2019.

The process of senescence (aging) is largely determined by the action of wild-type genes. For most organisms, this does not reflect any adaptive function of senescence, but rather evolutionary effects of declining selection against genes with deleterious effects later in life. To understand aging requires an account of how evolutionary mechanisms give rise to pathogenic gene action and late-life disease, that integrates evolutionary (ultimate) and mechanistic (proximate) causes into a single explanation. A well-supported evolutionary explanation by G.C. Williams argues that senescence can evolve due to pleiotropic effects of alleles with antagonistic effects on fitness and late-life health (antagonistic pleiotropy, AP). What has remained unclear is how gene action gives rise to late-life disease pathophysiology. One ultimate-proximate account is T.B.L. Kirkwood’s disposable soma theory. Based on the hypothesis that stochastic molecular damage causes senescence, this reasons that aging is coupled to reproductive fitness due to preferential investment of resources into reproduction, rather than somatic maintenance. An alternative and more recent ultimate-proximate theory argues that aging is largely caused by programmatic, developmental-type mechanisms. Here ideas about AP and programmatic aging are reviewed, particularly those of M.V. Blagosklonny (the hyperfunction theory) and J.P. de Magalhães (the developmental theory), and their capacity to make sense of diverse experimental findings is described.

Sumudu transform of the Dixon elliptic function with non zero modulus a ≠ 0 for arbitrary powers sm^N(x,a) ; N ≥ 1 ; sm^N(x,a)cm(x,a) ; N ≥ 0 and sm^N(x,a)cm²(x,a) ; N ≥ 0 is given by product of Quasi C fractions. Next by assuming denominators of Quasi C fraction to 1 and hence applying Heliermann correspondance relating formal power series (Maclaurin series of Dixon elliptic functions) and regular C fraction, Hankel determinants are calculated and showed by taking a = 0 gives the Hankel determinants of regular C fraction. The derived results were back tracked to the Laplace transform of sm(x,a) ; cm(x,a) and sm(x,a)cm(x,a).

Generalized distributions have been studied a lot recently because of their flexibility and relia-bility in modeling lifetime data. The two-parameter exponentially-modified logistic distribution is a new flexible modified distribution that was introduced recently. It is regarded as a strong competitor for widely used classical symmetrical and non-symmetrical distributions such as normal, logistic, lognormal, and log-logistic. In this study, the unknown parameters of the dis-tribution are estimated using the maximum likelihood method. The Grey Wolf Optimization al-gorithm, which is a new meta-heuristic algorithm, is applied in order to solve the nonlinear like-lihood equations of the study model. The performance of the Grey Wolf Optimization method is compared to that of the other meta-heuristic algorithms used in this study, which include the Whale Optimization Algorithm, the Sine Cosine Algorithm, and the Particle Swarm Optimiza-tion Algorithm. The efficiencies of maximum likelihood estimates for all algorithms are com-pared via an extensive Monte-Carlo simulation study. The likelihood estimates for the location α and scale β parameters of the exponentially-modified logistic distribution developed with the Grey Wolf Optimization algorithm are the most efficient among others, according to simulation findings. Four real datasets are analyzed to show the flexibility of this distribution.

The entry and exit method is a simple and practical method to determine the critical slip surface of slope. However, it has the drawback of sacrificing computational efficiency to improve search accuracy. To solve this problem, this paper proposes an improved entry and exit method to search for the critical slip surface. Based on the random fields generated by using Karhunen–Loève expansion method, the simplified Bishop’s method combined with the improved entry and exit method is used to determine the critical slip surface and its corresponding minimum factor of safety. Then, the failure probability is calculated by conducting Monte Carlo simulation. Two examples are reanalyzed to verify the accuracy and efficiency of the proposed method. Meaningful comparisons are made to demonstrate the calculating accuracy and calculating efficiency of the improved entry and exit method in searching for the minimum safety factor of slope, based on which the effect of the reduced searching range on slope reliability was explored. The results indicate that the proposed method provides a practical tool for evaluating the reliability of slopes in spatially variable soils. It can greatly improve the computational efficiency in relatively high-computational accuracy of slope reliability analysis.

In the field of pavement engineering, the determination of the mechanical characteristics is one of the essential process for reliable material design and highway sustainability. Early determination of mechanical characteristics of pavement is highly essential for road and highway construction and maintenance. Tensile strength (TS), compressive strength (CS) and flexural strength (FS) of roller compacted concrete pavement (RCCP) are very crucial characteristics as they are necessitated for many data from mixture proportions as input variables. In this research, the classification-based regression models named Random Forest (RF), M5rule model tree (M5rule), M5prime model tree (M5p) and Chi-square Automatic Interaction Detection (CHAID) are developed for simulation of the mechanical characteristics of RCCP. A comprehensive and reliable dataset comprising 621, 326 and 290 data records for CS, TS and FS experimental cases extracted from several open sources over the literature. The mechanical properties are developed based on influential inputs combination that processed using Principle Component Analysis (PCA). The applied PCA method as feature selection is specified that volumetric/weighted content forms of experimental variables (e.g., coarse aggregate, fine aggregate, supplementary cementitious materials, water and binder) and specimens’ age are the most effective inputs to generate the better performances. Several statistical metrics are measured to evaluate proposed classification-based regression models. RF model revealed an optimistic classification capacity of the CS, TS and FS prediction of the RCCP in comparison with the CHAID, M5rule, and M5p models. The research is extended for the results verification using Monte-carlo model for the uncertainty and sensitivity of variables importance analysis. Overall, the proposed methodology indicated a reliable soft computing model that can be implemented for the material engineering construction and design.

Confidence interval of is an interval corresponding to a specified confidence and including the true value. It can be used to describe the precision of a statistical quantity and quantify its uncertainty. Although the principle of maximum entropy (POME) has been used for a variety of applications in hydrology, the confidence intervals of the POME quantile estimators have not been available. In this study, the calculation formulas of asymptotic variances and confidence intervals of quantiles based on POME for Gamma, Pearson type 3 (P3) and Extreme value type 1 (EV1) distributions were derived. Monte Carlo Simulation experiments were performed to evaluate the performance of derived formulas for finite samples. Using four data sets for annual precipitation at the Weihe River basin in China, the derived formulas were applied for calculating the variances and confidence intervals of precipitation quantiles for different return periods and the results were compared with those of the methods of moments (MOM) and of maximum likelihood (ML) method. It is shown that POME yields the smallest standard errors and the narrowest confidence intervals of quantile estimators among the three methods, and can reduce the uncertainty of quantile estimators

The properties of rubber are strongly influenced by the distribution of filler within the polymer matrix. Here we introduce a Monte Carlo-based morphology generator. The basic elements of our model are cubic cells, which, in the current version, can be either silica filler particles or rubber volume elements in adjustable proportion. The model allows the assignment of surface free energies to the particles according to whether a surface represents, for instance, 'naked' silica or silanised silica. The amount of silanisation is variable. We use a nearest-neighbour site-exchange Monte Carlo algorithm to generate filler morphologies, mimicking flocculation. Transmission electron micrographs (TEM) as well as small angle scattering (SAS) intensities can be calculated along the Monte Carlo trajectory. In this work we demonstrate the application of our morphology generator in terms of selected examples. We illustrate its potential as a tool for screening studies, relating interface tensions between the components to filler network structure as characterized by TEM and SAS.

A simple bounded stochastic motion is employed for modelling equity price dynamics under shocks, where the prices are constrained to lie between two positive bounds. While the stochastic process has an inaccessible boundary, the other boundary is quasi-bounded, implying that the boundary can be breached if the probability leakage condition is met. The quasi-boundedness of the process at the boundary can thus provide an indicator of possible risk of equities under price shocks. Empirical calibration of model parameters of the proposed process for equities can be performed easily due to the availability of an analytically tractable probability density function which generates fat-tailed distributions consistent with empirical observations. The volatility and mean-reversion of the S&P500 dynamics calibrated by the process is positively and negatively co-integrated respectively with the VIX index representing the level of market distress. The process captures the high likelihood of Hertz’s default about two months earlier, using only information until that point, and before the firm filed for Chapter 11 bankruptcy in May 2020 as a result of the Covid-19 pandemic. Empirical calibration of the process for GameStop’s stock price shows that the short squeeze of the stock occurred when the condition for breaching an upper boundary was met on 14 January 2021, i.e., about two weeks before major short sellers closed out their positions with significant losses. The trading volume of the stock was positively co-integrated with the probability leakage ratio.

The purpose of this work is to study quasi-stationary wave structure in the mid-latitude stratosphere and mesosphere (40–50°N) and its role in the formation of the annual ozone cycle. Geopotential height and ozone from Aura MLS data are used and winter climatology for January–February 2011–2020 is considered. More closely examined is the 10-degree longitude segment centered on Longfengshan Brewer station, China, and located in the region of the Aleutian Low influence associated with the quasi-stationary zonal maximum of total ozone. Annual and semi-annual oscillations in ozone were compared using units of ozone volume mixing ratio and concentration, as well as changes in ozone peak altitude and in time series of ozone at individual pressure levels between 316 hPa (9 km) and 0.001 hPa (96 km). The ozone maximum in the vertical profile is higher in volume mixing ratio (VMR) values than in concentration by about 15 km (5 km) in the stratosphere (mesosphere), in consistency with some previous studies. We found that the properties of the annual cycle are better resolved in the altitude range of the main ozone maximum: middle–upper stratosphere in VMR and lower stratosphere in concentration. Both approaches reveal SAO/AO-related changes in the of ozone peak altitudes in a range of 4–6 km during the year. In the lower-stratospheric ozone of the Longfengshan domain, an earlier development of the annual cycle takes place with a maximum in February and a minimum in August compared to spring and autumn, respectively, in zonal means. This is presumably due to the higher rate of dynamical ozone accumulation in the region of the quasi-stationary zonal ozone maximum. The “no-annual-cycle” transition layers are found in the stratosphere and mesosphere. These layers with undisturbed ozone volume mixing ratio throughout the year are of interest for more detailed future study.

Firstly, a new set of Quasi-Cubic Trigonometric Bernstein basis with two tension shape parameters is constructed, and we prove that it is an optimal normalized totally basis in the framework of Quasi Extended Chebyshev space. And the Quasi-Cubic Trigonometric Bézier curve is generated by the basis function and the cutting algorithm of the curve are given, the shape features (cusp, inflection point, loop and convexity) of the Quasi-Cubic Trigonometric Bézier curve are analyzed by using envelope theory and topological mapping; Next we construct the non-uniform Quasi-Cubic Trigonometric B-spline basis by assuming the linear combination of the optimal normalized totally positive basis have partition of unity and continuity, and its expression is obtained. And the non-uniform B-spline basis is proved to have totally positive and high-order continuity. Finally, the non-uniform Quasi Cubic Trigonometric B-spline curve and surface are defined, the shape features of the non-uniform Quasi-Cubic Trigonometric B-spline curve are discussed, and the curve and surface are proved to be continuous.

Emerging semiconductor devices often lack accurate analytical models. The same is usually true of any devices working under extreme conditions like cryogenic temperatures. The usual workaround involves the use of approximation models, usually based on lookup tables or neural networks individually fitted to measurement data. In the case of experimental devices or ones working under extreme conditions, the number of units available for measurement is limited. As a result, the number of approximation-model instances is too small to enable a statistical simulation of even middle-sized circuits, which is a necessary step in integrated-circuit design since it provides the realistic picture of the circuit’s behavior in the presence of manufacturing process variations. Approximation models using structure parameters as inputs do exist in the literature, but are only useful if the end user knows the statistical distributions of those parameters, which is not usually the case. We propose a technique based on generative machine learning, namely the variational autoencoder, that uses only a small sample of devices to capture the essential features of their I–V curves under process variations and subsequently generates an arbitrary number of similarly disturbed curves. The model trained on as few as 20 instances per device type is shown to precisely reproduce the distributions of period and power consumption of a ring oscillator.

The assessment of existing bridges constitutes a serious challenge for the researchers since their safety evaluation is a complex task due to several sources of uncertainty affecting the parameters identified in the knowledge process. In particular, a high level of uncertainties may arise in the case of prestressing systems requiring special investigation strategies and experimental technologies. The paper aims to investigate how the knowledge uncertainties influence the safety assessment of existing post-tensioned bridges. The semi-probabilistic approach is used, and uncertainties related to the knowledge process are described considering the number of tests, the confidence level assumed for the estimation of the problem parameters, and the measurement errors of in-situ tests. Results concerning a typical post-tensioned bridge are reported and the variability of the outcomes is discussed and compared with a reference case of “perfect knowledge”. The main results of the study confirm the general robustness of the semi-probabilistic Eurocodes safety format and provide an overview of dispersion due to different choices about the number of tests and confidence level. The measurement errors may lead to significant under-estimation of the computed structural capacity with respect to the reference one.

Heterogeneous systems of limited capacity have general applications in manufacturing, but also for logistic or service systems, due to the differences in server or workstation performance or work assignment, in close relationship with system flexibility, where saturation and blocking are ordinary situations of systems with high demand and limited capacity, so accurate loss quantification is essential for performance evaluation. Multi-class systems of limited capacity have been studied much less than parallel homogeneous systems (Erlang models). In this context, accurate models for parallel heterogeneous ordered-entry systems are developed: without any prior queue, M/Mi/c/c, and with a k capacity queue M/Mi/c/c+k. These new matrix models give an exact state formulation, and their accuracy is verified through discrete event simulation and comparison with literature results. Also, the effect of queue capacity is studied in relationship with the pattern of service rates. Next, the heterogeneous recirculating system model is also developed with good approximation results. Finally, the proposed models are applied to evaluate systems with non-exponential service times, through a new hybrid methodology by combining the Markovian model and the Monte Carlo Method (MCM) for normal or lognormal service times that also yield useful good approximations to the simulated system.

The concrete sector is known for its significant contribution to the CO2 emissions. There are two main contributing factors for this situation: the large amount of concrete consumed per year on the planet and the high CO2 released from Portland cement manufacture, the key binding agent in concrete. To face the consequent sustainability issues, diverse strategies have been explored on the carbon capture and storage potential of cementitious materials. This paper addresses the potential of storing CO2 in concrete during the curing stage, applied to the precast Portuguese industry. To this purpose, it was assumed that CO2 will become a waste that will require an outlet in the future, considering that carbon capture will become mandatory in many industries. This work concluded that the net benefit in terms of carbon retention is positive for the process of storing carbon in concrete during the curing stage. More specifically, it was demonstrated that the additional emissions from the introduction of this new operation are only 10% of the stored amount, returning a storage potential of 76 000 tonnes of CO2 yearly. Moreover, the overall net reduction in the concrete life cycle averages 9.4% and 8.8% for precast elements and only non-structural elements, respectively. When a low cement dosage strategy is coupled with carbonation curing technology, the overall carbon net reduction is estimated to be 45%.

Aiming at the problem of hypersonic morphing vehicle avoiding no-fly zones and reaching the target, an improved predictor-corrector guidance method is proposed. Firstly, the aircraft motion model and the constraint model are established. Then, the basic algorithm is given, the Q-learning method is used to design the attack angle and sweep angle scheme to ensure that the aircraft can fly over the low-altitude zones. The B-spline curve is used to design the location of flight path points and the bank angle scheme is designed according to the predictor-corrector method, so that the aircraft can fly around to avoid high-altitude zones. Next, Monte Carlo reinforcement learning(MCRL) method is used to improve predictor-corrector method and Deep Neural Network(DNN) is used to fit reward function. The improved method can generate trajectory with better performance. Simulation results verify the effectiveness of the proposed algorithm.

Sleep deficiency impacts the quality of life and may have serious health consequences in the long run. Questionnaire-based subjective assessment of sleep deficiency has many limitations. On the other hand, objective assessment of sleep deficiency is challenging. In this study, we propose a polysomnography-based mathematical model for computing baseline sleep deficiency severity score and then investigated the estimation of sleep deficiency severity using features available only from wearable sensor data including heart rate variability and single-channel electroen-cephalography for a dataset of 500 subjects. We used Monte-Carlo Feature Selection (MCFS) and inter-dependency discovery for selecting the best features and removing multi-collinearity. For developing the Regression model we investigated both the frequentist and the Bayesian ap-proaches. An Artificial Neural Network achieved the best performance of RMSE = 5.47 and an R-squared value of 0.67 for sleep deficiency severity estimation. The developed method is com-parable to conventional methods of Functional Outcome of Sleep Questionnaire and Epworth Sleepiness Scale for assessing the impact of sleep apnea on sleep deficiency. Moreover, the results pave the way for reliable and interpretable sleep deficiency severity estimation using a wearable device in Smart Health.

We review and assess the classic SIR and SEIR epidemiological models regarding possible applications to the COVID-19 pandemic. In spite of numerous more complicated models, we show how the qualitative features of the solution to the SIR and SEIR models continue to provide valuable public health insights in some scenarios. Using estimated COVID-19 data as of this date, the SEIR model shows that if it were possible to reduce R₀ from 2.5 to 1.25 through social distancing and other measures, the maximum fraction of the population that would become infected at any particular time would drop from 17% to 4%, provided that all of the model assumptions are satisfied. Finally, we compare the classic SIR model with a recent stochastic model with favorable results. Since this comparison underscores the importance of underlying connectivity assumptions, we conclude with Monte-Carlo simulations with specific connectivity that reproduce the classical SIR model with standard incidence.

We estimated the daily intakes of chlorpyrifos and cypermethrin via ingestion of indoor dust and outdoor soils using the Stochastic Human Exposure and Dose Simulation Model on a probabilistic approach for Taiwan’s homes. Variable information for the daily intake estimation, such as concentration, ingestion rate, body weight, was adopted from previous studies. Monte Carlo simulation was performed with 1,000,000 iterations to simulate a single daily intake, which was shown in terms of percentage of the Acceptable Daily Intake (ADI) of either insecticide. The daily intakes were minimal with a 99% probability; at 99.9th percentile, however, the total intakes leaped to 13.1% and 20.0% of the respective ADIs of chlorpyrifos and cypermethrin. The sensitivity analysis indicated that concentration was the most determinant variable. Compared to the data of daily intakes via dietary ingestion of vegetables derived from a previous study, the estimated intakes by this study were considerable at the highest percentile, which referred to insecticide residues few days after insecticide application. Consequently, the non-dietary ingestion exposure to either insecticide was negligible in most cases; nevertheless, for those with indoor insecticide applications, the daily intakes for young children could be of concern. Frequently home cleaning is recommended to reduce the exposure.

Bayesian analysis for masked data under competing risk frameworks is studied for the purpose of assessing the impact of covariates on the hazard functions when the failure time is exactly observed for some subjects but only known to lie in an interval of time for the remaining subjects. Such data, known as partly interval-censored data, usually result from periodic inspection. Dirichlet and Gamma processes are assumed as priors for masking probabilities and baseline hazards. The Markov Chain Monte Carlo (MCMC) technique is employed for the implementation of the Bayesian approach. The effectiveness of the proposed model is tested through numerical studies, including simulated and real data sets.

Induction heating (IH) technology is widely recognized and utilized in residential applications due to its high efficiency and safe operating characteristics. Resonant inverter circuits are widely used in IH systems because of their high efficiency and ability to perform soft-switching. Among the various resonant inverters used in IH systems, the single-switch quasi-resonant (SSQR) inverter topology is typically preferred for low-cost and low-output-power applications. Despite its cost advantage, the SSQR topology has a relatively narrow soft-switching range, which can be unstable depending on the electrical parameters of the load and the resonant converter circuit. Accurately determining the capacitance value of the resonant capacitor and the inductance value of the induction coil, which are the key circuit elements of the SSQR induction cooker, is crucial for designing a reliable, efficient, and durable cooking system. In other words, there exists a critical relationship between the resonant converter circuit parameters, load characteristics, and safe operating conditions. Additionally, when considering closed-loop control methods used for power control and safety, selecting appropriate resonant circuit elements becomes vital in ensuring both reliable and efficient operation. This paper focuses on a novel and simplified design method for the SSQR inverter utilized in household appliances. The proposed method and its advantages in terms of the safe operating area of the switch are theoretically investigated and verified through simulations and prototype circuits.

Tick-borne encephalitis virus (TBEV) is a pathogenic, enveloped, positive-stranded RNA virus in the family Flaviviridae. Structural studies of flavivirus virions have primarily focused on mosquito-borne species with only one cryo-electron microscopy (cryo-EM) structure of a tick-borne species published. Here, we present a 3.3 Å cryo-EM structure of the TBEV virion of the Kuutsalo-14 isolate, confirming the overall organisation of the virus. We observe conformational switching of the peripheral and transmembrane helices of M protein, which can explain the quasi-equivalent packing of the viral proteins and highlights their importance in stabilizing the membrane protein arrangement in the virion. The residues responsible for the M protein inter-actions are highly conserved in TBEV but not in the structurally studied Hypr strain, nor in mosquito-borne flaviviruses. These interactions may compensate for the lower number of hydrogen bonds between E proteins in TBEV compared to the mosquito-borne flaviviruses. The structure reveals two lipids bound in the E protein, which are important for virus assembly. The lipid pockets are comparable to those recently described in mosquito-borne Zika, Spondweni, Dengue, and Usutu viruses. Our results thus advance the understanding of tick-borne flavivirus architecture and virion-stabilising interactions.

A system of simultaneous multi-variable nonlinear equations can be solved by the Newton’s method with local q-quadratic convergence if the Jacobian is analytically available. If this is not the case, then quasi-Newton methods with local q-superlinear convergence give solutions by approximating the Jacobian in some way. Unfortunately, the quasi-Newton condition (secant equation) doesn’t completely specify the Jacobian approximate in multi-dimensional case, so its full-rank update is not possible with classic variants of methods. The suggested new iteration strategy (“T-Secant”) allows full-rank update of the Jacobian approximate in each iteration by determining two independent approximates for the solution. They are used to generate a set of new independent trial approximates, then the Jacobian approximate can fully be updated. It is shown, that the T-Secant approximate is in the vicinity of the classic quasi-Newton approximate, providing that the solution is evenly surrounded by the new trial approximates. The suggested procedure increases the super-linear convergence of the secant method (φS = 1.618...) to super-quadratic (φT = φS + 1 = 2.618...) and the quadratic convergence of the Newton-method (φN = 2) to cubic (φT = φN + 1 = 3) in one dimensional case. The Broyden-type efficiency (mean convergence rate) of the suggested method in multi-dimensional case is an order higher than the efficiency of other classic low-rank update quasi-Newton methods as shown by numerical examples on a Rosenbrock-type test-function with up to 1000 variables. The geometrical representation (hyperbolic approximation) in single variable case helps explaining the basic operations and a vector-space description is also given in multi-variable case.

The paper proposes an approach to determine the optimal number of medical specialists in the region. According to the author's theoretical model, in order to maximize public welfare, the marginal contribution of the last physician recruited to the growth of the public utility function should be equal to the marginal cost of attracting him and providing conditions for his work. To empirically assess the contribution of physicians to the number of lives saved, the CVD mortality rate is modeled using the instrumental variable method. On average, a 1% increase in the number of cardiologists is associated with a 0.6% decrease in cardiovascular mortality. At the level of provision with cardiologists in the amount of 1 per 100 thousand people, their marginal contribution to the number of lives saved is not less than 129 per 100 thousand people, with a further decrease by 5 people with the increase in the level of provision by 1 unit. The use of the obtained results will increase the validity of managerial decisions and in determining the optimal level of doctors when choosing between alternative possibilities of spending money on hiring doctors of different profiles or other expenses, especially in the case of limited resources.

The pore-throat radius of shale oil reservoir is extremely small, and it is difficult to accurately obtain the absolute permeability and oil-water two-phase relative permeability of actual oil reservoir through conventional core experiments. However, these parameters are very important for reservoir numerical simulation. In this paper, a method for characterizing flow parameters based on pore network model that considers differential pressure flow and diffusion flow is proposed. Firstly, a digital core was reconstructed using focused ion beam scanning electron microscopy (FIB-SEM) from the Gulong shale reservoir in the Songliao Basin, China, and a pore network model was extracted. Secondly, quasi-static single-phase flow and two-phase flow equations considering diffusion were established in the pore network model. Finally, pore throat parameters, absolute permeability, and oil-water two-phase permeability curves were calculated, respectively. The results show that the pore throat distribution of Gulong shale reservoir is mainly concentrated in the nanometer scale, the mean pore radius is 87 nm, the mean throat radius is 41 nm, and the mean coordination number is 3.97; The calculated permeability considering diffusion is 0.000124mD, which is approximately twice the permeability calculated without considering diffusion; The irreducible water saturation of the Gulong shale reservoir is approximately 0.26, and the residual oil saturation is approximately 0.73. The method proposed in this paper can provide an important approach for characterizing the flow parameters of similar shale oil reservoirs.

A quasi-two-dimensional (quasi-2D) modelling approach is introduced to mimic transverse mixing of an inflow into a river from one of its banks, either an industrial outfall or a tributary. The concentrations of determinands in the inflow vary greatly from those in the river, leading to very long mixing lengths in the river downstream of the inflow location. Ideally, a two-dimensional (2D) model would be used on a small scale to capture the mixing of the two flow streams. However, for large-scale applications of several hundreds of kilometres of river length, such an approach demands too many computational resources and too much computational time, especially if the application will at some point require ensemble input from climate-change scenario data. However, a one-dimensional (1D) model with variables varying in the longitudinal flow direction but averaged across the cross-sections is too simple of an approach to capture the lateral mixing between different flow streams within the river. Hence, a quasi-2D method is proposed in which a simplified 1D solver is still applied but the discretisation of the model setup can be carried out in such a way as to enable a 2D representation of the model domain. The quasi-2D model setup also allows secondary channels and side lakes in floodplains to be incorporated into the discretisation. To show proof-of-concept, the approach has been tested on a stretch of the lower Athabasca River in Canada flowing through the oil sands region between Fort McMurray and Fort MacKay. A dye tracer and suspended sediments are the constituents modelled in this test case.

Current industrial trends bring new challenges in energy absorbing systems. Polymer materials as the traditional packaging material seem to be promising due to their low weight, structure and production price. Based on the review, the linear low-density polyethylene material was identified as the most promising material for absorbing impact energy. The current paper addresses the identification of the material parameters and the development of a Constitutive material model to be used in future design by virtual prototyping. The paper deals with the experimental measurement of the stress-strain relations of the linear low-density polyethylene under static and dynamic loading. The quasi-static measurement is realized in two perpendicular principal directions and is supplemented by a test measurement in the 45 degrees direction, i.e. exactly between the principal directions. The quasi-static stress-strain curves are analyzed as an initial step for dynamic strain rate dependent material behavior. The dynamic response is tested in the drop tower using a spherical impactor hitting the flat material multi-layered specimen at two different energy levels. The strain rate dependent material model is identified by optimizing the static material response obtained in the dynamic experiments. The material model is validated by the virtual reconstruction of the experiments and by comparing the numerical results to the experimental ones.

Solidification microstructure is formed under high cooling rates and temperature gradients in powder-based additive manufacturing. In this study, a non-equilibrium multi-phase field method (MPFM), which was based on a finite interface dissipation model proposed by Steinbach et. al., coupled with a CALPHAD database was developed for a multicomponent Ni alloy. A qua-si-equilibrium MPFM was also developed for comparison. Two-dimensional equiaxed micro-structural evolution for the Ni (Bal.)–Al–Co–Cr–Mo–Ta–Ti–W–C alloy was performed at various cooling rates. The temperature–γ fraction profiles obtained under 10^5 K/s using non- and qua-si-equilibrium MPFMs were in good agreement with each other. Over 10^6 K/s, the differences between non- and quasi-equilibrium methods grew as the cooling rate increased. The non-equilibrium solidification was strengthened over a cooling rate of 10^6 K/s. Colum-nar-solidification microstructural evolution was performed under cooling rates from 5×10^5 K/s to 1×10^7 K/s at various temperature gradient values under the constant interface velocity (0.1 m/s). The results showed that as the cooling rate increased, the cell space decreased in both methods, and the non-equilibrium MPFM agreed well with experimental measurements. Our results show that the non-equilibrium MPFM can simulate solidification microstructure in powder bed fusion additive manufacturing.

We study in this work a dilute granular gas immersed in a thermal bath made of smaller particles with nonnegligible masses as compared with the granular ones. The kinetic theory for this system is developed and described by an Enskog--Fokker--Planck equation for the one-particle velocity distribution function. Granular particles are assumed to have inelastic and hard interactions, losing energy in collisions as accounted by a constant coefficient of normal restitution. The interaction with the thermal bath is based on a nonlinear drag force plus a white-noise stochastic force. To get explicit results of the temperature aging and steady states, Maxwellian and first Sonine approximations are developed. The latter takes into account the coupling of the excess kurtosis with the temperature. Theoretical predictions are compared with direct simulation Monte Carlo and event-driven molecular dynamics simulations. While good results for the granular temperature are obtained from the Maxwellian approximation, a much better agreement, especially as inelasticity and drag nonlinearity increase, is found when using the first Sonine approximation. The latter approximation is, additionally, crucial to account for memory effects like Mpemba and Kovacs-like ones.

We study in this paper the time evolution of stock markets using a statistical physics approach. We consider an ensemble of agents who sell or buy a good according to several factors acting on them: the majority of the neighbors, the market ambiance, the variation of the price and some specific measure applied at a given time. Each agent is represented by a spin having a number of discrete states q or continuous states, describing the tendency of the agent for buying or selling. The market ambiance is represented by a parameter T which plays the role of the temperature in physics: low T corresponds to a calm market, high T to a turbulent one. We show that there is a critical value of T, say T_c, where strong fluctuations between individual states lead to a disordered situation in which there is no majority: the numbers of sellers and buyers are equal, namely the market clearing. The specific measure, by the government or by economic organisms, is parameterized by $H$ applied on the market at the time t₁ and removed at the time t₂. We have used Monte Carlo simulations to study the time evolution of the price as functions of those parameters. In particular we show that the price strongly fluctuates near T_c and there exists a critical value H_c above which the boosting effect remains after H is removed. Our model replicates the stylized facts in finance (time-independent price variation), volatility clustering (time-dependent dynamics) and persistent effect of a temporary shock. The second party of the paper deals with the price variation using a time-dependent mean-field theory. By supposing that the sellers and the buyers belong to two distinct communities with their characteristics different in both intra-group and inter-group interactions, we find the price oscillation with time. Results are shown and discussed.

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this mini-review, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions and, on other hand, highly simplified structure-based Elastic Network Models of protein flexibility. In contrast to classical Molecular Dynamics the modeling strategies discussed here allow quite accurate modeling of much larger systems and longer time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Social recommendation is almost as the integration of the business platform and social platform, and gradually become a top in recommendation system. Social recommendation algorithm solves the problem of cold start and data sparseness for traditional commodity, while the internal structure of the relationship graph in social relations has not been fully excavated. This paper proposes two models of Micro Relation Transfer Model and Macro Relation Transfer Model of social relations, and applies the social relations transfer models into the social recommendation system. A relationship graph is built from the relationship between customers on the Internet. Micro Relation Transfer Model establishes the transfer activation function by calculating the relationship between the two customers using the similarity of interests set. Micro Relation Transfer Model spreads the relationship of friends by calculating the proportion of common neighbors held by the customer's social relations. In order to effectively control the transmission range and effect of social relations graph, we introduce pruning algorithm based on Monte Carlo Decision Tree convergence algorithm. The experimental results show that SRRTC algorithm enhances the success rate and stability significantly.

Background: Pediatric Rapid Response Teams (RRTs) are expected to significantly lower pediatric mortality in healthcare settings. This study evaluates RRTs' effectiveness in decreasing cardiac arrests and unexpected Pediatric Intensive Care Unit (PICU) admissions. Methods: A quasi-experimental study (2014-2017) at King Abdulaziz University Hospital, Jeddah, Saudi Arabia, involved 3,261 pediatric inpatients, split into pre-intervention (1,604) and post-intervention (1,657) groups. The intervention entailed pediatric specialists collaboratively managing acute decompensations following set protocols. Results: RRT was activated 471 times, primarily for respiratory distress (29.30%), sepsis (22.30%), clinical anxiety (13.80%), and hematological abnormalities (6.7%). Family concerns triggered 0.1% of activations. Post-RRT implementation, unplanned PICU admissions significantly reduced (RR = 0.552, 95% CI 0.485-0.628, P < 0.0001), and non-ICU cardiac arrests were eliminated (RR = 0). Patient care improvement was notable, with a -9.61 coefficient for PICU admissions (95% CI: -12.65 to -6.57, P < 0.001) and a -1.641 coefficient for non-ICU cardiac arrests (95% CI: -2.22 to -1.06, P < 0.001). Sensitivity analysis showed mixed results for PICU admissions, while zero-inflation Poisson analysis confirmed a reduction in non-ICU arrests. Conclusions: The deployment of pediatric RRTs is associated with fewer unexpected PICU admissions and non-ICU cardiopulmonary arrests, indicating improved PICU management. Further research using robust scientific methods is necessary to conclusively determine RRTs' clinical benefits.

To improve the understanding of soil–pile interaction of laterally loaded piles, the quasi-static test of pile group foundation in silty clay with 5 foundation piles arranged in cross shape is carried out. The direct measurement-based p–y curves at different depths showed nonlinear behavior, in which the initial stiffness and ultimate soil reaction increased as the depth increased. The loading displacement amplitude corresponding to the failure of each foundation pile in the high axial compression reinforced concrete pile group foundation is given as follows: the front-end nappe foundation pile (#3 pile) is 6mm, the middle pile is 16mm, the front-end nappe foundation pile is 20mm, and the loading displacement amplitude of side pile is 24mm. In silty clay, due to the existence of pile group effect, the edge effect has little effect on the horizontal foundation bed coefficient of foundation pile in pile group, while the influence of masking effect on the horizontal foundation bed coefficient of each foundation pile in pile group is not related to the position of pile, but only to the existence of masking effect, and the influence of different columns of masking effect on the horizontal foundation bed coefficient is the same.

Abstract: A new method for constructing exact solutions of the General Relativity equations for a dusty matter with fractal property is proposed. This method allows to find the solution of the GR-equations in terms of matter velocities : the connection coefficients and the Ricci tensor of space-time are expressed in terms of matter velocities; the metric tensor and the matter density are found as functions of velocity from the GR-equations. The connection coefficients and the Ricci tensor are invariant with respect to the discrete scaling transformation of velocity , where is constant. Therefore, the found solution can be used to simulate the fractal properties of the cosmic web in terms of matter velocities. This solution includes isotropic and anisotropic distributions of matter density. In an isotropic case, there is a class of exact solutions including both the well-known Friedmann’s solution and a solution with a periodic distribution of the matter density in space. This last solution may be used to simulate the quasi-periodic distribution of matter in the cosmic web. It is possible that, the cosmic web and its fractal properties are the space-time primary properties. These properties are described with a deformation tensor of the space-time.

This paper deals with a non-standard finite difference scheme defined on a quasi-uniform mesh for approximate solutions of the Magneto-Hydro Dynamics (MHD) boundary layer flow of an incompressible fluid past a flat plate for a wide range of the magnetic parameter. We show how to improve the obtained numerical results via a mesh refinement and a Richardson extrapolation. The obtained numerical results are favourably compared with those available in the literature.

Waterbomb structures are origami-inspired deformable structural components used in new types of robots. They have a unique radially deployable ability that enables robots to better adapt to their environment. In this paper, we propose a series of new waterbomb structures with square, rectangle, and parallelogram base units. Through quasi-static axial and radial compression experiments and numerical simulations, we prove that the parallelogram waterbomb structure has a twist displacement mode along the axial direction. Compared with the square waterbomb structure, the proposed optimal design of the parallelogram waterbomb structure reduces the critical axial buckling load-to-weight ratio by 55.4% and increases the radial stiffness-to-weight ratio by 67.6%. The significant increase in the radial stiffness-to-weight ratio of the waterbomb structure and decrease in the critical axial buckling load-to-weight ratio make the proposed origami pattern attractive for practical robotics applications.

In this work we have critically reviewed the processes in high-temperature sublimation growth of graphene in Ar atmosphere using enclosed graphite crucible. Special focus is put on buffer layer formation and free charge carrier properties of monolayer graphene and quasi-freestanding monolayer garphene on 4H-SiC. We show that by introducing Ar at different temperatures, TAr one can shift to higher temperatures the formation of the buffer layer for both n-type and semi-insulating substrates. A scenario explaining the observed suppresed formation of buffer layer at higher TAr is proposed and discussed. Increased TAr is also shown to reduce the sp3 hybridization content and defect densities in the buffer layer on n-type conductive substrates. Growth on semi-insulating substrates results in ordered buffer layer with significantly improved structural properties, for which TAr plays only a minor role. The free charge density and mobility parameters of monolayer graphene and quasi-freestanding monolayer graphene with different TAr and different environmental treatment conditions are determined by contactless terahertz optical Hall effect. An efficient annealing of donors on and near the SiC surface takes place in intrinsic monolayer graphene grown at 2000∘C, and which is found to be independent of TAr. Higher TAr leads to higher free charge carrier mobility parameters in both intrinsically n-type and ambient p-type doped monolayer graphene. TAr is also found to have a profound effect on the free hole parameters of quasi-freestanding monolayer graphene. These findings are discussed in view of interface and buffer layer properties in order to construct a comprehensive picture of high-temperature sublimation growth and provide guidance for growth parameters optimization depending on the targeted graphene application.

The independent set, IS, on a graph $G=(V,E)$ is $V^{*}\subseteq V$ such that no two vertices in $V^{*}$ have an edge between them. The MIS problem on G seeks to identify an IS with maximum cardinality, i.e. MIS. $V^{*}\subseteq V$ is a vertex cover, i.e. VC of $G=(V,E)$ if every $e\in E$ is incident upon at least one vertex in $V^{*}$ . $V^{*}\subseteq V$ is dominating set, DS, of $G=(V,E)$ if $\forall v\in V$ either $v\in V^{*}$ or $\exists u\in V^{*}$ and $(u,v)\in E$ . The MVC problem on G seeks to identify a vertex cover with minimum cardinality, i.e. MVC. Likewise, MCV seeks a connected vertex cover, i.e. VC which forms one component in G, with minimum cardinality, i.e. MCV. A connected DS, CDS, is a DS that forms a connected component in G. The problems MDS and MCD seek to identify a DS and a connected DS i.e. CDS respectively with minimum cardinalities. MIS, MVC, MDS, MCV and MCD on a general graph are known to be NP-complete. Polynomial time algorithms are known for bipartite graphs, chordal graphs, cycle graphs, comparability graphs, claw-free graphs, interval graphs and circular arc graphs for some of these problems. We introduce a novel graph class, layered graph, where each layer refers to a subgraph containing at most some k vertices. Inter layer edges are restricted to the vertices in adjacent layers. We show that if $k=\Theta (log\mid V\mid )$ then MIS, MVC and MDS can be computed in polynomial time and if $k=O\left({(log\mid V\mid )}^{\alpha }\right)$ , where $\alpha <1$ , then MCV and MCD can be computed in polynomial time. If $k=\Theta \left({(log\mid V\mid )}^{1+ϵ}\right)$ , for $ϵ>0$ , then MIS, MVC and MDS require quasi-polynomial time. If $k=\Theta (log\mid V\mid )$ then MCV, MCD require quasi-polynomial time. Layered graphs do have constraints such as bipartiteness, planarity and acyclicity.

This paper aims to demonstrate how the detection of communities for epidemiological networks can be used to reduce the dimensionality of optimal control problems. A range of planted partition stochastic blockmodels are generated, where the underlying communities in the model vary in detectability. A high-performance computing workflow is then used to simulate and estimate reduced model dynamics for the networks, which in turn is used to find optimal control solutions.

Shale gas and coalbed methane are energy sources that partly or even mainly consist of methane stored in an adsorbed state in the pores of the organic-rich rock and coal seams. In this study, the graphene nanoslit model is employed to model the nanometer slit pores in shale and coal. The Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) modelling methods are used to investigate the mechanisms of adsorption and displacement of methane in the slit pore. It was found that the CVFF forcefield resulted in the largest adsorption amount, while the PCFF forcefield resulted in the least. The COMPASS and COMPASS II force fields led to similar results. As the width of the slit pore increases, the adsorption amount of gas molecules increases, and the number density profile of adsorbed methane molecules alters from a single adsorption layer to multi-adsorption layers. The minimum slit pore width at which methane molecules can penetrate the slit pore was found to be 0.7 nm. Moreover, it is demonstrated that by lowering the temperature, the adsorption rate of the methane increases since the adsorption is an exothermic process. Enhancing methane recovery was investigated by the injection of gases such as CO2 and N2 to displace the adsorbed methane. The comparison of adsorption isotherms of gas molecules provided the following order in terms of the amount of adsorption: CO2&gt;CH4&gt;N2.

The application of the sensors optical fiber in the areas of scientific instrumentation and industrial instrumentation is very attractive due to its numerous advantages. In the industry of civil engineering for example, quasi-distributed sensors made with optical fiber are used for reliable strain and temperature measurements. Here, a quasi-distributed sensor in the frequency domain is discussed. The sensor consists of a series of low-finesse Fabry-Perot interferometers where each Fabry-Perot interferometer acts as a local sensor. Fabry-Perot interferometers are formed by pairs of identical low reflective Bragg gratings imprinted in a single mode fiber. All interferometer sensors have different cavity length, provoking the frequency-domain multiplexing. The optical signal represents the superposition of all interference patterns which can be decomposed using the Fourier transform. The frequency spectrum is analyzed and sensor´s properties were defined. Following, a quasi-distributed sensor was numerically simulated. Our sensor simulation considers sensor properties, signal processing, noise system and instrumentation. The numerical results show the behavior of resolution vs. signal-to-noise ratio. From our results, the Fabry-Perot sensor has high resolution and low resolutions. Both resolutions are conceivable because the FDPA algorithm elaborates two evaluations of Bragg wavelength shift

Integrating photovoltaic (PV) systems plays a pivotal role in the global shift toward renewable energy, offering significant environmental benefits. In many instances, utility companies pay for both energy and peak demand. PV systems can effectively reduce energy consumption during their operational periods, while, their contribution during peak hours, particularly in the early morning and late evening, is often limited. To address this issue, battery storage systems are utilized for storing energy during off-peak hours and releasing it during peak times. However, finding the optimal size of PV-battery systems remains a challenge. This study proposes a novel statistical model to optimize PV-battery system sizing for enhanced peak demand reduction. The model aims to flatten 95% of peak demand days up to a certain demand threshold, ensuring consistent energy supply and financial benefit for utility companies. A straightforward and effective search methodology is employed to determine the optimal system sizes. Additionally, the model effectiveness is rigorously tested through a modified Monte Carlo simulation coupled with time series clustering to generate various scenarios to assess performance under different conditions. The results indicate that the optimal PV-battery system successfully flattens 95% of daily peak demand up to 2000 kW, yielding a financial benefit of $812,648 over 20 years.

In this work we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three dimensional space defined by energy, temperature and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of large fermion systems with numerical simulation. As an example, we present energies for up to 50 noninteracting fermions and up to 20 interacting fermions at all temperatures and show good agreement with the analytical result for noninteracting case.

Electricity distribution network operators face several operational constraints in the provision of safe and reliable power given that investments for network area reinforcement must be commensurate with improvements in network reliability. This paper provides an integrated approach for assessing the impact of different operational constraints on distribution network reliability by incorporating component lifetime models, time-varying component failure rates as well as the monetary cost of customer interruptions in an all-inclusive probabilistic methodology that applies a time-sequential Monte Carlo simulation. A test distribution network based on the Roy Billinton test system is modelled to investigate the system performance when overloading limits are exceeded as well as when preventive maintenance is performed. Standard reliability indices measuring the frequency and duration of interruptions, and the energy not supplied are complemented with a novel monetary reliability index. The comprehensive assessment includes not only average indices but also their probability distributions to adequately describe the risk of customer interruptions. Results demonstrate the effectiveness of this holistic approach as the impacts of operational decisions are assessed from both reliability and monetary perspectives. This informs network planning decisions through optimum investments and consideration of customer outage costs.

In this paper we provide some important and significant observations on invo-regular rings appeared in Ann. Univ. Mariae Curie-Sklodowska Sect. A Mathematica (2018). In addition we exhibit that strongly invo-regular rings, invo-regular rings and quasi invo-regular rings all coincide with the well known and well characterized tripotent rings and all these rings further coincide with the well known unit regular rings.

Rubble stone masonry walls are widely diffused in most of the cultural and architectural heritage of historical cities. The mechanical response of such material is rather complicated to predict due to their composite nature. Vertical compression tests, diagonal compression tests, and shear-compression tests are usually adopted to experimentally investigate the mechanical properties of stone masonries. However, further tests are needed for the safety assessment of these ancient structures. Since the relation between normal and shear stresses plays a major role in the shear behavior of masonry joints, governing the failure mode, triplet test configuration was here investigated. First, the experimental tests carried out at the laboratory (LPMS) of the University of L'Aquila on stone masonry specimens were presented. Then, the triplet test was simulated by using the Total Strain Crack Model, which reflects all the ultimate states of quasi-brittle material such as cracking, crushing and shear failure. The goal of the numerical investigation was to evaluate the shear mechanical parameters of the masonry sample, including strength, dilatancy, normal and shear deformations. Furthermore, the effect of (i) confinement pressure and (ii) bond behavior at the sample-plates interfaces were investigated, showing that they can strongly influence the mechanical response of the walls.

Statistical modeling of lifetime data plays an essential role in a wide range of practical fields, such as health and engineering. There have been a lot of studies done to develop statistical models that can better describe health data than traditional models. For the first time, we pioneer a novel family of continuous probability distributions called the generalized odd beta prime generalized (GOBP-G) family of distributions. The cumulative distribution and probability density functions of the new family are presented. A new generalization of the Weibull distribution called "generalized odd beta prime-Weibull" (GOBPW) is proposed using the pioneered GOBP-G family. The mixture representations of the new distribution are defined and derived. Some formal statistical properties of the GOBPW distribution, such as the moments, moment generating function, incomplete moments, information generating function, entropies, stress-strength function, quantile function, and order statistics, are derived. The estimation of the parameters of the proposed distribution is evaluated using the maximum likelihood estimation approach. Different cancer disease data sets, such as the bladder, head and neck, acute bone, and blood cancers, are used to illustrate the applicability and usefulness of the new model and were compared using several statistical accuracy measures with that of well-established extended Weibull distributions, which are the beta modified Weibull distribution, Kumaraswamy modified Weibull distribution, gamma generalized modified Weibull distribution, gamma log-logistic Weibull distribution, and beta log-logistic Weibull distribution. The results show that the proposed model gives better results than the competitive models. This study could guide the relevant stakeholders in choosing a suitable statistical model for the health data instead of relying on traditional models to enhance decision-making.

The requirements for efficient replacement of fossil fuel, combined with the growing energy crisis, put hydrogen production into focus. Efficient and cost-effective electrocatalysts are needed for H2 production, and novel strategies for their discovery must be developed. Here we utilized Kinetic Monte Carlo (KMC) simulations to demonstrate that hydrogen evolution reaction (HER) can be boosted via hydrogen spillover to the support when the catalyst surface is largely covered by adsorbed hydrogen under operating conditions. Based on the insights from KMC, we synthesized a series of reduced graphene oxide-supported catalysts and compared their activities towards HER in alkaline media with that of corresponding pure metals. For Ag, Au, and Zn, the support effect is negative, but for Pt, Pd, Fe, Co, and Ni, the presence of the support enhances HER activity. The HER volcano, constructed using calculated hydrogen binding energies and measured HER activities, shows a positive shift of the strong binding branch. This work demonstrates the possibilities of metal|support interface engineering for producing effective HER catalysts and provides general guidelines for choosing novel catalyst-support combinations for electrocatalytic hydrogen production.

We outline and test a new methodology for genuine simulation of stochastic processes with any dependence and any marginal distribution. We reproduce time dependence with a generalized, time symmetric or asymmetric, moving-average scheme. This implements linear filtering of non-Gaussian white noise, with the weights of the filter determined by analytical equations in terms of the autocovariance of the process. We approximate the marginal distribution of the process, irrespective of its type, using a number of its cumulants, which in turn determine the cumulants of white noise in a manner that can readily support the generation of random numbers from that approximation, so that it be applicable for stochastic simulation. The simulation method is genuine as it uses the process of interest directly without any transformation (e.g. normalization). We illustrate the method in a number of synthetic and real-world applications with either persistence or antipersistence, and with non-Gaussian marginal distributions that are bounded, thus making the problem more demanding. These include distributions bounded from both sides, such as uniform, and bounded form below, such as exponential and Pareto, possibly having a discontinuity at the origin (intermittence). All examples studied show the satisfactory performance of the method.

The planning and decision-making for a distributed energy supply concept in complex actor structures like in districts calls for the approach to be highly structured. Here, a strategy with strong use of energetic simulations is developed, the core elements are presented and research gaps are identified. The exemplary implementation is shown using the case study of a new district on the former Oldenburg airbase in northwestern Germany. The process is divided into four consecutive phases, which are carried out with different stakeholder participation and use of different simulation tools. Based on a common objective, a superstructure of the applicable technologies is developed. Detailed planning is then carried out with the help of a multi-objective optimal sizing algorithm and Monte Carlo based risk assessment. The process ends with the operating phase, which is to guarantee a further optimal and dynamic mode of operation. The main objective of this publication is present the core elements of the planning processes and decision-making framework based on the case study and to find and identify research gaps that will have to be addressed in the future.

The present work deals with the development of finite element methodology for obtaining the stress distributions in thick cylindrical HK40 stainless steel pipe that carry high temperature fluids. The material properties and loading are assumed to be random variables. Thermal stresses that are generated along radial, axial and tangential directions are computed generally using analytical expressions which are very complex. To circumvent such an issue, the probability theory and mathematical statistics have been applied to many engineering problems which allows to determine the safety both quantitatively and objectively based on the concepts of reliability. Monte Carlo simulation methodology is used to study the probabilistic characteristics of thermal stresses which is used for estimating the probabilistic distributions of stresses against the variations arising due to material properties and load. A 2-D Probabilistic finite element code is developed in MATLAB and the deterministic solution is compared with ABAQUS solutions. The values of stresses that are obtained from the variation of elastic modulus are found to be low as compared to the case where the load alone is varying. The probability of failure of the pipe structure is predicted against the variations in internal pressure and thermal gradient. These finite element framework developments are useful for the life estimation of piping structures in high temperature applications and subsequently quantifying the uncertainties in loading and material properties.

Previous studies generally reckon that the full converter-base wind power generation (FCWG) is a ’decoupled’ power source from the grid, which hardly participates in electromechanical oscillations. However, it is found recently that strong interaction could be induced which might incur severe resonance incidents in electromechanical dynamic timescale. In this paper, the participation of FCWG in electromechanical dynamics is extensively investigated, and particularly, an unusual transition of electromechanical oscillation mode (EOM) is uncovered for the first time. The detailed mathematical models of open-loop and closed-loop power systems are firstly established, and modal analysis is employed to quantify the FCWG participation in electromechanical dynamics, with two new mode identification criteria, i.e., FCWG dynamics correlation ratio (FDCR) and quasi-electromechanical loop correlation ratio (QELCR). On this basis, the impact of different wind penetration levels and controller parameter settings on the participation of FCWG is investigated. It is revealed that if an FOM has a similar oscillation frequency to the system EOMs, there is a high possibility to induce strong interactions between FCWG dynamics and system electromechanical dynamics of the external power systems. In this circumstance, an interesting phenomenon may occur that an EOM may be dominated by FCWG dynamics, and hence is transformed into a quasi-EOM, which actively involves the participation of FCWG quasi-electromechanical state variables.

In this paper, we study the effect of decreasing aperture diameter on signal transmission from the scintillator to the PMT. The apertures have been inserted directly between the scintillator and the photomultiplier. Three different aperture diameters have been used over a wide energy range.A newly developed fast digital spectrometer has been utilised for experimental measurement in the radial channel of the VR-1 research reactor. The detector signals with and without inserted apertures have been measured and evaluated. The aim of our research was to determine the effect of the apertures on the spectral quality, energy range and the evaluated physical quantities.

In this paper, we will use the maximum likelihood estimation (MLE) and the Bayes methods to perform estimation procedures for the reliability of stress-strength ℜ=P[X>Y] based on independent adaptive progressive censored samples that were taken from the Chen distributions. An approximate confidence interval of ℜ is constructed using a variety of classical techniques, such as the normal approximation of the MLE, the normal approximation of the log-transformed MLE, and the percentile bootstrap (Boot-p) procedure. The Bayesian estimation of ℜ is obtained based on the balanced loss function, which comes in two versions here, the symmetric balanced squared error (BSE) loss function and the asymmetric balanced linear exponential (BLINEX) loss function. When estimating ℜ using the Bayesian approach, all the unknown parameters of the Chen distributions are assumed to be independently distributed and to have informative gamma priors. Further, a mixture of Gibbs sampling algorithm and Metropolis-Hastings algorithm is used to compute the Bayes estimate of ℜ and the associated highest posterior density credible interval. In the end, simulation research are used to assess the general overall performance of the proposed estimators and a real data set is provided to exemplify the theoretical results.

The main issue in this work is to study limit functions necessary for the reliability assessment of the structural steel with the use of the relative entropy apparatus. This would be done using a few different mathematical theories relevant to this relative entropy, namely those proposed by Bhattacharyya, Kullback-Leibler, Jeffreys, and Hellinger. Probabilistic analysis in the presence of uncertainty in material characteristics would be delivered using three different numerical strategies – Monte-Carlo simulation, stochastic perturbation method as well as the semi-analytical approach. All these methods are based on the Weighted Least Squares Method approximations of the structural response functions versus the given uncertainty source and they would allow efficient determination of the first two probabilistic moments of the structural responses including stresses, displacements, and strains. The entire computational implementation would be delivered using the Finite Element Method system ABAQUS and computer algebra program MAPLE, where relative entropies, as well as polynomial response functions, would be determined. This study would demonstrate that the relative entropies may be efficiently used in reliability assessment close to the widely engaged First Order Reliability Method (FORM).

Bayesian nonlinear mixed effects models for data in the form of continuous, repeated measurements from a population, also known as Bayesian hierarchical nonlinear models, are a popular platform for analysis when interest focuses on individual specific characteristics and relevant uncertainty quantification. Due to the limitation of computational power, this framework was relatively dormant until the late 1980s, but in recent years, the statistical research community saw vigorous development of new methodological and computational techniques for these models, the emergence of software, and wide application of the models in numerous industrial and academic fields. This article presents an overview of the formulation, interpretation, and implementation of Bayesian nonlinear mixed effects models and surveys recent advances and applications.

In this paper, quasi-static axial compression tests are performed on the nodal Moso bamboos to study the size effect on energy absorption of the bamboos and the damage pattern of the multiple bamboo columns. Experimental results show that under the same moisture content, growth age and growing environment, the specific energy absorption (SEA) of the test samples increases with the increase of the out-diameter and thickness of the bamboo columns, indicating that size effect exists for energy absorption of the Moso Bamboos. For the multiple bamboo columns, there are mainly three failure modes for the constituent single bamboo columns: splitting above the node, splitting below the node and splitting through the node. Also, the tensile tests are conducted on three kinds of dog-bone shaped bamboo samples to investigate the macroscopic tensile fracture mode in the longitudinal direction of Moso bamboos. Results show that there is no direct relationship between the fracture pattern and moisture content of the bamboos, as well as the growth age of the bamboos. However, the tensile loading rate and the shape of the dog-bone shaped bamboo sample could affect the macroscopic fracture pattern of the bamboos in some cases.

Demographic change is leading to the aging of German society. As long as the baby boom co-horts are still of working age, the working population will also age - and decline as soon as this baby boom generation gradually reaches retirement age. At the same time, there has been a trend towards increasing absenteeism (times of inability to work) in companies since the zero years, with the number of days of absence increasing with age. We present a novel stochastic forecast approach that combines population forecasting with forecasts of labor force participation trends, considering epidemiological aspects. For this, we combine a stochastic Monte Carlo-based cohort-component forecast of the population with projections of labor force participation rates and morbidity rates. This article examines the purely demographic effect on the economic costs associated with such absenteeism due to the inability to work. Under expected future employment patterns and constant morbidity patterns, absenteeism is expected by close to 5 percent by 2050 relative to 2020, associated with increasing economic costs of almost 3 percent. Our results illustrate how strongly the pronounced baby boom/ baby bust phenomenon determines demographic development in Germany in the midterm.

This study addresses the central question that how can a precise measurement of the gluon Parton Distribution Function (PDF) be attained, thereby enhancing our understanding of the intricate quark-gluon composition within nucleons? Through a dual approach involving Quantum Chromodynamics (QCD) and high-energy particle collisions, the research combines computational tools and experimental techniques, focusing on protons and neutrons. Scattering experiments and advanced detectors are utilized to probe nucleon structure. Theoretical frameworks like QCD equations and the Momentum Sum Rule investigate Parton Distribution Functions at specific momentum fractions and energy scales. The primary goal is to achieve an accurate measurement of the gluon PDF, crucial for unraveling strong force complexities and internal nucleon dynamics. The study's theoretical plots and equations illustrate hypothetical behaviors of gluon distribution, quark distribution, and splitting functions, contributing to a deeper understanding of fundamental particles.

Autonomous shuttles have been used as end mile solutions for smart mobility in smart cities. The urban driving conditions of smart cities with many other actors sharing the road and the presence of intersections have posed challenges to the use of autonomous shuttles. Round intersections are more challenging as it is more difficult to perceive the other vehicles in and near the intersection. Thus, this paper focuses on decision making of autonomous vehicles for handling round intersections. The round intersection is introduced first, followed by introductions of the Markov Decision Process (MDP), the Partially Observable Markov Decision Process (POMDP) and the Object Oriented Partially Observable Markov Decision Process (OOPOMDP) which are used for decision making with uncertain knowledge of the motion of the other vehicles. The Partially Observable Monte-Carlo Planning (POMCP) algorithm is used as the solution method and OOPOMDP is applied to decision making for autonomous vehicles in round intersections. Decision making is formulated first as a POMDP problem, and the penalty function is formulated and set accordingly. This is followed by improvement of decision making with policy prediction. Augmented objective state and policy-based state transition are introduced and simulations are used to demonstrate the effectiveness of the proposed method for collision free handling of round intersections by the ego vehicle.

The goal of the manuscript is to design a relatively best control structure for the noise suppression of a drone’s camera gimbal action. The gimbal’s movement can be simplified as a rest-to-rest reorientation system that can achieve the boundary result of a dynamic system. Six different control architectures are proposed and evaluated based on their ability to control the trajectory of the dynamic-system position and speed, their running time, and quadratic cost. The robustness of the control architecture to uncertainties in inertia and sensor noise is also analyzed. Monte Carlo figures are used to assess the performance of the six control systems. The conditions for applying different architectures are identified through this analysis. The analysis and experimental tests reveal the most suitable control of the drone’s camera gimbal rotation.

In this review, we trace the evolution of the quantum spin-wave theory treating non-collinear spin configurations. Non-collinear spin configurations are consequences of the frustration created by competing interactions. They include simple chiral magnets due to competing nearest-neighbor (NN) and next-NN interactions and systems with geometry frustration such as the triangular antiferromagnet and the Kagomé lattice. We review here spin-wave results of such systems and also systems with the Dzyaloshinskii-Moriya interaction. Accent is put on these non-collinear ground states which have to be calculated before applying any spin-wave theory to determine the spectrum of the elementary excitations from the ground states. We mostly show results from a self-consistent Green’function theory to calculate the spin-wave spectrum and the layer magnetizations at finite T in two and three dimensions as well as in thin films with surface effects. Some new unpublished results are also included. Analytical details and the validity of the method are shown and discussed.

The paper examines common elements between Lévai’s and Milson’s potentials obtained by Liouville transformations of two rational canonical Sturm-Liouville equations (RCSLEs) with even density functions which are exactly solvable via Jacobi polynomials in a real or accordingly imaginary argument. We refer to the polynomial numerators of the given rational density function as ‘tangent polynomial’ (TP) and thereby term the aforementioned potentials as ‘e-TP’ A special attention is given to the overlap between the two potentials along symmetric curves which represent two different forms of the Ginocchio potential exactly quantized via Gegenbauer and Masjed-Jamei polynomials respectively. Our analysis reveals that the actual interconnection between Lévai’s parameters for these two rational realizations of the Ginocchio potential is much more complicated than one could expect based on the striking resemblance between two quartic equations derived by Lévai for ‘averaged’ Jacobi indexes.

For the in-house validation of a droplet digital PCR method, a factorial experimental design was implemented. This design serves different purposes. On the one hand, it is an efficient design in relation to the workload involved in achieving a desirable level of reliability of variance estimates. On the other hand, it allows a partitioning of total variance into different components, thus providing information regarding the dominant sources of random variation. The statistical modelling reflects the actual measurement mechanism, establishing relationships between nominal target DNA copies per well, the range of variation of copy numbers per droplet, probability of detection values, and estimated numbers of copies.

A quantitative analysis of X-ray emission from an electron cyclotron resonance (ECR) plasma was performed to probe the spatial properties of electrons having energy for effective ionisation. A series of measurements were taken by INFN-LNS and ATOMKI, capturing spatially and spectrally resolved X-ray maps as well as volumetric emissions from argon plasma. Comparing the former with model generated maps (involving space-resolved phenomenological electron energy distribution function and geometrical efficiency calculated using ray-tracing MC routine) furnished information on structural aspects of the plasma. Similarly, fitting a model composed of bremsstrahlung and fluorescence to the volumetric X-ray spectrum provided valuable insight into the density and temperature of confined and lost electrons. The latter can be fed back to existing electron kinetics models for simulating more relevant energies, consequently improving theoretical X-ray maps and establishing the method as an excellent indirect diagnostic tool for warm electrons, required for both fundamental and applied research in ECR plasmas.

The classical Rotational Isomeric State (RIS) model, originally proposed by Flory, has been used to rationalize a wide range of physicochemical properties of neutral polymers. However, many weak polyelectrolytes of interest are able to regulate their charge depending on the conformational state of the bonds. Recently, it has been shown that the RIS model can be coupled with the Site Binding (SB) model, for which the ionizable sites can adopt two states: protonated or deprotonated. The resulting combined scheme, the SBRIS model, allows to analyse ionization and conformational equilibria on the same foot. In the present work this approach is extended to include pH-dependent electrostatic Long Range (LR) interactions, ubiquitous in weak polyelectrolytes at moderate and low ionic strengths. With this aim the original LR interactions are taken into account by defining effective Short Range (SR) and pH-dependent parameters, such as effective microscopic protonation constants and rotational bond energies. The new parameters are systematically calculated using variational methods. The machinery of statistical mechanics for SR interactions, including the powerful and fast transfer matrix methods, can then be applied. The resulting technique, to which we will refer as Local Effective Interaction Parameters (LEIP) method, is illustrated with a minimal model of a flexible linear polyelectrolyte containing only one type of rotating bonds. LEIP reproduces very well the pH dependence of the degree of protonation and bond probabilities obtained by semi-grand canonical Monte Carlo simulations, where LR interactions are taken explicitly into account. The reduction in the computational time in several orders of magnitude suggests that the LEIP technique could be useful in a range of areas involving linear weak polyelectrolytes, allowing direct fitting of the relevant physical parameters to the experimental quantities.

Keywords: liquid alloys, microstratification, dispersed particles, spontaneous dispersion, criteria for true (entropic) spontaneous dispersion, mechanisms of quasi-spontaneous dispersion, methods and results for assessing size of dispersed particles.

In this study, we present the concept of interval-valued fuzzy soft point and then introduce the notions of neighborhood and quasi-neighbourhood of it in interval-valued fuzzy soft topological spaces. Separation axioms in interval-valued fuzzy soft topology, so-called $q$-$T_{i}$ for $ i=0,1,2,3,4 ,$ is introduced and some of its basic properties are also studied.

In this paper a new double-nozzle technique for in-gas-jet laser resonance ionization spectroscopy is proposed. The functionality of this new technique we explored by detailed gas dynamic and Monte Carlo atom-trajectory simulations, which results are presented and discussed. The results of similar computer simulations for JetRIS setup (as a typical representative of the conventional nowadays in-gas-jet technique) are presented and discussed as well. The direct comparison of calculation results for proposed new technique with the conventional one shows that the double-nozzle technique has many advantages compared with that one used in the JetRIS setup at GSI for future high-resolution laser spectroscopic study of heaviest elements. To fully implement the proposed new technique to all existing (or under construction) setups for in-gas-jet laser resonance ionization spectroscopy, it will be enough just replace the used their supersonic nozzle by the miniature double-nozzle device described in the paper.

This work addresses the role of charge regulation (CR) and the associated fluctuations in the conformational and mechanical properties of weak polyelectrolytes. Due to CR, changes in the pH-value modifies the average macromolecular charge and conformational equilibria. A second effect is that, for a given average charge per site, fluctuations can alter the intensity of the interactions by means of correlation between binding sites. We investigate both effects by means of Monte Carlo simulations at constant pH-value, so that the charge is a fluctuating quantity. Once the average charge per site is available, we turn off the fluctuations by assigning the same average charge to every site. A constant charge MC simulation is then performed. We make use of a model which accounts for the main fundamental aspects of a linear flexible polyelectrolyte i.e. proton binding, angle internal rotation, bond stretching and bending. Steric excluded volume and differentiated treatment for short-range and long-range interactions are also included in the model. It can be regarded as a kind of "minimal'' model in the sense that contains a minimum number of parameters but still preserving the atomistic detail. It is shown that, if fluctuations are activated, gauche state bond probabilities increase, and the persistence length decreases, so that the polymer becomes more folded. Macromolecular stretching is also analyzed in presence of CR (the charge depends on the applied force) and without CR (the charge is fixed to the value at zero force). The analysis of the low force scaling behavior concludes that Pincus exponent becomes pH-dependent. Both with and without CR, a transition from 1/2 at high pH-values (phantom chain) to 3/5 to low pH-values (Pincus regime), is observed. Finally, the intermediate force stretching regime is investigated. It is found that CR induces a moderate influence in the force-extension curves and persistence length (which in this force regime becomes force-dependent). It is thus concluded that the effect of CR on the stretching curves is mainly due to changes in the average charge at zero force. It is also found that, for the cases studied, the effect of steric excluded volume is almost irrelevant compared to electrostatic interactions.

Rapid developments in the field of underwater photogrammetry have given scientists1the ability to produce accurate 3-dimensional (3D) models which are now increasingly used in the representation and study of local areas of interest. This paper addresses the lack of systematic analysis of 3D reconstruction and navigation fusion strategies, as well as associated error evaluation of models produced at larger scales in GPS-denied environments using a monocular camera (often in deep-sea scenarios). Based on our prior work on automatic scale estimation of Structure from Motion (SfM)-based 3D models using laser scalers, an automatic scale accuracy framework is presented. The confidence level for each of the scale error estimates is independently assessed through the propagation of the uncertainties associated with image features and laser spot detections using a Monte Carlo simulation. The number of iterations used in the simulation was validated through the analysis of the final estimate behaviour. To facilitate the detection and uncertainty estimation of even greatly attenuated laser beams, an automatic laser spot detection method, mitigating the effects of scene texture, was developed, with the main novelty of estimating the uncertainties based on the recovered characteristic shapes of laser spots with radially decreasing1 intensities. The effects of four different reconstruction strategies resulting from the combinations of Incremental/GlobalSfM, and thea priori/a posterioriuse of navigation data were analyzed using two distinct survey scenarios captured during the SUBSAINTES 2017 cruise (doi: 10.17600/17001000). The study demonstrates that surveys with multiple overlaps of non-sequential images result in a nearly identical solution regardless of the strategy (SfM or navigation fusion), while surveys with weakly connected sequentially acquired images are prone to produce broad-scale deformation (doming effect) when navigation is not included in the optimization. Thus the scenarios with complex survey patterns substantially benefit from using multi-objective BA navigation fusion. In all cases, the errors in the models are inferior to 5%, with errors often being around 1%. The effects of combining data from multiple surveys were also evaluated. The introduction of additional vectors in the optimization of multi-survey problems successfully accounted for offset changes present in the underwater USBL-based navigation data and thus minimize the effect of contradicting navigation priors. Our results also illustrate the importance of collecting a multitude of evaluation data at different locations and moments during the survey.

This article presents a comprehensive study on enhancing grid resilience through advanced forecasting and optimization techniques in the context of power outages. Power outages pose significant challenges to modern societies, affecting various sectors such as industries, households, and critical infrastructures. The research combines statistical analysis, machine learning algorithms, and optimization methods to address this issue to develop a holistic approach for predicting and mitigating power outage events. The proposed methodology involves the use of Monte Carlo simulations in MATLAB for future outage prediction, Long Short-Term Memory (LSTM) networks for forecasting solar irradiance and load profiles, and a hybrid LSTM-Particle Swarm Optimization (PSO) model to improve accuracy. Furthermore, the role of Battery State of Charge (SoC) in enhancing system resilience is explored. The study also assesses the techno-economic advantages of a grid-tied microgrid integrated with solar panels and batteries over conventional grid systems. The results highlight the potential of the proposed approach in strengthening grid resilience, reducing downtime, and fostering sustainable energy utilization.

The aim of this study is to obtain the Bayes estimators and the maximum likelihood estimators (MLEs) for the unknown parameters of the Rayleigh Weibull (RW) distribution based on progres-sive type-II censored samples. The approximate Bayes estimators are calculated using the idea of Lindley and Tierney-Kadane's approximation method under the squared-error loss function when the Bayes estimators are not handed in explicit forms. In this study, the approximate Bayes esti-mates are compared with the maximum likelihood estimates in the aspect of the estimated risks (ERs) using Monte Carlo simulation. In addition, the coverage probabilities of the parametric bootstrap estimates are calculated. Real lifetime data sets belonging to the cancer types as bladder cancer, head and neck cancer, and leukemia are used to illustrate the emprical results belonging to the approximate Bayes estimates, the maximum likelihood estimates, and the parametric bootstrap intervals.

This paper is a further development of our “Proposal of a new double-nozzle technique for in-gas-jet laser resonance ionization spectroscopy” published in journal Atoms earlier in this year. Here, we propose to equip the proposed double-nozzle technique with the RF-only funnel and RF-buncher placed in the gas-jet chamber at 70 mm distance downstream the double-nozzle exit. It allows for highly effective extraction into vacuum heavy ion beams, produced in two-steps laser resonance ionization in the argon supersonic jet. We explored the operation of this new full version of the double-nozzle technique through detailed gas dynamic and then Monte Carlo trajectory simulations, which results are presented and discussed. In particularly, we have shown by calculations that more than 80% of all nobelium-254 neutral atoms, extracted by argon flow from the gas stopping cell through the double-nozzle into the supersonic jet can be extracted into vacuum in a form of pulsed ion beam having low transverse and longitudinal emittance.s dynamic and Monte Carlo simulations

Creating self-learning algorithms, developing deep neural networks and improving other methods that "learn" for various areas of human activity is the main goal of the theory of machine learning. It helps to replace the human with a machine, aiming to increase the quality of production. The theory of artificial neural networks, which already have replaced the humans in problems of detection of moving objects, recognition of images or sounds, time series prediction, big data analysis and numerical methods remains the most dispersed branch of the theory of machine learning. Certainly, for each area of human activity it is necessary to select appropriate neural network architectures, methods of data processing and some novel tools from applied mathematics. But the universal problem for all these neural networks with specific data is the achieving the highest accuracy in short time. Such problem can be resolved by increasing sizes of architectures and improving data preprocessing, where the accuracy rises with the training time. But there is a possibility to increase the accuracy without time growing, applying optimization methods. In this survey we demonstrate existing optimization algorithms of all types, which can be used in neural networks. There are presented modifications of basic optimization algorithms, such as stochastic gradient descent, adaptive moment estimation, Newton and quasi-Newton optimization methods. But the most recent optimization algorithms are related to information geometry, for Fisher-Rao and Bregman metrics. This approach in optimization extended the theory of classic neural networks to quantum and complex-valued neural networks, due to geometric and probabilistic tools. There are provided applications of all introduced optimization algorithms, what delighted many kinds of neural networks, which can be improved by including any advanced approaches in minimization of the loss function. Afterwards, we demonstrated ways of developing optimization algorithms in further researches, engaging neural networks with progressive architectures. Classical gradient based optimizers can be replaced by fractional order, bilevel and, even, gradient free optimization methods. There is a possibility to add such analogues in graph, spiking, complex-valued, quantum and wavelet neural networks. Besides the usual problems of image recognition, time series prediction, object detection, there are many are other tasks for modern theory of machine learning, such as solving problem of quantum computations, partial differential and integro-differential equations, stochastic processes and Brownian motion, making decisions and computer algebra.

The local and global order in dense packings of linear, semi-flexible polymers of tangent hard spheres are studied by employing extensive Monte Carlo simulations at increasing volume fractions. Chain stiffness is controlled by a tunable harmonic potential for the bending angle whose intensity dictates the rigidity of the polymer backbone as a function of the bending constant and equilibrium angle. The studied angles range from acute to obtuse ones, reaching the limit of rod-like polymers. We analyze how packing density and chain stiffness affect the ability of chains to self-organize at the local and global levels. The former corresponds to crystallinity as quantified by the Characteristic Crystallographic Element (CCE) norm descriptor, while the latter is computed through the scalar orientational order parameter. In all cases, we identify the critical volume fraction for the phase transition and gauge the established crystal morphologies, developing a complete phase diagram as a function of packing density and equilibrium bending angle. A plethora of structures is obtained, ranging from random hexagonal closed packed morphologies of mixed character and almost perfect face centered cubic (FCC) and hexagonal close-packed (HCP) crystals at the level of monomers, to nematic mesophases, with prolate and oblate mesogens at the level of chains. For rod-like chains, hysteresis is observed between the establishment of long-range nematic order and crystallization, while for right-angle chains both transitions are synchronized. A comparison is also provided against the analogous packings of monomeric and fully flexible chains of hard spheres.

In the present work we revise and extend the Characteristic Crystallographic Element (CCE) norm, an algorithm used to simultaneously detect radial and orientational similarity of computer-generated structures with respect to specific reference crystals and local symmetries. Based on the identification of point group symmetry elements, the CCE descriptor is able to gauge local structure with high precision and finely distinguish between competing morphologies. As test cases we use computer-generated monomeric and polymer systems of spherical particles interacting with the hard-sphere and square-well attractive potentials. We demonstrate that the CCE norm is able to detect and differentiate, between others, among: hexagonal close packed (HCP), face centered cubic (FCC), hexagonal (HEX) and body centered cubic (BCC) crystals as well as non-crystallographic fivefold (FIV) local symmetry in bulk 3-D systems; triangular (TRI), square (SQU) and honeycomb (HON) crystals, as well as pentagonal (PEN) local symmetry in thin films of one-layer thickness (2-D systems). The descriptor is general and can be applied to identify the symmetry elements of any point group for arbitrary atomic or particulate system in two or three dimensions, in the bulk or under confinement.