charmm force field develeopment based on Monte Carlo simulation

charmm-mcsa/README.md

@@ -0,0 +1,43 @@
+Here's a README file based on the provided code:
+```
+# Monte Carlo Reweight Dihedral Force Constant
+
+This Python script implements Monte Carlo reweighting of dihedral force constants, primarily for molecular dynamics simulations. The code has been updated for Python 3 compatibility and improved practices.
+
+## Usage
+
+1. **Run the Script**: Execute the script using Python 3 with the required command-line arguments.
+   
+   ```bash
+   python3 monte_carlo_reweight.py <run> <thread> <jobid> <f>
+   ```
+
+   - `<run>`: Integer representing the run number.
+   - `<thread>`: Number of threads to use for parallel execution.
+   - `<jobid>`: Integer representing the job ID.
+   - `<f>`: String representing the filename.
+
+2. **Input Files**: Ensure the necessary input files are available, including parameters (`parm.parm`), QM data (`<f>.qm`), and CHARMM scripts (`charmm.<jobid>.sh`).
+
+3. **Output**: The script generates parameter files with optimized parameters (`optimized.<f>.<run>.str`) and intermediate parameter files during each iteration of the Monte Carlo simulation.
+
+## Dependencies
+
+- Python 3
+- `mcsa_module` (Assumed to be updated for Python 3 compatibility)
+- CHARMM (for molecular dynamics simulations)
+
+## Details
+
+- The script performs Monte Carlo simulations to optimize dihedral force constants based on a target function, using NMR data.
+- It utilizes multi-threading to speed up the optimization process.
+- Intermediate parameter files (`tmp.<jobid>.<f>.<run>.<thread>.<istep>.str`) are generated for each iteration of the simulation.
+- The best parameter set (`optimized.<f>.<run>.str`) is saved with the lowest RMSE value.
+- Execution time is displayed at the end of the script.
+
+## License
+
+This project is licensed under the MIT License - see the [LICENSE](LICENSE) file for details.
+```
+
+This README provides information on how to use the script, its dependencies, details about the implementation, and licensing information. Feel free to adjust or expand it further based on your needs.

reweight-trj-helix/README.md

@@ -0,0 +1,32 @@
+Sure, here's a simple README file based on the provided code:
+
+```
+# Monte Carlo Reweight Dihedral Force Constant
+
+This Python script performs Monte Carlo reweighting of dihedral force constants using RMSD with NMR J-coupling as a target function. It is designed to optimize parameters in molecular dynamics simulations.
+
+## Usage
+
+1. **Run the Script**: Execute the script using Python (version 2 or 3) with appropriate command-line arguments.
+   
+   ```bash
+   python monte_carlo_reweight.py
+   ```
+
+2. **Input Files**: The script requires input files containing information about dihedral angles, parameters, NMR data, and phi/psi angles. These files should be provided interactively during script execution.
+
+3. **Parameters**: Set the number of Monte Carlo steps (`nstep`) and the temperature (`tempr0`) in the script file.
+
+4. **Output**: The script generates parameter files (`reweight.<run>.str`) and energy files (`reweight.<run>.ene`) in the specified directory.
+
+## Dependencies
+
+- Python 2 or 3
+- Standard Python libraries (math, random, string, sys, copy)
+
+## License
+
+This project is licensed under the MIT License - see the [LICENSE](LICENSE) file for details.
+```
+
+This README provides a brief overview of the script, instructions on usage, dependencies, and licensing information. Feel free to expand or customize it further based on your needs.