The paper introduces the new version of the Mathematica package SARAH, which now supports non-supersymmetric (non-SUSY) models in addition to supersymmetric (SUSY) models. Key features include: 1. **Model Definition**: Enhanced support for defining non-SUSY models, including the ability to define component fields and global symmetries. 2. **Renormalization Group Equations (RGEs)**: Calculation of two-loop RGEs for general quantum field theories, including gauge dependence and kinetic mixing effects. 3. **Output Files**: Generation of model files for various tools such as FeynArts, CalcHep/CompHep, WHIZARD, and UFO format. 4. **SPHeno Integration**: Generation of source code for SPHeno to calculate the mass spectrum, precision observables, and decay widths. 5. **Vevacious Support**: Generation of model files for Vevacious to check for the global minimum of the one-loop effective potential. 6. **Lie Group Handling**: Improved handling of Lie groups, including support for non-$SU(N)$ gauge groups and non-fundamental representations. 7. **Global Symmetries**: Support for $Z_N$ and $U(1)$ global symmetries, with checks for charge conservation and gauge invariance. 8. **Numerical Evaluation**: Generation of files for numerical evaluation of RGEs using Mathematica's NDSolve. The paper also discusses the validation of the new features against independent Python code PyR8TE and provides examples of model files and commands for using SARAH 4.The paper introduces the new version of the Mathematica package SARAH, which now supports non-supersymmetric (non-SUSY) models in addition to supersymmetric (SUSY) models. Key features include: 1. **Model Definition**: Enhanced support for defining non-SUSY models, including the ability to define component fields and global symmetries. 2. **Renormalization Group Equations (RGEs)**: Calculation of two-loop RGEs for general quantum field theories, including gauge dependence and kinetic mixing effects. 3. **Output Files**: Generation of model files for various tools such as FeynArts, CalcHep/CompHep, WHIZARD, and UFO format. 4. **SPHeno Integration**: Generation of source code for SPHeno to calculate the mass spectrum, precision observables, and decay widths. 5. **Vevacious Support**: Generation of model files for Vevacious to check for the global minimum of the one-loop effective potential. 6. **Lie Group Handling**: Improved handling of Lie groups, including support for non-$SU(N)$ gauge groups and non-fundamental representations. 7. **Global Symmetries**: Support for $Z_N$ and $U(1)$ global symmetries, with checks for charge conservation and gauge invariance. 8. **Numerical Evaluation**: Generation of files for numerical evaluation of RGEs using Mathematica's NDSolve. The paper also discusses the validation of the new features against independent Python code PyR8TE and provides examples of model files and commands for using SARAH 4.