Managing a portfolio whose value closely tracks an index by trading only in a subset of the index constituents involves an NP-hard optimisation problem. In the prior literature, it has been suggested to solve this problem using sequential Monte Carlo (SMC, a.k.a. particle filter) methods. However, this literature does not take transaction costs into account, although transaction costs are the primary motivation for attempting to replicate the index by trading in a subset, rather than the full set of index constituents. This paper extends the SMC approach to index tracking to allow for proportional transaction costs and implements this extended method on empirical data for a variety stock indices. In addition to providing a more practically useful tracking strategy by allowing for transaction costs, we find that including a penalty for transaction costs in the optimisation objective can actually lead to better tracking performance.

Sword-duels are interesting, from a game-theoretical perspective, as they correspond to sudden-death interactions between players (i.e., defeat means removal from the game), and where victory/defeat depend on not just skill, but also luck. Analyzing probabilities of victory streaks, given a certain level of self and others' skill, is thus relevant for any application domain that can be modeled in the same way. This paper takes inspiration from Miyamoto Musashi, famously undefeated for 61 duels, and implements a Markov-chain Monte-Carlo simulation approach to evaluate this scenario. Results suggest that a 61 victory streak can be probabilistically observed when skill level is roughly 6.5 times that of the average duelist, and with 95\% confidence when skill level is roughly 1000 times that of the average duelist. More generally, this paper provides a method for determining chances of victory streaks in game-theoretical sudden-death encounters, when both "skill" and "luck" contribute to the outcome of the encounter. Specific scenarios can be modeled by modifying the utilized Markov chain and adjusting sampled distributions as required.

This paper compares the finite sample performance of three non-parametric threshold estimators via Monte Carlo method. Our results show that the finite sample performance of the three estimators is not robust to the relative position of the threshold level along the distribution of threshold variable, especially when a structural change occurs at the tail part of the distribution.

This paper concentrates on computing the efficiency of heliostat field systems. By conducting a comprehensive analysis of the mathematical model for optical efficiency, we have effectively simplified the complex three-dimensional illumination problem of heliostat fields into a two-dimensional issue. Our work primarily investigates the calculation methods for the average optical efficiency and output power of heliostat. Traditionally, calculating shadow and blocking efficiency, key factors affecting the optical performance of mirror fields, has posed a significant technical challenge. Through innovative improvements to the Monte Carlo simulation technique, this study has successfully reduced the spatial dimensionality of the problem. We employed a series of novel coordinate transformations to simulate light tracing on a two-dimensional plane, thereby efficiently calculating the efficiency of shadows and obstructions. Not only have we solved the computational challenges of shadow and obstruction efficiency, but our unique algorithm design also allows us to precisely obtain the overall average optical efficiency of the heliostat field while enhancing computational efficiency.

After 2010, the source model of the microSelectron HDR Afterloader System was slightly modified from the previous model. Granero et al. named the modified source model “mHDR-v2r (revised model mHDR-v2)” and the previous model “mHDR-v2”. They concluded that the dosimetric differences arising from the dimensional changes between the mHDR-v2 and mHDR-v2r designs were negligible at almost all locations (within 0.5 % for r ≥ 0.25 cm), the two-dimensional anisotropy function difference between the two sources is found 2.1 % at r = 1.0 cm when compared with the results of the other experimental group. To confirm this difference, we performed a full Monte Carlo simulation without the energy fluence approximation. This is useful near the radiation source where charged-particle equilibrium does not hold. The two-dimensional anisotropy function of the TG-43U1 dataset showed a few percent difference between the mHDR-v2r and mHDR-v2 sources. There was no agreement in the immediate vicinity from the source (0.10 cm and 0.25 cm), when compared to Granero et al. in mHDR-v2r sources. The differences in these two-dimensional anisotropy functions were identified.

Analyzing network robustness against random failures or malicious attacks is a critical research issue in network science as it helps to enhance the robustness of beneficial networks or efficiently disintegrate harmful networks. Most studies commonly neglect the impact of the attack success rate (ASR) and assume that attacks on the network will always be successful. However, in real-world scenarios, an attack may not always succeed. This paper proposes a novel robustness measure called Robustness-ASR (RASR), which utilizes mathematical expectations to assess network robustness when considering the ASR of each node. To efficiently compute the RASR for large-scale networks, a parallel algorithm named PRQMC is presented, which leverages randomized quasi-Monte Carlo integration to approximate the RASR with a faster convergence rate. Additionally, a new attack strategy named HBnnsAGP is introduced to better assess the lower bound of network RASR. Finally, the experimental results on 6 representative real-world complex networks demonstrate the effectiveness of the proposed methods compared with the state-of-the-art baselines.

Flexible implementations of Monte Carlo Tree Search (MCTS), combined with domain specific knowledge and hybridization with other search algorithms, can be a very powerful for the solution of problems in complex planning. We introduce mctreesearch4j, a standard MCTS implementation written as a standard JVM library following key design principles of object oriented programming. We define key class abstractions allowing the MCTS library to flexibly adapt to any well defined Markov Decision Process or turn-based adversarial game. Furthermore, our library is designed to be modular and extensible, utilizing class inheritance and generic typing to standardize custom algorithm definitions. We demon- strate that the design of the MCTS implementation provides ease of adaptation for unique heuristics and customization across varying Markov Decision Process (MDP) domains. In addition, the implementation is reasonably performant and accurate for standard MDP’s. In addition, via the implementation of mctreesearch4j, the nuances of different types of MCTS algorithms are discussed.

An iterative outlier elimination procedure based on hypothesis testing, commonly known as Iterative Data Snooping (IDS) among geodesists, is often used for the quality control of the modern measurement systems in geodesy and surveying. The test statistic associated with IDS is the extreme normalised least-squares residual. It is well-known in the literature that critical values (quantile values) of such a test statistic cannot be derived from well-known test distributions, but must be computed numerically by means of Monte Carlo. This paper provides the first results about Monte Carlo-based critical value inserted to different scenarios of correlation between the outlier statistics. From the Monte Carlo evaluation, we compute the probabilities of correct identification, missed detection, wrong exclusion, overidentifications and statistical overlap associated with IDS in the presence of a single outlier. Based on such probability levels we obtain the Minimal Detectable Bias (MDB) and Minimal Identifiable Bias (MIB) for the case where IDS is in play. MDB and MIB are sensitivity indicators for outlier detection and identification, respectively. The results show that there are circumstances that the larger the Type I decision error (smaller critical value), the higher the rates of outlier detection, but the lower the rates of outlier identification. For that case, the larger the Type I Error, the larger the ratio between MIB and MDB. We also highlight that an outlier becomes identifiable when the contribution of the measures to the wrong exclusion rate decline simultaneously. In that case, we verify that the effect of the correlation between the outlier statistics on the wrong exclusion rates becomes insignificant from a certain outlier magnitude, which increases the probability of identification.

We numerically study surface models defined on hexagonal disks with a free boundary. 2D surface models for planer surfaces have recently attracted interest due to the engineering applications of functional materials such as graphene and its composite with polymers. These 2D composite meta-materials are strongly influenced by external stimuli such as thermal fluctuations if they are sufficiently thin. For this reason, it is very interesting to study the shape stability/instability of thin 2D materials against thermal fluctuations. In this paper, we study three types of surface models including Landau-Ginzburg (LG) and Helfirch-Polyakov models defined on triangulated hexagonal disks using the parallel tempering Monte Carlo simulation technique. We find that the planer surfaces undergo a first-order transition between the smooth and crumpled phases in the LG model and continuous transitions in the other two models. The first-order transition is relatively weaker compared to the transition on spherical surfaces already reported. The continuous nature of the transition is consistent with the reported results, although the transitions are stronger than that of the reported ones.

Monte Carlo Filtering (MCF) is one of the methods of Experimental Statistical Energy Analysis (E-SEA), which allows the correction of a negative LF (Loss Factor). In this article, a modification of the MCF method, called DESA (Diagonal Extension of the Search Area), is proposed. The technique applies a non-uniform extension of the search area when generating a population of normalized energy matrices. The degree of expansion of the search area is controlled by the Diagonal Penalty Factor (DPF). The authors demonstrated the method's effectiveness on a system that could not be identified in several frequency bands by the classical MCF method. After applying DESA, it was possible to fill in the problematic bands that were missing CLF and DLF values. The paper also proposes a way to minimize the errors introduced by using overly high DPF values.

The rational fabrication of low-dimensional materials with a well-defined topology and functions is an incredibly important aspect of nanotechnology. In particular, the on-surface synthesis (OSS) methods based on the bottom-up approach enable facile construction of sophisticated molecular architectures unattainable by traditional methods of wet chemistry. Among such supramolecular constructs especially interesting are surface-supported metal-organic networks (SMONs), composed of low-coordinated metal atoms and π-aromatic bridging linkers differing in size, shape and functionalization scheme. In this work, the lattice Monte Carlo (MC) simulation technique was used to extract the chemical information encoded in a family of Y-shaped positional isomers coadsorbed with trivalent metal atoms on a flat metallic surface with (111) geometry. For this purpose the investigated tritopic molecules/metal atoms were carefully mapped on a triangular lattice as interconnected/isolated segments with appropriatelly embedded anisotropic interactions. Depending on the intramolecular distribution of active centers (functional groups) within the simulated molecular bricks, we observed the metal-directed self-assembly of two-dimensional (2D) openwork patterns, aperiodic mosaics and metal-organic ladders. The obtained theoretical findings could be especially relevant for the scanning tunneling microscopy (STM) experimentalists interested in a surface-assisted construction of complex nanomaterials stabilized by directional coordination bonds.

In this paper the authors present a method which facilitates computationally efficient parameter estimation of dynamical systems from a continuously growing set of measurement data. It is shown that the proposed method, which utilises Sequential Monte Carlo samplers, is guaranteed to be fully parallelisable (in contrast to Markov chain Monte Carlo methods) and can be applied to a wide variety of scenarios within structural dynamics. Its ability to allow convergence of one's parameter estimates, as more data is analysed, sets it apart from other sequential methods (such as the particle filter).

Dynamic hedging has been adopted by many insurance companies to mitigate the financial risks associated with variable annuity guarantees. In order to simulate the performance of dynamic hedging for variable annuity products, insurance companies rely on nested stochastic projections, which is highly computationally intensive and often prohibitive for large variable annuity portfolios. Metamodeling techniques have recently been proposed to address the computational issues. However, it is difficult for researchers to obtain real datasets from insurance companies to test metamodeling techniques and publish the results in academic journals. In this paper, we create synthetic datasets that can be used for the purpose of addressing the computational issues associated with the nested stochastic valuation of large variable annuity portfolios. The runtime used to create these synthetic datasets would be about 3 years if a single CPU were used. These datasets are readily available to researchers and practitioners so that they can focus on testing metamodeling techniques.

