Split supersymmetry under GUT and current dark matter constraints

There are four types of two-Higgs doublet models under a discrete \(Z_2\) symmetry imposed to avoid tree-level flavor-changing neutral current, i.e. type-I, type-II, type-X, and type-Y models. We investigate the possibility to discriminate the four models in the light of the flavor physics data, including \(B_s\) – \(\bar{B}_s\) mixing, \(B_{s,d} \rightarrow \mu ^+ \mu ^-\) , \(B\rightarrow \tau \nu \) and \(\bar{B} \rightarrow X_s \gamma \) decays, the recent LHC Higgs data, the direct search for charged Higgs at LEP, and the constraints from perturbative unitarity and vacuum stability. After deriving the combined constraints on the Yukawa interaction parameters, we have shown that the correlation between the mass eigenstate rate asymmetry \(A_{\Delta \Gamma }\) of \(B_{s} \rightarrow \mu ^+ \mu ^-\) and the ratio \(R=\mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{exp}/ \mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{SM}\) could be a sensitive probe to discriminate the four models with future precise measurements of the observables in the \(B_{s} \rightarrow \mu ^+ \mu ^-\) decay at LHCb.     

Optimal N-Point Configurations on the Sphere: “Magic” Numbers and Smale’s 7th Problem

This paper inquires into the concavity of the map \(N\mapsto v_s(N)\) from the integers \(N\ge 2\) into the minimal average standardized Riesz pair-energies \(v_s(N)\) of \(N\) -point configurations on the sphere \(\mathbb {S}^2\) for various \(s\in \mathbb {R}\) . The standardized Riesz pair-energy of a pair of points on \(\mathbb {S}^2\) a chordal distance \(r\) apart is \(V_s(r)= s^{-1}\left( r^{-s}-1 \right) \) , \(s \ne 0\) , which becomes \(V_0(r) = \ln \frac{1}{r}\) in the limit \(s\rightarrow 0\) . Averaging it over the \(\left( \begin{array}{c} N\\ 2\end{array}\right) \) distinct pairs in a configuration and minimizing over all possible \(N\) -point configurations defines \(v_s(N)\) . It is known that \(N\mapsto v_s(N)\) is strictly increasing for each \(s\in \mathbb {R}\) , and for \(s<2\) also bounded above, thus “overall concave.” It is (easily) proved that \(N\mapsto v_{-2}^{}(N)\) is even locally strictly concave, and that so is the map \(2n\mapsto v_s(2n)\) for \(s<-2\) . By analyzing computer-experimental data of putatively minimal average Riesz pair-energies \(v_s^x(N)\) for \(s\in \{-1,0,1,2,3\}\) and \(N\in \{2,\ldots ,200\}\) , it is found that the map \(N\mapsto {v}_{-1}^x(N)\) is locally strictly concave, while \(N\mapsto {v}_s^x(N)\) is not always locally strictly concave for \(s\in \{0,1,2,3\}\) : concavity defects occur whenever \(N\in {\mathcal {C}}^{x}_+(s)\) (an \(s\) -specific empirical set of integers). It is found that the empirical map \(s\mapsto {\mathcal {C}}^{x}_+(s),\ s\in \{-2,-1,0,1,2,3\}\) , is set-theoretically increasing; moreover, the percentage of odd numbers in \({\mathcal {C}}^{x}_+(s),\ s\in \{0,1,2,3\}\) is found to increase with \(s\) . The integers in \({\mathcal {C}}^{x}_+(0)\) are few and far between, forming a curious sequence of numbers, reminiscent of the “magic numbers” in nuclear physics. It is conjectured that these new “magic numbers” are associated with optimally symmetric optimal-log-energy \(N\) -point configurations on \(\mathbb {S}^2\) . A list of interesting open problems is extracted from the empirical findings, and some rigorous first steps toward their solutions are presented. It is emphasized how concavity can assist in the solution to Smale’s \(7\) th Problem, which asks for an efficient algorithm to find near-optimal \(N\) -point configurations on \(\mathbb {S}^2\) and higher-dimensional spheres.     

Some implications of inclusive electro- and neutrino production data

We systematically exploit the reported data on \(F_2^{\gamma p} ,F_2^{\gamma n} ,\sigma ^{vN} ,\sigma ^{\bar vN} ,\left\langle {xy} \right\rangle _{vN} ,\left\langle {xy} \right\rangle _{\bar vN} ,\left\langle {1 - y} \right\rangle _{vN} \) and \(\left\langle {1 - y} \right\rangle _{\bar vN} \) in order to test various versions of the quark parton model and to obtain further predictions.     