The aim of this study is to present a method for developing fragility curves for soil liquefaction that are consistent with the seismic hazard using Monte Carlo simulation. This method can account for all the uncertainties and variabilities of the input parameters. The seismic parameters of earthquake magnitude (M) and the associated peak ground acceleration (PGA) are considered together for the liquefaction assessment. The liquefaction potential and damages obtained by this method are more realistic. A case study is performed using data from a sand-boil site in Yuanlin, Changhua county, where liquefaction occurred during the Chi-Chi earthquake in Taiwan in 1999. The results show that the liquefaction potential index, IL, the post-liquefaction settlement, St, and the liquefaction probability index, PW, are all suitable parameters for evaluating the liquefaction damages. The fragility curves for soil liquefaction developed by this method can support the performance-based earthquake engineering (PBEE) approach, provide guidance for liquefaction evaluation to the ‘Taiwan Earthquake Loss Estimation System’-TELES, and serve as a basis for scenario simulation and an earthquake early warning system for liquefaction damages.

Accurate treatment planning in radiotherapy essentially decreases secondary damage to healthy tissue surrounding the tumor. Due to plans to use a direct, highly collimated, narrow beam with high intensity to treat small area tumors, researchers have studied microbeam radiation therapy extensively. Using a synchrotron beam as the radiation source may help to limit secondary damage, but treatment planning using computerized simulations and dosimetry is still necessary to achieve optimal results. For this purpose, PDA-gel dosimeters were developed and their sensitivity around a 150 keV induced synchrotron X-ray radiation beam was examined via Monte Carlo simulations using the EGS5 code system. The microbeam development is now at the animal study stage. In this study, we simulate the irradiation of a rat's brain. The simulation results obtained spectra for two types of PDA-gel dosimeters that were compared with the spectrum obtained in a modelized brain tumor of a rat. Additionally, percentage depth dose curves were calculated for the brain tissue and the two gels. Correction equations for the dosimeters were obtained from the dose-difference plots. For further references, these equations can be used to calculate the actual dose in a brain tumor in a rat. The Monte Carlo simulations demonstrate that PDA-gel dosimeters can be used as a reference technique in planning treatment using synchrotron irradiations and in other medical treatment planning methods.

We report the results of the characterization of local Monte Carlo (MC) dynamics of an equilibrium bond fluctuation model polymer matrix (BFM), in time interval typical for MC simulations of non-linear optical phenomena in host-guest systems. The study contributes to the physical picture of the dynamical aspects of quasi-binary mosaic states characterized previously in the static regime. The polymer dynamics was studied at three temperatures (below, above and close to the glass transition), using time-dependent generalization of the static parameters which characterize local free volume and local mobility of the matrix. Those parameters play the central role in the kinetic MC model of host-guest systems. The analysis was done in terms of the probability distributions of instantaneous and time-averaged local parameters. The main result is the characterization of time scales characteristic of various local structural processes. Slowing down effects close to the glass transition are clearly marked. The approach yields an elegant geometric criterion for the glass transition temperature. A simplified quantitative physical picture of the dynamics of guest molecules dispersed in BFM matrix at low temperatures offers a starting point for stochastic modeling of host-guest systems.

Nowadays, the uses of laser and optics in the medical areas are extremely vivid, especially low-level laser therapy. The light with the wavelength of 633 nm to 1200 nm could penetrate and propagate deep in biological tissue. To develop the low-level laser therapy device, optimizing light delivery is critical to accurately stimulate the biological effects inside the biological tissue. Nevertheless, each form of the tissues at each zone on the body had various refractive optic, absorption, scattering, and anisotropy coefficients. This paper describes the simulation results of low-level laser propagation from skin surface at the lower spine, the knee, the femur and the prostate gland with four wavelengths (633 nm, 780 nm, 850 nm, and 940 nm) by the Monte Carlo method. These simulation results are the base for developing the low-level laser therapy device, that could be used in clinical for treating the fracture, knee osteoarthritis, spinal degeneration, and benign prostatic hypertrophy.

The total enthalpies of the 16 different spin configurations that can be realized in the unit cell of the archetype spin crossover complex [Fe(phen)₂(NCS)₂] (phen=1,2-phenanthroline) were calculated applying periodic density functional theory combined with the Hubbard model and the Grimme-D2 dispersion correction (DFT+U/D2). The obtained enthalpy differences between the individual spin configurations were used to determine spin couplings of an Ising-like model and subsequent Monte Carlo simulations for this model allowed to estimate the phenomenological interaction parameter Г of the Slichter-Drickamer model, which is commonly used to describe the cooperativity of the spin transition. The calculational procedure described here, which led to an estimate of about 3 kJ/mol⁻¹ for Г, in good agreement with experiment, may be used to predict from first principles how modifications of spin crossover complexes can change the character of the spin transition from gradual to abrupt and vice versa.

Estimates of extreme water level return periods in river systems are crucial for hydraulic engineering design and planning. Recorded historical water level data of Taiwan’s rivers are not long enough for traditional frequency analyses when predicting extreme water levels for different return periods. In this study, the integration of a one-dimensional flash flood routing hydrodynamic model with the Monte Carlo simulation was developed to predict extreme water levels in the Danshuei River system of northern Taiwan. The numerical model was calibrated and verified with observed water levels using four typhoon events. The results indicated a reasonable agreement between the model simulation and observation data. Seven parameters, including the astronomical tide and surge height at the mouth of the Danshuei River and the river discharge at five gauge stations, were adopted to calculate the joint probability and generate stochastic scenarios via the Monte Carlo simulation. The validated hydrodynamic model driven by the stochastic scenarios was then used to simulate extreme water levels for further frequency analysis. The design water level was estimated using different probability distributions in the frequency analysis at five stations. The design high-water levels for a 200-year return period at Guandu Bridge, Taipei Bridge, Hsin-Hai Bridge, Da-Zhi Bridge, and Chung-Cheng Bridge were 2.90 m, 5.13 m, 6.38 m, 6.05 m, and 9.94 m, respectively. The estimated design water levels plus the freeboard are proposed and recommended for further engineering design and planning.

Up to now, many studies have been performed on particle radiations before or during earthquakes (EQs). In our previous study, with the help of piezoelectricity relationships and the elastic energy formula, the MCNPX simulation code was applied to find the amount of created atomic/ nuclear particles, the dominant interactions; and the energy of the particles for various sizes of quartz and granite blocks. In this study, using the MCNPX simulation code, we have estimated the flux of the particles (created from under-stressed granitic rocks) at different distances from the EQ hypocenter inside the fractures, filled with air, water, and CO2. It was found that inside a water-filled fracture, the particles do not show the flux far from the EQ hypocenter, but inside the gases like air and CO2 with the normal condition density, different types of particles can have a flux far from the source (more than a kilometer) and they might reach themselves to the surface in the case that the EQ hypocenter is very shallow (0- 5 km). However, for deep EQs, it seems that the most detected nuclear particles on the surface should pass via the vacuum-filled fractures and reach themselves to the surface. Moreover, it was concluded that the more density of the fracture’s filling fluid, the less distance that the particles can have a flux.

With the growth of air traffic demand in busy airspace, there is an urgent need for airspace sectorization to increase air traffic throughput and ease the pressure on controllers. The purpose of this paper is to develop a method framework that can perform airspace sectorization automatically, reasonably, which can be used as an advisory tool for controllers as an automatic system, especially for eliminating irregular sector shapes generated by simulated annealing algorithm (SAA) based on region growth method. Two graph cutting method, dynamic Monte Carlo method by changing location of flexible vertices (MC-CLFV) and Monte Carlo method by radius changing (MC-RC) were developed to eliminating irregular sector shapes generated by SAA in post-processing. The experimental results show that the proposed method framework of AS can automatically and reasonably generate sector design schemes that meet the design criteria. Our methodology framework and software can provide assistant design and analysis tools for airspace planners to design airspace, improve the reliability and efficiency of airspace design, and reduce the burden of airspace planners. In addition, this lays the foundation for reconstructing airspace with more intelligent method.

Every firm is devoted to polishing and fine-tuning its maintenance strategy. More exactly, when we develop robust and efficient preventive and conditional maintenance approaches, our objective is to combine performance and robustness to increase maintenance decision-making. This, in turn, leads to a reduction in system failures and a subsequent drop in related expenses. This impulse drives the origin of this essay: to provide fresh criteria that allow us to improve and optimize maintenance techniques. Furthermore, our study has been realistically utilized via the deployment of condition-based maintenance (CBM) techniques, which are commonly regarded as a critical lever that allows firms to gain a competitive advantage in the context of the fourth industrial revolution. Combining both perfect and imperfect maintenance actions into CBM strategies appears to be an effective technique to obtain high economic performance and robustness in the maintenance strategies of the companies. Following this idea, we propose and optimize in this paper a set of criteria enabling the combined assessment of the mean economic performance and the resilience of various kinds of maintenance techniques. The advantage of the proposed criterion is that it adapts to different types of maintenance strategies and provides access to a simple and relevant evaluation model. Thanks to this new criterion, we will be able to choose the better maintenance strategies that have more performance and robustness between these strategies for a system with more or less stable behavior. The Monte Carlo Method is incorporated into the comparative assessment of maintenance methods, adding a layer to the evaluation of the performance and robustness of each adaptation in decision-making.

This work analyzes the severity and risk associated with Automatic Dependent Surveillance-Broadcast (ADS-B) message drop out in Detect and Avoid (DAA) function of Unmanned Aircraft System (UAS). Performance assessment of the Universal Access Transceiver (UAT) ADS-B message implies that in some cases ADS-B fails to update within a specified update interval, which is referred to as drop out in this work. ADS-B is a fundamental surveillance sensor for both class 1 and class 2 DAA systems. Message loss or drop out has been found as one of the common limitations of the ADS-B system. The key feature of this study is incorporating the update rate of real ADS-B data transmitted from the manned aircraft. The data were received from the Grand Forks International Airport, North Dakota. Monte Carlo method has been adopted to resolve encounter scenarios in the presence of drop out. The change in the alert triggered by the UAS DAA in the presence of ADS-B drop out has been investigated. Furthermore, the risk matrices are created to quantify the associated risk with drop out affected alerts. Simulation results depict that both the duration of drop out and DAA look-ahead time affect the alert-triggering function of UAS. With a small look- ahead window and longer duration of drop out, the number of warning alerts increases. Also, alerts are affected more during an overtaking encounter than that of a head-to-head encounter. A system-level analysis is also carried out to recognize the potential reasons behind the ADS-B drop out.

Breaching and liquefaction are two types of coastal flow slides proposed in the literature. While their final results are similar, liquefaction and breaching failure are two different processes, with breaching failure characterized by a different progression and sand movement mechanisms. In densely packed sand, breaching causes delayed sand grain discharge and negative excess pore pressures. Liquefaction, when a mass of soil suddenly behaves like a liquid, flow out over overly gentle slopes. Flow slides typically go unnoticed until the slope fails above ground because they start below the water. Without modern technology, diving equipment, risk assessment, and probabilistic and sensitivity assessments, flow slides inherently challenging. The authors developed a new sensitivity index to detect the susceptibility of coastal flow slides. It was established a sophisticated hybrid model that accommodates flow slides in sync with unpredictable variables employed in this new sensitivity index. Innovative hybrid model has three distinctive models. The Hybrid Hydrodynamic Model addresses wave, wind, current, climate change, and sediment transfer. The Monte Carlo Simulation analyses sensitivity, and the Bayesian Network calculates joint probability of coastal flow slide parameters of this new index that includes environmental characteristics, including climate change. With the assistance of these three models, researchers aim at a) presenting a methodology for coupling coastal flow slide projections with outcomes; b) distinguishing two flow slides; c) examining variables influencing flow slides. The use of such a hybrid model and risk index offers a robust and computationally efficient approach to evaluating the critical angle slope for breaching failure.

Long-range interactions are relevant for a large variety of quantum systems in quantum optics and condensed matter physics. In particular, the control of quantum-optical platforms promises to gain deep insights in quantum-critical properties induced by the long-range nature of interactions. From a theoretical perspective, long-range interactions are notoriously complicated to treat. Here, we give an overview of recent advancements to investigate quantum magnets with long-range interactions focusing on two techniques based on Monte Carlo integration. First, the method of perturbative continuous unitary transformations where classical Monte Carlo integration is applied within the embedding scheme of white graphs. This linked-cluster expansion allows to extract high-order series expansions of energies and observables in the thermodynamic limit. Second, stochastic series expansion quantum Monte Carlo which enables calculations on large finite systems. Finite-size scaling can then be used to determine physical properties of the infinite system. In recent years, both techniques have been applied successfully to one- and two-dimensional quantum magnets involving long-range Ising, XY, and Heisenberg interactions on various bipartite and non-bipartite lattices. Here, we summarise the obtained quantum-critical properties including critical exponents for all these systems in a coherent way. Further, we review how long-range interactions are used to study quantum phase transitions above the upper critical dimension and the scaling techniques to extract these quantum critical properties from the numerical calculations.

Computationally modelling a nuclear reactor startup core for a benchmark against the existing models is highly desirable for an independent assessment informing nuclear engineers and energy policymakers. This work presents a startup core model of the UK’s first Evolutionary Pressurised Water Reactor (EPR) based on Monte Carlo simulations of particle collisions using Serpent 2, a continuous-energy Monte Carlo reactor physics burnup code. Coupling between neutronics and thermal-hydraulic conditions with the fuel depletion is incorporated into the multi-dimensional branches, obtaining the thermal flux and fission rate (power) distributions radially and axially from the three dimensional (3D) single assembly level to a 3D full core. Shannon entropy is employed to characterise the convergence of the fission source distribution, with 3 billion neutron histories tracked by parallel computing. Source biasing is applied for the variance reduction. Benchmarking the proposed Monte Carlo 3D full-core model against the traditional deterministic transport computation suite used by the UK Office for Nuclear Regulation (ONR), a reasonably good agreement within statistics is demonstrated for the safety-related reactivity coefficients, which creates trust in the EPR safety report.

X-ray fluorescence is largely employed in the measurement of the thickness of coatings. Despite of its diffusion, the task is not straightforward because of the complex physics involved that results in high dependence on matrix effects. Thickness quantification is in practice accomplished using the Fundamental Parameters approach, adjusted with empirical measurements of standards with known composition and thickness. This approach has two major drawbacks: i) there are no standards for any possible coating and coating architecture and ii) even relying on standards, the quantification of unknown samples requires the precise knowledge of the matrix nature (e.g., in case of multilayer coatings the thickness and the composition of each underlayer). In this work, we describe a semiquantitative approach to coatings thickness measurement based on the construction of calibration curves through simulated XRF spectra built with Monte Carlo simulations. Simulations have been performed with the freeware software XMI-MSIM. We have assessed the accuracy of the methods by comparing the results with those obtained by i) XRF thickness determination with standards and ii) FIB-SEM cross-sectioning. Then we evaluated which parameters are critical in this kind of indirect thickness measurements.

The presented results of simulations take into account the optical parameters of the selected sea region (from literature data on the southern Baltic Sea) and two optically extreme types of crude oil (from historical data) which exist in the form of a highly diluted oil-in-water emulsion (10 ppm). The spectral index was analyzed based on the results of modelling the radiance reflectance distribution for almost an entire hemisphere of the sky (zenith angle from 0 to 80o). The spectral index was selected and is universal for all optically different types of oil (wavelengths 650 and 412 nm). The possibility of detecting pollution in the conditions of the wavy sea surface (as a result of wind of up to 10 m/s) was studied. It has been also shown that if the viewing direction is close to a direction perpendicular to the sea surface, observations aimed at determining the spectral index are less effective than observation under the zenith angle of incidence of sunlight for all azimuths excluding the direction of sunlight specular reflection.

New setting is introduced to study quasi-degree and quasi-co-degree arising from co-neighborhood. quasi-degree and quasi-co-degree is about a vertex which are applied into the setting of neutrosophic graphs. . The structure of set is studied and general results are obtained. Also, some classes of neutrosophic graphs namely path-neutrosophic graphs, cycle-neutrosophic graphs, complete-neutrosophic graphs and star-neutrosophic graphs, complete-bipartite-neutrosophic graphs and complete-multipartite-neutrosophic graphs are investigated in the terms of a vertex which is called either quasi-degree or quasi-co-degree. Neutrosophic number is reused in this way. It’s applied to use the type of neutrosophic number in the way that, three values of a vertex are used and they’ve same share to construct this number to compare with other vertices. Summation of three values of vertex makes one number and applying it to a comparison. This approach facilitates identifying vertices which form quasi-degree and quasi-co-degree. Quasi-degree is a value of a vertex which is maximum amid all values of vertices which are neighbors to a fixed vertex. Quasi-co-degree is a value of an edge which is maximum amid all values of edges which are neighbors to a fixed vertex but corresponded vertex is representative for this notion. Using different values which are related to a vertex inspire us to focus on edge and vertices which are corresponded to a fixed vertex. The notion of neighborhood is used to collect either vertices are titled neighbors or edges are incident to fixed vertex. In both settings, some classes of well-known neutrosophic graphs are studied. Some clarifications for each result and each definitions are provided. Using fixed vertex has key role to have these notions in the form of vertex or edge. The value of an edge has eligibility to call quasi-co-degree but the value of a vertex has eligibility to call quasi-degree. Some results get more frameworks and perspective about these definitions. The way in that, two vertices have connection together, open the way to define neighborhood and co-neighborhood. The maximum values in neighborhood and co-neighborhood introduces quasi-degree and quasi-co-degree, respectively. New name is chosen from degree. Since amid all vertices with different degrees, one vertex is chosen. In other words, one vertex is fixed and its degree turns out quasi-degree where two degrees could be assigned to a vertex. Degree of edges and degree of vertices. The number of edges which are incident to the vertex and the number of vertices which are neighbors to the vertex. Degree and co-degree are the notions which are transformed to use in quasi-style. Two neutrosophic values introduce two neutrosophic vertices separately in each settings. These notions are applied into neutrosophic graphs as individuals but not family of them as drawbacks for these notions. Finding special neutrosophic graphs which are well-known, is an open way to purse this study. Some problems are proposed to pursue this study. Basic familiarities with graph theory and neutrosophic graph theory are proposed for this article.

New setting is introduced to study co-neighborhood, neutrosophic t-neighborhood, neutrosophic quasi-vertex set, neutrosophic quasi-order, neutrosophic neighborhood, neutrosophic co-t-neighborhood, neutrosophic quasi-edge set, neutrosophic quasi-size, Neutrosophic number, neutrosophic co-neighborhood, co-neutrosophic number, quasi-number and quasi-co-number. Some classes of neutrosophic graphs are investigated.

The technology for determining the point’s coordinates on the earth’s surface using the global navigation satellite system (GNSS) is becoming the norm along with ground-based methods. In this case, determining coordinates does not cause any particular difficulties. However, to identify normal heights using this technology with a given accuracy, special research is required. The fact is that satellite determinations of geodetic heights (h) over an ellipsoid surface differ from ground-based measurements of normal height (HN) over a quasi-geoid surface by a certain value called quasi-geoid height or height anomaly (ζ). In relation to determining heights of a certain territory, the concept of geoid height (N) is usually operated when dealing with a geoid model. In this work, geodetic and normal heights are determined for 5 control points in three different regions in Lebanon, where measurements are carried out using GNSS technology and geometric levelling. The obtained quasi-geoid heights are compared with geoid heights derived from the global Earth model EGM2008. The results obtained showed that, in the absence of gravimetric data, the combination of global Earth model data, geometric levelling for selected areas and satellite determinations allows the creation of a high accurate altitude network for mountainous areas.

Analytical solutions for radiation dominated phase of Quasi Steady State Cosmology in Friedmann-Robertson-Walkar models are obtained. We find that matter density is positive in all the cases (k = 0,-1,1). The nature of Hubble parameter (H) in [0,2] is discussed. The deceleration parameter (q) is marginally less than zero indicating accelerating universe. The scale factor (S) is graphically shown with time. The model represents oscillating universe between the above mentioned limits. Because of bounce in QSSC, the maximum density phase is still matter dominated. The models represent singularity free model. We also find that the models have event horizon i.e. no observer beyond the proper distance rH can communicate each other in FRW mdels for radiation dominated phase in the frame work of QSSC. The FRW models are special classes of Bianchi type I, V, IX space-times with zero, negative and positive curvatures respectively. Initially i.e. at  = 0, the model represents steady model. We have tried to show how a good fit can be obtained to the observations in the framework of QSSC during radiation dominated phase.

We prove a quasi-pseudometric version of the strong Ekeland principle proved by P. Georgiev and T. Suzuki.

This paper introduces an innovative Particle Swarm Optimization (PSO) Algorithm in-corporating Sobol and Halton Random number samplings. It evaluates the enhanced PSO's performance against conventional PSO employing Monte Carlo Random Number Samplings. The comparison involves assessing the algorithms across nine benchmark problems and the renowned Travelling Salesman Problem (TSP). The re-sults reveal consistent enhancements achieved by the enhanced PSO utilizing Sob-ol/Halton samplings across the benchmark problems. Particularly noteworthy are the substantial improvements demonstrated by the enhanced PSO in solving the Travel-ling Salesman Problem. These findings underscore the efficacy of employing Sobol and Halton random number generation methods in enhancing algorithm efficiency.

With the advent of submesoscale O(1km) permitting basin-scale ocean simulations, the seasonality in the mesoscale O(50km) eddies with kinetic energies peaking in summer has been commonly attributed to submesoscale eddies feeding back onto the mesoscale via an inverse energy cascade under the constraint of stratification and Earth’s rotation. In contrast, by running a 101-member, seasonally forced, three-layer quasi-geostrophic (QG) ensemble configured to represent an idealized double-gyre system of the subtropical and subpolar basin, we find that the mesoscale kinetic energy shows a seasonality consistent with the summer peak without resolving the submesoscales; by definition, a QG model only resolves small Rossby number dynamics (O(Ro)≪1) while as submesoscale dynamics are associated with O(Ro)∼1. Here, by quantifying the Lorenz cycle of the mean and eddy energy, defined as the ensemble mean and fluctuations about the mean respectively, we propose a different mechanism from the inverse energy cascade by which the stabilization and strengthening of the western-boundary current during summer due to increased stratification leads to a shedding of stronger mesoscale eddies from the separated jet. Conversely, the opposite occurs during the winter; the separated jet destablizes and results in overall lower mean and eddy kinetic energies despite the domain being more susceptible to baroclinic instability from weaker stratification.