Confinement, average forces, and the Ehrenfest theorem for a one-dimensional particle

The topics of confinement, average forces, and the Ehrenfest theorem are examined for a particle in one spatial dimension. Two specific cases are considered: (i) A free particle moving on the entire real line, which is then permanently confined to a line segment or ‘a box’ (this situation is achieved by taking the limit V ₀?→?∞ in a finite well potential). This case is called ‘a particle-in-an-infinite-square-well-potential’. (ii) A free particle that has always been moving inside a box (in this case, an external potential is not necessary to confine the particle, only boundary conditions). This case is called ‘a particle-in-a-box’. After developing some basic results for the problem of a particle in a finite square well potential, the limiting procedure that allows us to obtain the average force of the infinite square well potential from the finite well potential problem is re-examined in detail. A general expression is derived for the mean value of the external classical force operator for a particle-in-an-infinite-square-well-potential, $\hat{F}$ . After calculating similar general expressions for the mean value of the position ( $\hat{X}$ ) and momentum ( $\hat{P}$ ) operators, the Ehrenfest theorem for a particle-in-an-infinite-square-well-potential (i.e., $\mathrm{d}\langle\hat{X}\rangle/\mathrm{d}t=\langle\hat{P}\rangle/M$ and $\mathrm{d}\langle\hat{P}\rangle/\mathrm{d}t=\langle\hat{F}\rangle$ ) is proven. The formal time derivatives of the mean value of the position ( $\hat{x}$ ) and momentum ( $\hat{p}$ ) operators for a particle-in-a-box are re-introduced. It is verified that these derivatives present terms that are evaluated at the ends of the box. Specifically, for the wave functions satisfying the Dirichlet boundary condition, the results, $\mathrm{d}\langle\hat{x}\rangle/\mathrm{d}t=\langle\hat{p}\rangle/M$ and $\mathrm{d}\langle\hat{p}\rangle/\mathrm{d}t=\mathrm{b.t.}+\langle\hat{f}\rangle$ , are obtained where b.t. denotes a boundary term and $\hat{f}$ is the external classical force operator for the particle-in-a-box. Thus, it appears that the expected Ehrenfest theorem is not entirely verified. However, by considering a normalized complex general state that is a combination of energy eigenstates to the Hamiltonian describing a particle-in-a-box with v(x)?=?0 ( $\Rightarrow\hat{f}=0$ ), the result that the b.t. is equal to the mean value of the external classical force operator for the particle-in-an-infinite-square-well-potential is obtained, i.e., $\mathrm{d}\langle\hat{p}\rangle/\mathrm{d}t$ is equal to $\langle\hat{F}\rangle$ . Moreover, the b.t. is written as the mean value of a quantity that is called boundary quantum force, f _B. Thus, the Ehrenfest theorem for a particle-in-a-box can be completed with the formula $\mathrm{d}\langle\hat{p}\rangle/\mathrm{d}t=\langle{{f_\mathrm{B}}}\rangle$ .     

Antineutrino science in KamLAND

The primary goal of KamLAND is a search for the oscillation of \({\bar{\nu }}_\mathrm{e}\) ’s emitted from distant power reactors. The long baseline, typically 180 km, enables KamLAND to address the oscillation solution of the “solar neutrino problem” with \({\bar{\nu }}_{e} \) ’s under laboratory conditions. KamLAND found fewer reactor \({\bar{\nu }}_{e} \) events than expected from standard assumptions about \(\overline{\nu }_e\) propagation at more than 9 \(\sigma \) confidence level (C.L.). The observed energy spectrum disagrees with the expected spectral shape at more than 5 \(\sigma \) C.L., and prefers the distortion from neutrino oscillation effects. A three-flavor oscillation analysis of the data from KamLAND and KamLAND + solar neutrino experiments with CPT invariance, yields \(\Delta m_{21}^2 \) = [ \(7.54_{-0.18}^{+0.19} \times \) 10 \(^{-5}\) eV \(^{2}\) , \(7.53_{-0.18}^{+0.19} \times \) 10 \(^{-5}\) eV \(^{2}\) ], tan \(^{2}\theta _{12}\) = [ \(0.481_{-0.080}^{+0.092} \) , \(0.437_{-0.026}^{+0.029} \) ], and sin \(^{2}\theta _{13}\) = [ \(0.010_{-0.034}^{+0.033} \) , \(0.023_{-0.015}^{+0.015} \) ]. All solutions to the solar neutrino problem except for the large mixing angle region are excluded. KamLAND also demonstrated almost two cycles of the periodic feature expected from neutrino oscillation effects. KamLAND performed the first experimental study of antineutrinos from the Earth’s interior so-called geoneutrinos (geo \({\bar{\nu }}_{e} \) ’s), and succeeded in detecting geo \({\bar{\nu }}_{e} \) ’s produced by the decays of \(^{238}\) U and \(^{232}\) Th within the Earth. Assuming a chondritic Th/U mass ratio, we obtain \(116_{-27}^{+28} {\bar{\nu }}_{e}\) events from \(^{238}\) U and \(^{232}\) Th, corresponding a geo \({\bar{\nu }}_{e}\) flux of \(3.4_{-0.8}^{+0.8}\times \) 10 \(^{6}\) cm \(^{-2}\)  s \(^{-1}\) at the KamLAND location. We evaluate various bulk silicate Earth composition models using the observed geo \({\bar{\nu }}_{e} \) rate.     

Pesin’s Entropy Formula for Systems Between C^1 and C^{1+\alpha }

In this article we give a new observation of Pesin’s entropy formula, motivated from Mañé’s proof of (Ergod Theory Dyn Sys 1:95–102, 1981). Let \(M\) be a compact Riemann manifold and \(f:\,M\rightarrow M\) be a \(C^1\) diffeomorphism on \(M\) . If \(\mu \) is an \(f\) -invariant probability measure which is absolutely continuous relative to Lebesgue measure and nonuniformly-H \(\ddot{\text {o}}\) lder-continuous(see Definition 1.1), then we have Pesin’s entropy formula, i.e., the metric entropy \(h_\mu (f)\) satisfies $$\begin{aligned} h_{\mu }(f)=\int \sum _{\lambda _i(x)> 0}\lambda _i(x)d\mu , \end{aligned}$$ where \(\lambda _1(x)\ge \lambda _2(x)\ge \cdots \ge \lambda _{dim\,M}(x)\) are the Lyapunov exponents at \(x\) with respect to \(\mu .\) Nonuniformly-H \(\ddot{\text {o}}\) lder-continuous is a new notion from probabilistic perspective weaker than \(C^{1+\alpha }.\)      

On strange quark scalar condensate in the nucleon

A well known difficulty with a large value of the σ term in πN scattering is analysed from positions of the QCD sum rules approach. The matrix element \(\left\langle {p\left| {\bar ss} \right|p} \right\rangle\) is related to flavour singlet correlation function of two quark condensates at zero momentum. The splittings \(\left\langle {p\left| {\bar uu - \bar ss} \right|p} \right\rangle\) and \(\left\langle {p\left| {\bar dd - \bar ss} \right|p} \right\rangle\) are calculated and turn to be in agreement withSU ₃ relations.     

Reunion Probabilities of N One-Dimensional Random Walkers with Mixed Boundary Conditions

In this work we extend the results of the reunion probability of \(N\) one-dimensional random walkers to include mixed boundary conditions between their trajectories. The level of the mixture is controlled by a parameter \(c\) , which can be varied from \(c=0\) (independent walkers) to \(c\rightarrow \infty \) (vicious walkers). The expressions are derived by using Quantum Mechanics formalism (QMf) which allows us to map this problem into a Lieb-Liniger gas (LLg) of \(N\) one-dimensional particles. We use Bethe ansatz and Gaudin’s conjecture to obtain the normalized wave-functions and use this information to construct the propagator. As it is well-known, depending on the boundary conditions imposed at the endpoints of a line segment, the statistics of the maximum heights of the reunited trajectories have some connections with different ensembles in Random Matrix Theory. Here we seek to extend those results and consider four models: absorbing, periodic, reflecting, and mixed. In all four cases, the probability that the maximum height is less or equal than \(L\) takes the form \(F_N(L)=A_N\sum _{\varvec{k}\in \Omega _{\text {B}}} \mathrm{e}^{-\sum _{j=1}^Nk_j^2}\mathcal {V}_N(\varvec{k})\) , where \(A_N\) is a normalization constant, \(\mathcal {V}_N(\varvec{k})\) contains a deformed and weighted Vandermonde determinant, and \(\Omega _{\text {B}}\) is the solution set of quasi-momenta \(\varvec{k}\) obeying the Bethe equations for that particular boundary condition.     