Systematic traders employ algorithmic strategies to manage their investments. As a result of the deterministic nature of such strategies, it is possible to determine their exact responses to any conceivable set of market conditions. Consequently, sensitivity analysis can be conducted to systematically uncover undesirable strategy behavior and enhance strategy robustness by adding controls to reduce exposure during periods of poor performance / unfavorable market conditions or increase exposure during periods of strong performance / favorable market conditions. In this study, we formulate both a simple systematic trend-following strategy (i.e., trading model) to simulate investment decisions, and a market model to simulate the evolution of instrument prices. We then map the relationship between market model parameters under various conditions and strategy performance. We focus, in particular, on identifying the performance impact of changes in both serial dependence in price variability and changes in the trend. The long-range serial dependence of the true range worsens performance of the simple classic trend-following strategy. During periods of strong performance, the dispersion of trading outcomes increases significantly as long-range serial dependence increases.

A quantum particle constrained between two high potential barriers provides a paradigmatic example of a system sustaining quasi bound (or resonance) states. When the system is prepared in one of such quasi bound states, the wave function approximately maintains its shape but decays in time in a nearly exponential manner radiating into the surrounding space, the lifetime being of the order of the reciprocal of the width of the resonance peak in the transmission spectrum. Naively, one could think that adding more lateral barriers would preferentially slow down or prevent the quantum decay since tunneling is expected to become less probable and because of quantum backflow induced by multiple scattering processes. However, this is not always the case and in the early stage of the dynamics quantum decay can be accelerated (rather than decelerated) by additional lateral barriers, even when the barrier heights are arbitrarily large. The decay acceleration originates from resonant tunneling effects and is associated to large deviations from an exponential decay law. We discuss such a counterintuitive phenomenon by considering the hopping dynamics of a quantum particle on a tight-binding lattice with on-site potential barriers.

In this study, we present more than 8 years of observations of the quasi-16-day wave (Q16DW) in the mesosphere and lower thermosphere (MLT) wind at middle latitudes observed by the Mengcheng (33.4°N, 116.5°E) meteor radar. The long-term variation in amplitudes calculated from the data between April 2014 and December 2022 shows enhanced wave activity during winter and early spring (near equinox) and suppressed wave activity during the summer. The Q16DWs are relatively weak in the meridional wind. During the winter months, the Q16DWs in the zonal component exhibit a burst below 85 km, and their amplitudes reach up to 10 m/s. In the early spring, the Q16DWs strengthen above 90 km with amplitudes in excess of 12 m/s. The phase differences between the zonal and meridional components of the Q16DW are, on average, slightly smaller than 90°, suggesting the existence of orthogonal relationships between them. During strong bursts, the periods of the Q16DW in winter range between 15 and 18 d, whereas in winter, the periods tend to be more diffuse. The wintertime Q16DW is amplified, on average, when the zonal wind shear peaks, suggesting that barotropic instability may be one source of Q16DW. Q16DW amplitudes exhibit considerable interannual variability; however, a relationship between the 11-year solar cycle and the Q16DW is not found.

The purpose of this work is to study Rossby wave parameters in total ozone over Arctic in 2000–2021. We consider the averages in the January–March period, when stratospheric trace gases (including ozone) in sudden stratospheric warming events are strongly disturbed by planetary waves. To characterize the wave parameters, we analyzed ozone data at the latitudes of 50° N (the sub-vortex area), 60° N (the polar vortex edge) and 70° N (inner region of the polar vortex). Total ozone column (TOC) measurements during 22-year time interval were used from Total Ozone Mapping Spectrometer (TOMS) / Earth Probe and Ozone Mapping Instrument (OMI) / Aura satellite observations. The total ozone zonal distribution and variations in the parameters of the Fourier spectral components with zonal wave numbers m = 1–5 are presented. Daily and interannual variations in TOC, amplitudes and phases of spectral wave components, and linear trends of the quasi-stationary wave 1 (QSW1) amplitudes are discussed. The positive TOC peaks inside the vortex in 2010 and 2018 alternate with negative ones in 2011 and 2020. The latter TOC anomalies correspond to severe depletion of stratospheric ozone over the Arctic in the strong vortex conditions due to anomalously low activity of planetary waves. Variations in TOC in sub-vortex region exhibit the statistically significant negative trend –4.8±5.4 DU decade–1 in QSW1 amplitude, while the trend is statistically insignificant at the vortex edge region due to increased TOC variability. Processes associated with polar vortex dynamics are discussed, including quasi-stationary vortex asymmetry and quasi-circumpolar migration of the wave-1 phase at the vortex edge.

The Boltzmann-Gibbs (BG) entropy has been used in a wide variety of problems for more than a century. It is well known that BG entropy is extensive, but for certain systems such as those dictated by long-range interactions, the entropy must be non-extensive. Tsallis entropy possesses non-extensive characteristics, which is parametrized by a variable q (q = 1 being the classic BG limit), but unless q is determined from microscopic dynamics, the model remains but a phenomenological tool. To this date very few examples have emerged in which q can be computed from first principles. This paper shows that the space plasma environment, which is governed by long-range collective electromagnetic interaction, represents a perfect example for which the q parameter can be computed from micro-physics. By taking the electron velocity distribution function measured in the heliospheric environment into account, and considering them to be in quasi equilibrium state with electrostatic turbulence known as the quasi-thermal noise, it is shown that the value corresponding to q = 9/13 = 0.6923 may be deduced. This prediction is verified against observation made by spacecraft, and it is shown to be in excellent agreement.

Enhancement in chemisorption is one of the active research areas in carbon materials. To remedy the thermally degraded chemisorption at high temperatures, we report a comprehensive study of kink structures in free-standing monoatomic carbon nanowires upon heating. Our Monte Carlo simulation considers multi-monoatomic carbon chains laterally interacted by the Van der Waals force. Our study reveals that the carbon nanowires maintain their linearity regardless of chain length at low temperatures, but this is not the case at high temperatures. Disordered kink structure is observed in the short carbon chains, especially above the Peierls transition temperature. The severe kink structure might increase the possibility of attaching negatively charged atoms, thereby contributing to the development of next-generation materials for chemisorption at high temperatures. We have also provided an important inspection that any physical property of the finite-length carbon chain predicted by ab-initio calculation should reconsider the atomic rearrangement due to the thermal instability at high temperatures.

Mobile robots need an environmental perception ability in order to interact with the surrounding environment. In this paper we present the G-SLAM method, where the map is consisted of a cloud of scattered points in the continuous space and each point is accompanied by an obstacle existence probability. On the other hand the robot's pose (and trajectory) is estimated by an particle filters while the cloud of points is estimated by an adaptive recursive algorithm which is presented. Main feature of the presented method is the adaptive repositioning of the scattered points and their convergence around obstacles. In addition to our previous publications, In this paper a mathematical derivation of the method is presented and real case scenarios are discussed, presented and compared with other methods.

The practice of spatial econometrics revolves around a weighting matrix, which is often supplied by the user on previous knowledge. This is the so called $\mathbf{W}$ issue. Probably, the aprioristic approach is not the best solution although, nowadays, there few alternatives for the user. Our contribution focuses on the problem of selecting a $\mathbf{W}$ matrix from among a finite set of matrices, all of them considerer appropriate for the case. We develop a new and simple method based on the Entropy corresponding to the distribution of probability estimated for the data. Other alternatives, which are common in current applied work, are also reviewed. The paper includes a large Monte Carlo to calibrate the effectiveness of our approach compared to the others. A well-known case study is also included.

In many research fields, statistical probability models are often used to analyze real-world data. However, data from many fields, such as the environment, economics, and health care, do not necessarily fit traditional models. New empirical models need to be developed to improve their fit. In this paper, we explore a further extension of the quasi-Lindley model. Maximum likelihood, least square error, Anderson-Darling algorithm, and expectation-maximization algorithm are four techniques for estimating the parameters under study. All techniques provide accurate and reliable estimates of the parameters. However, the mean square error of the expectation maximization approach was lower. The usefulness of the proposed model was demonstrated by analyzing a dynamical systems data set, and the analysis shows that it outperforms the other models in all statistical models considered.

In this article, we analyze the thermodynamics, thermal fluctuations and quasi-normal modes of a torus-like charge AdS black hole. The parameters of a torus-like charged AdS black hole are its mass, electric charge, and cosmological constant. We used a torus-like charge AdS black hole to evaluate thermodynamics through thermal fluctuations and found that it significantly affected the corrected thermodynamic potentials. Black holes of larger radii are more stable and significant. The effect of thermal fluctuation has been investigated on the other thermodynamic potentials such as the Hemohltz free energy, Gibbs free energy, enthalpy, and specific heat, using the concept of corrected entropy. Finally, we use null geodesics to deal with quasi-normal modes. The null geodesics provide the photon sphere angular velocity and Lyapunov exponent, which are the real and imaginary components of the quasi-normal modes in the eikonal limit. With inceasing charge and photon radius, angular velocity reaches its maximum, similar to Lyapunov's exponent.

Mesoscale eddies, the weather system of the oceans, although being on the scales of O(20-100km), have a disproportionate role in shaping the mean jets such as the separated Gulf Stream in the North Atlantic Ocean, which is on the scale of O(1000km) in the along-jet direction. With the increase in computational power, we are now able to partially resolve the eddies in basin-scale and global ocean simulations, a model resolution often referred to as mesoscale permitting. It is well known, however, that due to grid-scale numerical viscosity, mesoscale permitting simulations have less energetic eddies and consequently weaker eddy feedback onto the mean flow. In this study, we run a quasi-geostrophic model at mesoscale resolving resolution in a double gyre configuration and formulate a deterministic parametrization for the eddy rectification term of potential vorticity (PV), namely, the eddy PV flux divergence. We have moderate success in reproducing the spatial patterns and magnitude of eddy kinetic and potential energy diagnosed from the model. One novel point about our approach is that we account for non-local eddy feedbacks onto the mean flow by solving the eddy PV equation prognostically in addition to the mean flow. In return, we are able to parametrize the variability in total (mean+eddy) PV at each time step instead of solely the mean PV. A closure for the total PV is beneficial as we are able to account for both the mean state and extreme events.

The influence of the grain structure on the tensile deformation strength is studied for precipitation-strengthened Cu-0.2%Zr at 673 K. Subgrains and grains are formed by ECAP and annealing. The fraction of high-angle boundaries increases with prestrain. Subgrains and grains coarsen during deformation. This leads to softening in the quasi-stationary state. The initial quasi-stationary state of severely predeformed, ultrafine-grained material exhibits relatively high rate-sensitivity at relatively high stresses. This is interpreted as result of the stress dependences of the quasi-stationary subgrain size and the volume fraction of subgrain-free grains.