Is the CNMSSM more credible than the CMSSM?

“Post-sphaleron baryogenesis”, a fresh and profound mechanism of baryogenesis accounts for the matter–antimatter asymmetry of our present universe in a framework of Pati–Salam symmetry. We attempt here to embed this mechanism in a non-SUSY SO(10) grand unified theory by reviving a novel symmetry breaking chain with Pati–Salam symmetry as an intermediate symmetry breaking step and as well to address post-sphaleron baryogenesis and neutron–antineutron oscillation in a rational manner. The Pati–Salam symmetry based on the gauge group \(\mathrm{SU}(2)_L \times \mathrm{SU}(2)_{R} \times \mathrm{SU}(4)_C\) is realized in our model at \(10^{5}\) – \(10^{6}\)  GeV and the mixing time for the neutron–antineutron oscillation process having \(\Delta B=2\) is found to be \(\tau _{n-\bar{n}} \simeq 10^{8}\) – \(10^{10}\)  s with the model parameters, which is within the reach of forthcoming experiments. Other novel features of the model include low scale right-handed \(W^{\pm }_R\) , \(Z_R\) gauge bosons, explanation for neutrino oscillation data via the gauged inverse (or extended) seesaw mechanism and most importantly TeV scale color sextet scalar particles responsible for an observable \(n\) – \(\bar{n}\) oscillation which may be accessible to LHC. We also look after gauge coupling unification and an estimation of the proton lifetime with and without the addition of color sextet scalars.     

Production of π0 mesons and charged hadrons in\bar v neon and ν neon charged current interactions

The average multiplicities of charged hadrons and of π⁺, π^? and π⁰ mesons, produced in \(\bar v\) Ne and νNe charged current interactions in the forward and backward hemispheres of theW ^±-nucleon center of mass system, are studied with data from BEBC. The dependence of the multiplicities on the hadronic mass (W) and on the laboratory rapidity (y _Lab) and the energy fraction (z) of the pion is also investigated. Special care is taken to determine the π⁰ multiplicity accurately. The ratio of average π multiplicities \(\frac{{2\left\langle {n_{\pi ^O } } \right\rangle }}{{[\left\langle {n_{\pi ^ + } } \right\rangle + \left\langle {n_{\pi ^ - } } \right\rangle ]}}\) is consistent with 1. In the backward hemisphere \(\left\langle {n_{\pi ^O } } \right\rangle \) is positively correlated with the charged multiplicity. This correlation, as well as differences in multiplicities between \(\mathop v\limits^{( - )} \) and \(\mathop v\limits^{( - )} \) , \(\mathop v\limits^{( - )} \) scattering, is attributed to reinteractions inside the neon nucleus of the hadrons produced in the initial \(\mathop v\limits^{( - )} \) interaction.     

Comprehensive Bayesian analysis of rare (semi)leptonic and radiative \varvec{B} decays

The available data on \(|\Delta B| = |\Delta S| = 1\) decays are in good agreement with the Standard Model when permitting subleading power corrections of about \(15\,\%\) at large hadronic recoil. Constraining new-physics effects in \(\mathcal {C}_{7}^{\mathrm {}}\) , \(\mathcal {C}_{9}^{\mathrm {}}\) , \(\mathcal {C}_{10}^{\mathrm {}}\) , the data still demand the same size of power corrections as in the Standard Model. In the presence of chirality-flipped operators, all but one of the power corrections reduce substantially. The Bayes factors are in favor of the Standard Model. Using new lattice inputs for \(B\rightarrow K^*\) form factors and under our minimal prior assumption for the power corrections, the favor shifts toward models with chirality-flipped operators. We use the data to further constrain the hadronic form factors in \(B\rightarrow K\) and \(B\rightarrow K^*\) transitions.     

Maximal Distance Travelled by N Vicious Walkers Till Their Survival