This research work deals with the modeling and control of hybrid photovoltaic (PV) - Wind micro-grid using Quasi Z-Source inverter. This inverter provides better buck/boost characteristics, able to regulate the phase angle output, less harmonic content, no requirement of the filter and has high power performance characteristics over conventional inverter as major benefits. A SEPIC converter as dc-dc switched power apparatus is employed for maximum power point tracking (MPPT) functions which provides high voltage gain throughout the process. Moreover, a modified power ratio variable step (MPRVS) based perturb & observe (P&O) method has been proposed in the PV MPPT action which forces the operating point close to maximum power point (MPP). Practical responses justify the performance of hybrid PV-Wind micro-grid with Quasi Z-Source inverter structure.

Ultrafast pump-probe spectroscopic studies allow for deep insights into the mechanisms and timescales of photophysical and photochemical processes. Extracting valuable information from these studies, such as reactive intermediates' lifetimes and coherent oscillation frequencies, is an example of the inverse problems of chemical kinetics. This article describes a consistent approach for solving this inverse problem that avoids the common obstacles of simple least-squares fitting that can lead to unreliable results. The presented approach is based on the regularized Markov Chain Monte Carlo sampling for the strongly nonlinear parameters, allowing for a straightforward solution of the ill-posed nonlinear inverse problem. The software to implement the described fitting routine is introduced and the numerical examples of its application are given. We will also touch on critical experimental parameters, such as the temporal overlap of pulses and cross-correlation time and their connection to the minimal reachable time resolution.

In this paper, we consider a class of bilevel programming problems (abbreviated as, BLPP). Exploiting the generalized approximate convexity assumptions, we investigate the relations among the solutions of approximate Minty (respectively, Stampacchia) type variational inequalities (abbreviated as, AMTVI (respectively, ASTVI)), and the local ϵ-quasi solutions of the BLPP. Moreover, by employing the generalized Knaster–Kuratowski–Mazurkiewicz (abbreviated as, KKM)-Fan’s lemma, we derive some existence results for the solutions of approximate variational inequalities (abbreviated as, AVI), namely, AMTVI and ASTVI. A non-trivial example is given to highlight the importance of the established results. To the best of our knowledge, there is no research paper available in the literature that establishes relationships between the AVI and the BLPP under the assumptions of generalized approximate convexity in terms of limiting subdifferentials.

Polymer-based composites are widely used in the automotive, security, aeronautical and space industries, to mention a few. This is because of their good mechanical properties, which are similar to those of metals but with the attraction of being lightweight. Kevlar® as a reinforcement is extensively used in the security industry owing to its good ballistic properties. This investigation presents the mechanical characterization based on in-plane quasi-static tensile testing of Kevlar® composite using a universal testing system. The methodology developed for the preparation of coupons is based on pressure, temperature and time. As a result of this work, elastic moduli (EL and ET), Poisson’s ratio (νLT), shear modulus (GLT) and strengths (XT, YT, S) were obtained. It is worth mentioning that there is scarce or no characterization of this material in the literature, and those studies that do characterize it do not present the coupon thermoforming conditions or the reasons of the specimen dimensions (width, length and thickness).

This work presents a new type of chaotic dynamical system originated from Dirichlet Eta function. Although the exposition describes only a bidimensional instance, it is certainly possible to construct higher dimensional such models. The novelty resides in the multitude of interesting characteristics concentrated in only one kind of dynamical system, mainly because of its simple abstract definition - it provokes dissipative, chaotic and itinerant behavior, and its paths typically display nonrecurring quasi-attractors and associate regions of ruin. In addition, it is discrete and nonautonomous, being additively driven by a simple real vector sequence. Also, the defining formulas have parameters which permit to activate bifurcatory effects, mainly related to setting the number of quasi-attractors and overall geometric complexity and diversity. In this fashion, it holds a certain similarity with classical itinerant chaotic systems, but it is a completely different object. One of its parameters is very influential and may be intuitively compared to a type of "potential energy", capable of changing the duration of sojourns etc. In this fashion, there is a potential opportunity to employ artificial inference techniques for designing specific structures on paths. Of course, the other parameters are also relevant, but not necessarily in the control of the number of quasi-attractors and the permanence periods near them. It is important to state that the described system was found during a previous investigation about certain features of the Riemann hypothesis, using the Dirichlet Eta function instead of the Riemann Zeta function - the components of the system under study here were part of certain series whose convergence had to be established, and that scenario produced this very positive side effect. The text also suggests an atypical, but important, potential application in philosophy: the representation of the Hindu concepts of Samsara and Moksha.

The fundamental theory of nuclear interactions, Quantum Chromodynamics (QCD), operates in terms of quarks and gluons at higher resolution. At low resolution the relevant degrees of freedom are nucleons. Two-nucleon Short-Range Correlations (SRC) help to interconnect these two descriptions. SRCs are temporary fluctuations of strongly interacting close pairs of nucleons. The distance between the two nucleons is comparable to their radii and their relative momenta are larger than the fermi sea level. According to the electron scattering experiments held in the last decade, SRCs have far-reaching impacts on many-body systems, the nucleon-nucleon interactions, and nuclear substructure. The modern experiments with ion beams and cryogenic liquid hydrogen target make it possible to study properties of the nuclear fragments after quasi-elastic knockout of a single nucleon or an SRC pair. The first SRC experiment at BM@N (2018) with a carbon beam has shown that detection of an intact 11B nucleus after interaction selects out the quasi-elastic knockout reaction with minimal contribution of initial- and final-state interactions. Also, 25 events of SRC-breakups showed agreement in SRC properties as knows from electron beam experiments. The analysis of the second measurement of SRC at BM@N held in 2022 with an improved setup is currently ongoing. The SRC project at JINR moved to a new experimental area in 2023, where the next measurement is being planned in terms of experimental setup and physics goals.

In this paper, a fuzzy logic-based energy management system and multi-port quasi-z-source converter are proposed that utilizes hybrid renewable energy sources. The system ensures that each energy source module can be used individually by employing fuzzy logic to define power modes. This approach also helps prevent switching losses resulting from extra switching of the source modules. Also, the proposed energy management does not have a mathematical model, so applicability is simple and it is suitable for different multiple input topologies. The Mamdani fuzzy inference system can be designed to capture the nonlinear behavior of the system because of the linguistic rules. Moreover, by adopting the quasi-z-source network to the end of the multi-port converter, the switching losses of the multi-port modules have been significantly reduced. The simulation results are obtained using MATLAB Simulink, and the experimental results are obtained by connecting the circuit to MATLAB Simulink via an Arduino Due.

In this work, the recently proposed cracking elements method (CEM) is used for simulating the damaging processes of structures with initial imperfections. CEM is built in the framework of conventional FEM which is formally like a special type of finite element. The disconnected piecewise cracks are used for representing crack paths. Taking the advantages of CEM that both the initiations and propagations of cracks can be naturally captured, we numerically study the uni-axial compression tests of specimens with multiple joints and fissures, where the cracks may propagate from the tips, or from some other unexpected positions. Though uni-axial compression tests are considered, mainly tensile damage criteria are used in the numerical model. On one hand, the results demonstrate the robustness and effectiveness of the CEM while on the other hand, some drawbacks of the present model are demonstrated, showing the future work.

The aim of this paper is to discuss the relations between various notions of sequential completeness and the corresponding notions of completeness by nets or by filters in the setting of quasi-metric spaces. We propose a new definition of right $K$-Cauchy net in a quasi-metric space for which the corresponding completeness is equivalent to the sequential completeness. In this way we complete some results of R.~A. Stoltenberg, Proc. London Math. Soc. \textbf{17} (1967), 226--240, and V.~Gregori and J.~Ferrer, Proc. Lond. Math. Soc., III Ser., \textbf{49} (1984), 36. A discussion on nets defined over ordered or pre-ordered directed sets is also included.

One of the main challenges in geodetic deformation analysis is to infer whether the geometry of some engineering structures or zones of natural hazards has changed from their initial state or not. The pillar that supports such an analysis is tightly based on the fundamentals of statistical hypothesis testing. The null hypothesis model indicates that no displacement has occurred. It is tested against a class of alternative models, which stipulate different displacement patterns. In this contribution, we present an innovative geodetic displacement detection which integrates combinatorial analysis and likelihood ratio tests into a sequential procedure for the case where the differences between observations of two epochs are in play. This framework is applied and investigated to two test scenarios: a synthetic and a real simulated trilateration network. From a statistical point of view, our approach is rigorous, because the alternative model can identify simultaneously more than one unstable point. In addition, the relationship between the unknown parameters and the observations is always linear, even if the problem manifests itself as non-linear. Consequently, we avoid a potential loss of statistical test power due to the model linearization. In addition, the proposed method controls the false positive rate efficiently. One of the problems that arise here is also related to the selection of the maximum number of points to be considered in the test procedure ( ). Here, we provide an innovative methodology based on rank computation of the design matrices to define , which can even be extended to the problem of outlier detection. Determining avoids the problem of having non-separable models in identifying unstable points. The algorithms and data are available in the repository.

In this study, we present different methods of estimating fuzzy reliability of a two-parameter Rayleigh distribution via the maximum likelihood estimator, median first-order statistics estimator, quartile estimator, L-moment estimator, and mixed Thompson-type estimator. The mean-square error MSE as a measurement for comparing the considered methods using simulation through different values for the parameters and unalike sample sizes is used. The results of simulation show that the fuzziness values are better than the real values for all sample sizes, as well as the fuzzy reliability at the estimation of the Maximum

In order to improve the charge collection efficiency, we simulated and experimentally tested the doping structure of the electron multiplication layer in EBCMOS. In this paper, we simulate the charge collection efficiency of EBCMOS under different doping methods by modeling the collisional scattering of electrons with solid atoms in semiconductor materials and combining the transport trajectories of electrons in the electron multiplication layer, the simulation results indicate that using a layered doping structure to reduce the thickness of the index heavily doped layer can effectively optimize the electric field distribution in the electron multiplier layer and reduce the recombination rate of electrons. The optimized doping structure has a significant effect on improving charge collection efficiency. Based on the simulation results, doped samples were prepared and tested. The test results showed that the charge collection efficiency obtained under the condition of a P-type silicon substrate thickness of 7μm and an index heavily doped layer thickness of 1μm was 72.65%, reducing the thickness of the index heavily doped layer to 0.1μm, the charge collection efficiency obtained can reach 86.27%, which proves that reducing the thickness of the index heavily doped layer can effectively improve the charge collection efficiency of EBCMOS devices.

The quantile estimation of extreme wind speed is needed in various environmental fields such as climatology, design of structures, renewable energy sources and agricultural operations. These calculations are crucial for the coding of wind speed. In this study, the required wind speed series of 16 stations in Khyber Pakhtunkhwa, Pakistan, was obtained from the NASA official website and measured in meters per second (m/s) at a 10-meter distance. A Regional Frequency Analysis of 16 AMWS stations was performed using L-moments. The quantile estimates of extreme wind speed are needed for various areas of interest using Regional Frequency Analysis (RFA) and extreme value theory. These calculations are crucial for the coding of wind speed. The data was taken from the NASA official website at a 10-meter distance and measured in meter per second (m/s). A Regional Frequency Analysis of AMWS using L-moments is performed utilizing wind speed data from sixteen sites (16) in Pakistan's Khyber Pakhtunkhwa province. There are no sites that are found to be discordant. The wards method is used to construct a homogenous region and make two homogenous regions from 16 sites. The heterogeneity test justifies that both clusters are homogeneous. The most appropriate probability distribution from the Generalized Normal (GNO), Generalized Logistic (GLO), Pearson Type-3 (P3), Generalized Pareto (GPA), and Generalized Extreme Value (GEV) distributions are chosen to calculate regional quantiles. According to the L-moments diagram and Z statistics, GEV for Cluster- Ι and GLO for Cluster- ΙΙ are the best suggestions from the others. Both clusters’ robustness is measured utilizing Relative Bias (RB) and Relative Root Mean Square Error (RRMSE). Overall, GEV distribution is fit for cluster-Ι, and the GLO distribution is fit for cluster-ΙΙ. Utilizing the site mean and median as index parameters, we can also find at-site quantiles from regional quantiles. The study’s quantile estimates can be employed in codified structural designs with policy consequences.

This research evaluates the inhibitory effect of L-amino acids (AAs) with different side chain lengths on Fe surface implementing Monte Carlo (MC) simulation. A quantitative and qualitative description of the adsorption behavior of AAs on the iron surface has been carried out. Calculations have shown that the absolute values of the adsorption energy of L-amino acids increase with side chain prolongation; they are also determined by the presence of heteroatoms. AAs from nonpolar and basic groups have the best adsorption ability to the iron surface, which indicates their highest inhibitory efficiency according to the results of MC simulation.

Muography is an expanding technique for the investigation of the internal structure of targets of interest in geophysics. The flux of high penetrating muons produced by primary cosmic rays is attenuated by traversing kilometer size objects like X-ray flux is attenuated through the human body. This gives the possibility to study the internal structure of volcanoes or underground cavities, e.g., starting from the measure of the muon flux transmission through the target. Many groups of researchers working with this technique have to face with high background level that afflicts the reconstruction of muon tracks near the horizontal direction. An important source of background is the scattering of particles near the detector that produces a wrong reconstruction of the incoming direction. An innovative technique based on Cherenkov radiation was investigated by means of Monte Carlo simulations developed in Geant4 toolkit and MATLAB. The results reported in this paper show that the presented technique is able to correctly distinguish the incoming direction of particles with an efficiency higher than 98%. This new kind of detector could represent an alternative for high resolution charged particle tracking also for other applications.

Evaluating new, promising organic molecules to make next-generation organic optoelectronic devices necessitates the evaluation of charge carrier transport performance through the semi-conducting medium. In this work, we utilize quantum chemical calculations (QCC) and kinetic Monte Carlo (KMC) simulations to predict the zero-field hole mobilities of ~100 morphologies of the benchmark polymer poly(3-hexylthiophene), with varying simulation volume, structural order, and chain-length polydispersity. Morphologies with monodisperse chains were generated previously using an optimized molecular dynamics force-field and represent a spectrum of nanostructured order. We discover that a combined consideration of backbone clustering and system-wide disorder arising from side-chain conformations are correlated with hole mobility. Furthermore, we show that strongly interconnected thiophene backbones are required for efficient charge transport. This definitively shows the role "tie-chains" play in enabling mobile charges in P3HT. By marrying QCC and KMC over multiple length- and time-scales, we demonstrate that it is now possible to routinely probe the relationship between molecular nanostructure and device performance.

The present work is concerned with the study of the generalized Langevin equation and its link to the physical theories of statistical mechanics and scale relativity. It is demonstrated that the form of the coefficients of Langevin equation depend critically on the assumption of continuity of the reconstructed trajectory. This in turn demands for the fluctuations of the diffusion term to be discontinuous in time. This paper further investigates the connection between the scale-relativistic and stochastic mechanics approaches, respectively, with the study of the Burgers equation, which in this case appears as a stochastic geodesic equation. By further demanding time reversibility of the drift the Langevin equation can also describe equivalent quantum-mechanical systems in a path-wise manner. The resulting statistical description obeys the Fokker-Plank equation of the probability density of the differential system, which can be readily estimated from Monte Carlo simulations of the random paths. Based on the Fokker-Plank formalism a new derivation of the transient probability densities is presented. Finally, stochastic simulations are compared to the theoretical results.

We develop an interactive likelihood (ILH) for sequential Monte-Carlo (SMC) methods for image-based multiple target tracking applications. The purpose of the ILH is to improve tracking accuracy by reducing the need for data association. In addition, we integrate a recently developed deep neural network for pedestrian detection along with the ILH with a multi-Bernoulli filter. We evaluate the performance of the multi-Bernoulli filter with the ILH and the pedestrian detector in a number of publicly available datasets (2003 PETS INMOVE, AFL, and TUD-Stadtmitte) using standard, well-known multi-target tracking metrics (OSPA and CLEAR MOT). In all datasets, the ILH term increases the tracking accuracy of the multi-Bernoulli filter.

Quasi-interpolation is a powerful tool for approximating functions using radial basis functions (RBFs) such as Gaussian kernels. This avoids solving large systems of equations as in RBF interpolation. However, quasi-interpolation with Gaussian kernels on compact intervals can have significant errors near the boundaries. This paper proposes a quasi-interpolation method with Gaussian kernels using Chebyshev points and boundary corrections to improve the approximation near the boundaries. The boundary corrections use a linear approximation of the function beyond the interval to estimate the truncation error and add correction terms. Numerical studies on test functions show that the proposed method reduces errors significantly near boundaries compared to quasi-interpolation without corrections, for both equally spaced and Chebyshev points. The convergence and accuracy with the boundary corrections are generally better with Chebyshev points compared to equally spaced points. The proposed method provides an efficient way to perform quasi-interpolation on compact intervals while controlling the boundary errors.

The objective of this study is to examine a nonparametric estimate , using the kernel approach, of the conditional distribution function of a scalar response variable that is given a random variable whose values take place in a separable real Hilbert space. The observations will be dependent on one another in a quasi-associated fashion. The pointwise practically perfect consistencies with rates of this estimator are established by us under some broad conditions. The study’s major objective is to investigate the convergence rate of the proposed estimator and its application in the convergence rate and asymptotic normality of the hazard function. The asymptotic normality of the developed estimator is established precisely. Simulation studies were conducted to investigate the behavior of the asymptotic property in the context of finite sample data.

A structural study of the 6√3 reconstruction conversion on 4H-SiC substrate surface in an Ar medium with a short sublimation annealing time into quasi-free epitaxial graphene by means of hydrogen intercalation between the upper layer of SiC and the adjacent reconstructed layer was carried out using the technique of reflection high-energy electron diffraction (RHEED) characterizing the surface of a crystalline material. A slight violation of the 6√3 reconstruction layer formation uniformity was found. The results of the study of the crystal structure of quasi-free graphene and single-layer graphene comprising a buffer layer formed on 4H-SiC in the traditional way in an Ar atmosphere without intercalation were compared.

We present set-theoretical solutions of the Yang-Baxter equation in BCK–algebras. Some solutions in BCK−algebras are not solutions in other structures (such as MV −algebras). Related to our investigations we also consider some new structures: Boolean coalgebras and a unified braid condition – quantum Yang-Baxter equation. Finally, we will see how poetry has accompanied the development / history of the Yang–Baxter equation.

Although nanotechnology is a new and rapidly growing area of science, the impact of nanomaterials on living organisms is unknown in many aspects. In this regard it is extremely important to perform toxicological tests, but complete characterization of all varying preparations is extremely laborious. The computational technique called quantitative structure-activity relationship, or QSAR, allows reducing the cost of time- and resource-consuming nanotoxicity tests. In this review, (Q)SAR cytotoxicity studies of the past decade are systematically considered. We regard here five classes of engineered nanomaterials (ENMs): metal oxides, metal-containing nanoparticles, multi-walled carbon nanotubes, fullerenes, and silica nanoparticles. Some studies reveal that QSAR models are better than classification SAR models, while other reports conclude that SAR is more precise than QSAR. The quasi-QSAR method appears to be the most promising tool, as it allows accurately taking experimental conditions into account. However, experimental artifacts are a major concern in this case

The wave packet consisting of two harmonic plane waves with the same frequencies, but with different wave vectors is considered. The dispersion relation of a packet is structurally similar to the dispersion relation of a relativistic particle with a nonzero rest mass. The possibility of controlling the group velocity of a quasi-monochromatic wave packet by varying the angle between the wave vectors of its constituent waves is discussed.

We investigate entropy generation in unsteady three-dimensional axisymmetric MHD nanofluid flow over a non-linearly stretching sheet. The flow is subject to thermal radiation and a chemical reaction. The conservation equations were solved using the spectral quasi-linearization method. The novelty of the work is in the study of entropy generation in three-dimensional axisymmetric MHD nanofluid and the choice of the spectral quasilinearization method as the solution method. The effects of Brownian motion and thermophoresis are also taken into account when the nanofluid particle volume fraction on the boundary in passively controlled. The results show that as the Hartman number increases, both the Nusselt number and the Sherwood number decrease whereas the skin friction increases. It is further shown that an increase in the thermal radiation parameter corresponds to a decrease in the Nusselt number. Moreover, entropy generation increases with the physical parameters.

Zeta potential refers to the electrokinetic potential present in colloidal systems, exerting significant influence on the diverse properties of nano-drug delivery systems. The impact of dielectric constant on the zeta potential and charge inversion of highly charged colloidal particles, immersed in a variety of solvents spanning from polar, such as water, to nonpolar solvents, and in the presence of multivalent salts, was investigated through primitive model Monte Carlo (MC) simulations. Zeta potential ξ is decreased as decreasing the dielectric constant of the solvent, and upon further increase in the salinity and the valency of the salt. At elevated levels of salt, the colloidal particles become overcharged, in all solvents. As a result, their apparent charge becomes opposite in sign to the stoichiometric charge. This reversal of charge intensifies until reaching a saturation point with further increase in salinity, as a result, their apparent charge becomes opposite in sign to the stoichiometric charge. This reversal of charge intensifies until reaching a saturation point with further increase in salinity.

The direct measurements of the magnetocaloric effect were performed in Heusler Ni44.4Mn36.2Sn14.9Cu4.5 alloy at cryogenic temperatures in magnetic fields up to 10 T. The maximum value of the inverse magnetocaloric effect in a field of 10 T was observed ∆Tad = –2.7 K in the vicinity of the 1st order magnetostructural phase transition at T0 = 117 K. The ab initio and Monte Carlo calculations were performed to discuss the effect Cu doping into Ni-Mn-Sn compound on the ground state structural and magnetic properties. It is shown that with increasing Cu content the martensitic transition temperature decreases and the Curie temperature of austenite slightly increases. In general, calculated transition temperatures, magnetization values are correlated well with the experimental ones.