We consider N Brownian particles moving on a line starting from initial positions \(\mathbf{{u}}\equiv \{u_1,u_2,\ldots u_N\}\) such that \(0 . Their motion gets stopped at time \(t_s\) when either two of them collide or when the particle closest to the origin hits the origin for the first time. For \(N=2\) , we study the probability distribution function \(p_1(m|\mathbf{{u}})\) and \(p_2(m|\mathbf{{u}})\) of the maximal distance travelled by the \(1^{\text {st}}\) and \(2^{\text {nd}}\) walker till \(t_s\) . For general N particles with identical diffusion constants \(D\) , we show that the probability distribution \(p_N(m|\mathbf{u})\) of the global maximum \(m_N\) , has a power law tail \(p_i(m|\mathbf{{u}}) \sim {N^2B_N\mathcal {F}_{N}(\mathbf{u})}/{m^{\nu _N}}\) with exponent \(\nu _N =N^2+1\) . We obtain explicit expressions of the function \(\mathcal {F}_{N}(\mathbf{u})\) and of the N dependent amplitude \(B_N\) which we also analyze for large N using techniques from random matrix theory. We verify our analytical results through direct numerical simulations.     

Line intensity measurements of methane’s \nu _3-band using a cw-OPO

We report on absolute line strength measurements of P(1), R(0) and R(1) singlet lines in the \(3.3\,\upmu\hbox {m}\,\nu _{3}\) (C–H stretching) band of methane \(^{12}\hbox {CH}_4\) at reference temperature \(T=296\)  K. Line strength measurements are performed at low pressure \((P \le 1\hbox { Torr})\) using direct absorption spectroscopy technique based on a widely tunable continuous-wave singly resonant optical parametric oscillator. The \(1\sigma \) overall accuracy in line strength determinations ranges between 7 and 8 % mostly limited by pressure and frequency measurements. A comparison with previous reported values is made. Our results show good agreement with the HITRAN 2012 database.     

Modeling of HOT (111) HgCdTe MWIR detector for fast response operation

The paper reports on photoelectrical performance of the mid-wave infrared (MWIR) (111) HgCdTe high operating temperature detector for the fast response conditions. Detector structure was simulated with software APSYS by Crosslight Inc. The detailed analysis of the time response as a function of device architecture and applied voltage was performed pointing out optimal working conditions. The time response of the MWIR HgCdTe detector with 50 % cut-off wavelength of \(\lambda _{c} \approx 5.3\, \upmu \hbox {m}\) at \(T = 200\)  K was estimated at the level of \(\tau _{s} \approx \) 2,500 ps for \(V = 100\)  mV and series resistance \(R_{Series} = 510\,\Omega \) . The series resistance’s reduction enables to reach \(\tau _{s}\approx 60\!-\!500\)  ps.     

Phase Transitions in Layered Systems

We consider the Ising model on \(\mathbb Z\times \mathbb Z\) where on each horizontal line \(\{(x,i), x\in \mathbb Z\}\) , called “layer”, the interaction is given by a ferromagnetic Kac potential with coupling strength \(J_{ \gamma }(x,y)={ \gamma }J({ \gamma }(x-y))\) , where \(J(\cdot )\) is smooth and has compact support; we then add a nearest neighbor ferromagnetic vertical interaction of strength \({ \gamma }^{A}\) , where \(A\ge 2\) is fixed, and prove that for any \(\beta \) larger than the mean field critical value there is a phase transition for all \({ \gamma }\) small enough.     

Passage Time Statistics in Exponential Distributed Time-Delay Models: Noisy Asymptotic Dynamics

The stochastic dynamics toward the final attractor in exponential distributed time-delay non-linear models is presented, then the passage time statistic is studied analytically in the small noise approximation. The problem is worked out by going to the associated two-dimensional system. The mean first passage time \(\left\langle t_{e}\right\rangle \) from the unstable state for this non-Markovian type of system has been worked out using two different approaches: firstly, by a rigorous adiabatic Markovian approximation (in the small mean delay-time \(\epsilon =\lambda ^{-1}\) ); secondly, by introducing the stochastic path perturbation approach to get a non-adiabatic theory for any \(\lambda \) . This first passage time distribution can be written in terms of the important parameters of the models. We have compared both approaches and we have found excellent agreement between them in the adiabatic limit. In addition, using our non-adiabatic approach we predict a crossover and a novel behavior for the relaxation scaling-time as a function of the delay parameter which for \(\lambda \ll 1\) goes as \(\left\langle t_{e}\right\rangle \sim 1/\sqrt{\lambda }\) .     

QCD analysis of forward neutron production in DIS

Observing light-by-light scattering at the large hadron collider (LHC) has received quite some attention and it is believed to be a clean and sensitive channel to possible new physics. In this paper, we study the diphoton production at the LHC via the process \({{pp}}\rightarrow {{p}}\gamma \gamma {{p}}\rightarrow {{p}}\gamma \gamma {{p}}\) through graviton exchange in the large extra dimension (LED) model. Typically, when we do the background analysis, we also study the double Pomeron exchange of \(\gamma \gamma \) production. We compare its production in the quark–quark collision mode to the gluon–gluon collision mode and find that contributions from the gluon–gluon collision mode are comparable to the quark–quark one. Our result shows, for extra dimension \(\delta =4\) , with an integrated luminosity \(\mathcal{L} = 200\,\mathrm{fb}^{-1}\) at the 14 TeV LHC, that diphoton production through graviton exchange can probe the LED effects up to the scale \({M}_{S}=5.06 (4.51, 5.11)\,\mathrm{TeV}\) for the forward detector acceptance \(\xi _1 (\xi _2, \xi _3)\) , respectively, where \(0.0015<\xi _1<0.5\) , \(0.1<\xi _2<0.5\) , and \(0.0015<\xi _3<0.15\) .     

Effective potential for local composite operators and the problem of factorization of multi-field condensates

A new formulation of effective potential for local composite operators is given. The two-fermion condensate \(\langle \bar \Psi \Psi \rangle \) and four-fermion condensate \(\langle \bar \Psi \Psi \bar \Psi \Psi \rangle \) are calculated simultaneously in the Gross-Neveu model up to next-to-the-leading order in 1/N expansion. It is shown that factorization \(\langle \bar \Psi \Psi \bar \Psi \Psi \rangle = C_1 \langle \bar \Psi \Psi \rangle ^2 \) holds only in theN→∞ limit and the non-factorized part of \(\langle \bar \Psi \Psi \bar \Psi \Psi \rangle \) contributed by the order-1/N terms is comparable to \(C_1 \langle \bar \Psi \Psi \rangle ^2 \) when takingN=3.     

Dilaton gravity, (quasi-) black holes,and scalar charge

Electrically charged dust is considered in the framework of Einstein–Maxwell–dilaton gravity with a Lagrangian containing the interaction term \(P(\chi )F_{\mu \nu }F^{\mu \nu }\) , where \(P(\chi )\) is an arbitrary function of the dilaton scalar field \(\chi \) , which can be normal or phantom. Without assumption of spatial symmetry, we show that static configurations exist for arbitrary functions \(g_{00} = \exp (2\gamma (x^{i}))\) ( \(i=1,2,3\) ) and \(\chi =\chi (\gamma )\) . If \(\chi = \mathrm{const}\) , the classical Majumdar–Papapetrou (MP) system is restored. We discuss solutions that represent black holes (BHs) and quasi-black holes (QBHs), deduce some general results and confirm them by examples. In particular, we analyze configurations with spherical and cylindrical symmetries. It turns out that cylindrical BHs and QBHs cannot exist without negative energy density somewhere in space. However, in general, BHs and QBHs can be phantom-free, that is, can exist with everywhere nonnegative energy densities of matter, scalar and electromagnetic fields.     

Semiclassical strings in supergravity PFT

In this paper, we introduce the bulk viscosity in the formalism of modified gravity theory in which the gravitational action contains a general function \(f(R,T)\) , where \(R\) and \(T\) denote the curvature scalar and the trace of the energy–momentum tensor, respectively, within the framework of a flat Friedmann–Robertson–Walker model. As an equation of state for a prefect fluid, we take \(p=(\gamma -1)\rho \) , where \(0 \le \gamma \le 2\) and a viscous term as a bulk viscosity due to the isotropic model, of the form \(\zeta =\zeta _{0}+\zeta _{1}H\) , where \(\zeta _{0}\) and \(\zeta _{1}\) are constants, and \(H\) is the Hubble parameter. The exact non-singular solutions to the corresponding field equations are obtained with non-viscous and viscous fluids, respectively, by assuming a simplest particular model of the form of \(f(R,T) = R+2f(T)\) , where \(f(T)=\alpha T\) ( \(\alpha \) is a constant). A big-rip singularity is also observed for \(\gamma <0\) at a finite value of cosmic time under certain constraints. We study all possible scenarios with the possible positive and negative ranges of \(\alpha \) to analyze the expansion history of the universe. It is observed that the universe accelerates or exhibits a transition from a decelerated phase to an accelerated phase under certain constraints of \(\zeta _0\) and \(\zeta _1\) . We compare the viscous models with the non-viscous one through the graph plotted between the scale factor and cosmic time and find that the bulk viscosity plays a major role in the expansion of the universe. A similar graph is plotted for the deceleration parameter with non-viscous and viscous fluids and we find a transition from decelerated to accelerated phase with some form of bulk viscosity.