Effective household and individual disaster preparedness can minimize physical harm and property damage during catastrophic events. To assess the risk and vulnerability of affected areas, it is crucial for relief agencies to understand the level of public preparedness. Traditionally, government agencies have employed nationwide random telephone surveys to gauge the public’s attitudes and actions towards disaster preparedness. However, these surveys may lack generalizability in certain affected locations due to low response rates or areas not covered by the survey. To address this issue and enhance the comprehensiveness of disaster preparedness assessments, we develop a framework that seamlessly integrates machine learning and simulation. Our approach leverages machine learning algorithms to establish relationships between public attitudes towards disaster preparedness and demographic characteristics. Using Monte Carlo simulation, we generate datasets that incorporate demographic information of the affected location based on government-provided demographic distribution data. The generated dataset is then input into the machine learning model to predict the disaster preparedness attitudes of the affected population. We demonstrate the effectiveness of our framework by applying it to Miami-Dade County, where it accurately predicts the level of disaster preparedness. With this innovative approach, relief agencies can have a clearer and more comprehensive understanding of public disaster preparedness.

Water resource management and flood forecasting are crucial societal and financial stakes requiring reliable predictions of flow parameters (depth, velocity), the accuracy of which is often limited by uncertainties in hydrodynamic numerical models. In this study, we assess the effect of two uncertainty sources, namely breach characteristics induced by overtopping and the roughness coefficient, on water elevations and inundation extent. A two-dimensional (2D) hydraulic solver was applied in a Monte Carlo integration framework to a reach of the Loire river (France) including about 300 physical parameters. Inundation hazard maps for different flood scenarios allowed highlighting the impact of the breach development chronology. A special attention was paid to proposing a relevant sensitivity analysis to exam the factors influencing the depth and extent of flooding. The spatial analysis of vulnerability area induced by a levee breach width exhibits that, with increasing the flood discharge, the rise of the parameter influence is accompanied by a more localized spatial effect. This argues for a local analysis to allow a clear understanding of the flood hazard. The physical interpretation, highlighted by a global sensitivity analysis, showed precisely the dependence of the flood simulation on the main factors studied, i.e. the roughness coefficients and the characteristics of the breaches.

Dermal water content is an important biophysical parameter in preserving skin integrity and preventing skin damage. Traditional electrical-based and open-chamber evaporimeters have several well-known limitations. In particular, such devices are costly, sizeable, and only provide arbitrary outputs. They also do not permit continuous and non-invasive monitoring of dermal water content, which can be beneficial for various consumer, clinical and cosmetic purposes. We report here on the design and development of a digital multi-wavelength optical sensor that performs continuous and non-invasive measurement of dermal water content. In-silico investigation on porcine skin was carried out using the Monte Carlo modelling strategy to evaluate the feasibility and characterise the sensor. Subsequently, an in-vitro experiment was carried out to evaluate the performance of the sensor and benchmark its accuracy against a high-end, broad band spectrophotometer. Reference measurements were made against gravimetric analysis. The results demonstrate that the developed sensor can deliver accurate, continuous, and non-invasive measurement of skin hydration through measurement of dermal water content. Remarkably, the novel design of the sensor exceeded the performance of the high-end spectrophotometer due to the important denoising effects of temporal averaging. The authors believe, in addition to wellbeing and skin health monitoring, the designed sensor can particularly facilitate disease management in patients presenting diabetes mellitus, hypothyroidism, malnutrition, and atopic dermatitis.

The residential sector in the United States is in need of comprehensive policy-making reforms that concurrently address housing affordability and environmental sustainability. This study investigates the feasibility of state-wide zero-energy affordable housing by analyzing historical data on climate, energy use, and solar system costs in the Commonwealth of Virginia. The hypothesis examined is that the net present cost of implementation of rooftop residential solar systems to achieve zero-energy Low-Income Housing Tax Credit (LIHTC) buildings is lower than the discounted present cost of energy of otherwise identical conventional buildings that run without renewable energy generation systems. The authors propose a generalizable framework for analyzing the feasibility of achieving region- or state-wide zero-energy LIHTC developments. To validate the framework, the authors employ a longitudinal sample of monthly energy use data from 2013-2016 obtained from 310 residential units of 15 LIHTC developments across the state. Based on statistical regression analysis, energy simulation, and simulation-based risk analysis, the authors find that the net present value of zero-energy LIHTC investments can be positive with a low risk. The investment value varies often depending on the zero-energy building definition, weather characteristics, retail price of electricity, and incentive rate. This study can help housing policymakers and industry professionals analyze and benchmark the feasibility of innovative zero-energy housing policies and projects.

The reliability analysis allows to estimate the system's probability of detecting and identifying outlier. Failure to identify an outlier can jeopardise the reliability level of a system. Due to its importance, outliers must be appropriately treated to ensure the normal operation of a system. The system models are usually developed from certain constraints. Constraints play a central role in model precision and validity. In this work, we present a detailed optical investigation of the effects of the hard and soft constraints on the reliability of a measurement system model. Hard constraints represent a case in which there exist known functional relations between the unknown model parameters, whereas the soft constraints are employed for the case where such functional relations can slightly be violated depending on their uncertainty. The results highlighted that the success rate of identifying an outlier for the case of hard constraints is larger than soft constraints. This suggested that hard constraints should be used in the stage of pre-processing data for the purpose of identifying and removing possible outlying measurements. After identifying and removing possible outliers, one should set up the soft constraints to propagate the uncertainties of the constraints during the data processing. This recommendation is valid for outlier detection and identification purpose.

In this paper we evaluate the effects of hard and soft constraints on the Iterative Data Snooping (IDS), an iterative outlier elimination procedure. Here, the measurements of a levelling geodetic network were classified according to the local redundancy and maximum absolute correlation between the outlier test statistics, referred to as clusters. We highlight that the larger the relaxation of the constraints, the higher the sensitivity indicators MDB (Minimal Detectable Bias) and MIB (Minimal Identifiable Bias) for both the clustering of measurements and the clustering of constraints. There are circumstances that increase the family-wise error rate (FWE) of the test statistics, increase the performance of the IDS. Under a scenario of soft constraints, one should set out at least three soft constraints in order to identify an outlier in the constraints. In general, hard constraints should be used in the stage of pre-processing data for the purpose of identifying and removing possible outlying measurements. In that process, one should opt to set out the redundant hard constraints. After identifying and removing possible outliers, the soft constraints should be employed to propagate their uncertainties to the model parameters during the process of least-squares estimation.

Ability of quantitative structure – property / activity relationships (QSPRs/QSARs) to serve for epistemological processes in natural sciences is discussed. Some weirdness of QSPR/QSAR state-of-art are listed. There are some contradictions in the research results in this area. Sometimes, these should be classified as paradoxes or weirdness. These points often are ignored. Here these are listed and briefly commented. In addition, hypothesises on the future evolution of the QSPR/QSAR theory and practice are suggested.

Land Surface Temperature (LST) is a basic parameter in energy exchange between the land and atmosphere and is frequently used in many sciences such as climatology, hydrology, agriculture, ecology, etc. LST time series data have usually deficient, missing and unacceptable data caused by the presence of clouds in images, presence of dust in atmosphere and sensor failure. In this study, Singular Spectrum Analysis (SSA) algorithm was used to resolve the problem of missing and outlier data caused by cloud cover. The region studied in the present research included an image frame of MODIS with horizontal number 22 and vertical number 05 (h22v05). This image involved a large part of Iran and Turkmenistan and Caspian Sea. In this study, MODIS LST sensor (MOD11A1) was used during 2015 with 1×1 Km spatial resolution and day/night LST data (daily temporal resolution). The results of the data quality showed that cloud cover caused 36.37% of missing data in the studied time series with 730 day/night LST images. Further, the results of SSA algorithm in reconstruction of LST images indicated the Root Mean Square Error (RMSE) of 2.95 K between the original and reconstructed data in LST time series in the study region. In general, the findings showed that SSA algorithm using spatio-temporal interpolation in LST time series can be effectively used to resolve the problem of missing data caused by cloud cover.

In this paper, a class of coupled competitive neural networks, which subject to disturbance and discontinuous activation functions, is concerned. To realize the fixed-time quasi-bipartite synchronization, an aperiodic intermittent controller is initially designed. Subsequently, by combining the fixed-time stability theory and nonsmooth analysis, several criteria are established to ensure the bipartite synchronization in fixed time. Moreover, synchronization error bounds and settling time estimates are provided. Finally, numerical simulations are presented to verify the main results.

The purpose of the present paper is to introduce and investigate new subclasses of analytic function class of bi univalent function defined in open unit disk connected with a linear q-convolution operator, which are associated with the quasi-subordination. We find coefficients estimate |h_2 |,|h_3 | for functions in these subclasses. Several known and new consequences of these results are also pointed out.

Currently, efficient utilization of low-grade thermal energy is a great challenge. Thermoelectricity is an extremely promising method of generating electrical energy from temperature differences. As a promising energy conversion technology, thermo-electrochemical cells (TECs) have attracted much attention in recent years for their ability to convert thermal energy directly into electricity with high thermal power. Within TECs, anions and cations gain and lose electrons, respectively at the electrodes using the potential difference between the hot and cold terminals of the electrodes by redox couples. Additionally, the anions and cations therein are constantly circulating and mobile via concentration diffusion and thermal diffusion, providing an uninterrupted supply of power to the exterior. This review article focuses mainly on the operation of TECs, recent ad-vances in redox couples, electrolytes, and electrodes. The outlook for optimization strategies re-garding TECs is also presented in this paper.

This study investigates the application of Unmanned Aerial Vehicles (UAVs) equipped with a Micasense RedEdge-MX multispectral camera for the estimation of Secchi Depth (SD) in inland water bodies. The research analyzed and compared five sun-glint correction methodologies — Hedley, Goodman, Lyzenga, Joyce, and threshold removed glint to model the SD values derived from UAV multispectral imagery, highlighting the role of reflectance accuracy and algorithmic precision in SD modeling. While Goodman's method showed a higher correlation (0.92) with in situ SD measurements, Hedley's method exhibited the smallest average deviation (0.65 m), suggesting its potential in water resource management, environmental monitoring, and ecological modeling. The study also underscored the Quasi-Analytical Algorithm (QAA) potential in estimating SD due to its flexibility to process data from various sensors without requiring in situ measurements, offering scalability for large-scale water quality surveys. The accuracy of SD measures, calculated using QAA, was related to variability in water constituents of coloured dissolved organic matter and the solar zenith angle. A practical workflow for SD acquisition using UAVs and multispectral data was proposed for monitoring inland water bodies.

In response to the fire protection design challenges faced by the integrated station-city engineering under the new model of station-city fusion, such as large building volume, diverse traffic types, functional integration, spatial complexity, and inapplicable standards, this paper selects Chongqing East Station as a typical engineering case. The fire safety engineering methods and performance-based fire protection design ideas were used. The paper summarized the fire protection difficulties, such as multifunctional combination construction, large fire protection distance, large fire protection zoning area, long evacuation distance, and design of evacuation (quasi) safety zones. Combining with the actual situation of the engineering, this paper analyzed the basic technical measures taken to address these design challenges. By adopting the method of fire safety engineering, the safety of the design scheme is verified, which can provide a solution approach for the current fire protection design challenges in the integrated station-city engineering.

The paper presents a new conception of quasi-distributed sensor for simultaneous measurement of several physical fields and the results of an experimental study of this sensor. A distinctive feature of the sensor defined by the sensitive RC elements connected in the original tree-like structure. The proposed is the structure of the sensor and measurement system as well as corresponding measurement algorithm. The high accuracy demonstrated by the sensor’s prototype gives possibility to effective use the proposed sensors in many technical and scientific applications.

To effectively estimate airport terminal area capacity and assess the maximum throughput the sector can achieve when the capacity is given, this research proposes an approach to assess the terminal area capacity from the viewpoint of service provision resources. Terminal area capacity is optimized based on the equilibrium of air traffic service resource supply and demand. The supply-demand nexus is examined in consideration of terminal area route structure, traffic flow characteristics, and safety regulations. A flight service probability matrix and a terminal area demand and supply service time model are constructed to quantify resource expenditure at varied capacity levels. A optimization model is then developed to apportion maximal capacity under resource limitations. Model computation and computation results demonstrate the deviation between estimated and amended capacities is under 0.3 flight sorties per hour. The outcomes are congruent with historical statistics, thereby validating the accuracy and reliability of the model proposed in this study. Given capacity parameters, the model can deduce the maximal aircraft quantity served concurrently in terminal areas during peak periods. These revelations indicate the submitted model furnishes theoretical foundation and reference for terminal area sector partition and traffic alerting.

Raman-based Distributed Temperature Sensing (DTS) is a valuable tool for field testing and validating heat transfer models in Borehole Heat Exchanger (BHE) and Ground Source Heat Pump (GSHP) applications. However, temperature uncertainty is rarely reported in the literature. In this paper, a new calibration method is proposed for single-ended DTS configurations, along with a method to remove fictitious temperature drifts due to ambient air variations. The methods are implemented for a Distributed Thermal Response Test case study in an 800 m deep coaxial BHE. The results show that the calibration method and temperature drift correction are robust and give adequate results, with a temperature uncertainty of ± 0.58, 0.74, 1.1 and 1.7 K found at depths of 200, 300, 500 and 800 m, respectively. The temperature uncertainty increases non-linearly with depth and is dominated by the uncertainty in the calibrated parameters for depths larger than 200 m. The paper also offers insights into thermal features observed during the DTRT, including a heat flux inversion along the borehole depth and the slow temperature homogenization under circulation. Distributed temperature measurements and their uncertainties are useful tools in determining the strengths and drawbacks of ground heat transfer models.

3D facial surface imaging is a useful tool in dentistry and in terms of diagnostics and treatment planning. Between-groups PCA (bgPCA) is a method that has been used to analyse shapes in biological morphometrics, although various “pathologies” of bgPCA have recently been proposed. Monte Carlo (MC) simulated datasets were created here in order to explore “pathologies” of multilevel PCA (mPCA), where mPCA with two levels is equivalent to bgPCA. The first set of MC experiments involved 300 uncorrelated normally distributed variables, whereas the second set of MC experiments used correlated multivariate MC data describing 3D facial shape. We confirmed previous results of other researchers that indicated that bgPCA (and so also mPCA) can give a false impression of strong differences in component scores between groups when there is none in reality. These spurious differences in component scores via mPCA reduced strongly as the sample sizes per group were increased. Eigenvalues via mPCA were also found to be strongly effected by imbalances in sample sizes per group, although this problem was removed by using weighted forms of covariance matrices suggested by the maximum likelihood solution of the two-level model. However, this did not solve problems of spurious differences between groups in these simulations, which was driven by very small sample sizes in one group here. As a “rule of thumb” only, all of our experiments indicate that reasonable results are obtained when sample sizes per group in all groups are at least equal to the number of variables. Interestingly, the sum of all eigenvalues over both levels via mPCA scaled approximately linearly with the inverse of the sample size per group in all experiments. Finally, between-group variation was added explicitly to the MC data generation model in two experiments considered here. Results for the sum of all eigenvalues via mPCA predicted the asymptotic amount for the total amount of variance correctly in this case, whereas standard “single-level” PCA underestimated this quantity.

This study evaluates eight-year ownership costs for battery electric vehicles (BEV) versus non-plugin hybrid vehicles using forecasting to estimate future electricity and conventional gasoline prices and incorporating these in a multiple design of experiments simulation. Results suggest that while electric vehicles are statistically dominant in terms of variable costs over an 8-year life-span, high-performance hybrid non-plugins achieve variable fuel costs nearly as good as low-performing electric vehicles (those attaining only 3 miles per kilowatt hour) and that these hybrid acquisition costs are (on average) lower yet the vehicles retain higher residual values. In general, the six smallest ownership costs are split evenly between hybrid and electric vehicles; however, inflation for conventional regular gasoline is estimated to outstrip inflation per kilowatt hour. Thus, non-plugin hybrid cars are likely to require considerably more advanced engineering to keep pace.

Quantum theory is the land of indiscernible things, of things that in some situations cannot be put apart. This contrasts in much with standard mathematics and classical logic, which were elaborated to deal with \textit{individuals}, discernible things. Many thinkers, among them the mathematician Yuri Manin, have proposed that a more general theory of `sets' than standard set theories (where sets are collections of discernible elements) should be elaborated to cope with quantum physics. Quasi-set theory was proposed with such an aim. In this paper, we consider some aspects of this theory and discuss the way we can use such a theory as a framework for constructing deformable (not rigid) structures where indiscernible things can be considered but with one additional and fundamental detail: these structures cannot be extended to rigid structures as it happens in the standard frameworks. Really, the way of dealing with indiscernible (or indistinguishable) elements within a `standard' framework such as a `classical' set theory (like the ZFC system) requires the confinement of the discussion to deformable (non-rigid) structures. But in ZFC every structure can be extended to a rigid structure (encompassing only the trivial automorphism -- the identity function), so that in the extended structure (or in the whole universe) one can realize that the supposed indiscernible objects are not indiscernible at all. Quasi-set theory is such that the universe of quasi-sets is deformable and in such a theory one can construct deformable structures which cannot be extended to rigid ones. The aim of this paper is to put these things in a clear way.

Variations in the implant thread shape and occlusal load behavior may result in significant changes in the biological and mechanical properties of dental implants and surrounding bone tissue. Most previous studies consider a single implant thread design, an isotropic bone structure, and a static occlusal load. However, the effects of different thread designs, bone material properties, and loading conditions are important concerns in clinical practice. Accordingly, the present study performs Finite Element Analysis (FEA) simulations to investigate the static, quasi-static and dynamic response of the implant and implanted bone material under various thread designs and occlusal loading directions (buccal-lingual, mesiodistal and apical). The simulations focus specifically on the von Mises stress, displacement, shear stress, compressive stress and tensile stress within the implant and the surrounding bone. The results show that the thread design and occlusal loading rate have a significant effect on the stress distribution and deformation of the implant and bone structure during clinical applications. Overall, the results provide a useful insight into the design of enhanced dental implants for an improved load transfer efficiency and success rate.

In the wake of the current SARS-CoV-2 pandemic devastating the world, it is imperative to elucidate the comparative genomics of geographically-diverse strains of this novel coronavirus to gain insights into its microevolution, pathogenesis and control. Here we explore the molecular nature, genome-wide frequency, and gene-wise distribution of mutations in three distinct datasets encompassing 68 SARS-CoV-2 RNA-genomes altogether. While phylogenomic analysis revealed parallelism between the evolutionary paths charted by distinct quasispecies clusters of the virus, occurrence of mutations across genomes was found to be non-random. Whereas deletion mutations are extremely scarce and insertions totally absent, of all the instances of single nucleotide substitution detected, the overwhelming majority were transition mutations with cytidine to uridine being the most prevalent type. Propensity of this transition could be attributed to hydrolytic deamination mediated by ultra-violet irradiation or bisulfite reagent, both of which find wide usage as sterilizer/disinfectant. Transversions, albeit few and predominated by the guanosine to uridine form, were found concentrated in loci encoding the structural proteins of the virus, so might confer versatile tissue-colonization potentials. Mutation frequency of the three distinct genome-sets ranged narrowly between 0.07-1.08 × 10^-4 nucleotides positions mutated per nucleotide aligned. Gene-wise mapping of the global mutations illuminated the highly conserved nature of the genes encoding the non-structural proteins Nsp7, Nsp8 (two essential cofactors of the viral RNA-dependent RNA-polymerase) and Nsp9 (Nsp8-interacting single-strand RNA-binding protein), plus the envelope protein E (involved in SARS-CoV-2 assembly, budding and pathogenesis). These mutation-free genomic loci and/or their protein products could be potent targets for future drug designing/targeting.

The asymptotic dynamics of the classes of rational difference equations (RDEs) of third order defined over the positive real-line as $$\displaystyle{x_{n+1}=\frac{x_{n}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-1}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-2}}{ax_n+bx_{n-1}+cx_{n-2}}}$$ and $$\displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-1}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-2}}}$$ is investigated computationally with theoretical discussions and examples. It is noted that all the parameters $a, b, c$ and the initial values $x_{-2}, x_{-1}$ and $x_0$ are all positive real numbers such that the denominator is always positive. Several periodic solutions with high periods of the RDEs as well as their inter-intra dynamical behaviours are studied.

Abstract: Development of the fast neutron radiography with higher penetration capability in Non-Destructive Testing (NDT) techniques is rapidly. The detection efficiency and the imaging spatial resolution are the main problem of the radiography detector for many fast neutron imaging systems. In this paper, a detector composed of pixelated EJ200 scintillator array，an 16 × 16 silicon-photomultipliers (SiPM) array, and a capacitive multiplexing network readout electronics is designed. Using Monte-Carlo simulations，the main parameters of the detector was optimization. Then the detector prototype was fabricated and tested at a 14 MeV D-T neutron source. The preliminary test results showed that the prototype has spatial resolution of 1.26 mm. Using pixelated scintillator structure，the conflict between spatial resolution and detect efficiency could improved. Moreover, SiPMs offer a flexibility for the expansion of application field of the imaging system due to their excellent photon detection performance, relatively low price and joint possibility for large area.

To solve the problem of a uneven load for a dual-input counter-rotating transmission system, this research is based on Small Displacement Torsor Jacobian-Torsor theory.Considering the influence of the local parallel chain of a gearbox on the calculation accuracy of the backlash,the processing method of a meshing tooth pair parallel chain is proposed. The probability distribution of the tolerance band of gearbox manufacturing and assembly errors is modeled, the three-dimensional tolerance analysis model of a dual-input counter-rotating gearbox is established, and the Monte Carlo method is used to solve the backlash of the output end. The output backlash is solved with the Monte Carlo method, and the theoretical design backlash is combined with the output backlash to derive the gearbox side clearance range. The gearbox shaft system is analyzed statically, the range of the equal load coefficient is calculated, and the equal load performance optimization method based on the matching relationship between the gear backlash and the elastic shaft torsional stiffness is proposed. The correctness of the three-dimensional tolerance analysis combined with the elastic deflection angle used to optimize the equal load performance of the dual-input counter-rotating transmission system is verified by comparison with the relevant test data.

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.