{"users":[{"id":23787,"username":"Sefika","name":"Şefika","avatar_template":"/letter_avatar_proxy/v4/letter/s/35a633/{size}.png","trust_level":0},{"id":18864,"username":"Michal_Husak","name":"Michal Hušák","avatar_template":"/user_avatar/matsci.org/michal_husak/{size}/19368_2.png","trust_level":1},{"id":3299,"username":"julian","name":"Julian Gale","avatar_template":"/letter_avatar_proxy/v4/letter/j/74df32/{size}.png","trust_level":3},{"id":22114,"username":"Emre_Genc","name":"Emre Genç","avatar_template":"/user_avatar/matsci.org/emre_genc/{size}/24488_2.png","trust_level":1},{"id":16298,"username":"Alexa","name":"Alexa","avatar_template":"/letter_avatar_proxy/v4/letter/a/e36b37/{size}.png","trust_level":1},{"id":25280,"username":"user17","name":"Franziska","avatar_template":"/letter_avatar_proxy/v4/letter/u/a9a28c/{size}.png","trust_level":0},{"id":25067,"username":"IanShuttleworth","name":"Ian Shuttleworth","avatar_template":"/user_avatar/matsci.org/ianshuttleworth/{size}/26982_2.png","trust_level":0},{"id":24852,"username":"Theo","name":"","avatar_template":"/letter_avatar_proxy/v4/letter/t/a87d85/{size}.png","trust_level":0},{"id":24343,"username":"Jared_Paulovich","name":"Jared Paulovich","avatar_template":"/user_avatar/matsci.org/jared_paulovich/{size}/22151_2.png","trust_level":0},{"id":24486,"username":"Grant_Roberts","name":"Grant Roberts","avatar_template":"/letter_avatar_proxy/v4/letter/g/c67d28/{size}.png","trust_level":0},{"id":24709,"username":"malahmer","name":"","avatar_template":"/letter_avatar_proxy/v4/letter/m/c89c15/{size}.png","trust_level":1},{"id":24503,"username":"Yongliang","name":"Yongliang YONG","avatar_template":"/letter_avatar_proxy/v4/letter/y/5e9695/{size}.png","trust_level":0},{"id":14000,"username":"jack1","name":"jack","avatar_template":"/letter_avatar_proxy/v4/letter/j/2bfe46/{size}.png","trust_level":1},{"id":20378,"username":"Stalminer","name":"Кирилл Романов","avatar_template":"/user_avatar/matsci.org/stalminer/{size}/22868_2.png","trust_level":1},{"id":13354,"username":"aephillips","name":"Anthony Phillips","avatar_template":"/user_avatar/matsci.org/aephillips/{size}/13831_2.png","trust_level":1},{"id":18099,"username":"Durjoy_Sarkar_Dhrubo","name":"Durjoy Sarkar Dhrubo","avatar_template":"/user_avatar/matsci.org/durjoy_sarkar_dhrubo/{size}/19853_2.png","trust_level":2},{"id":16560,"username":"Abid","name":"ullah","avatar_template":"/user_avatar/matsci.org/abid/{size}/24169_2.png","trust_level":1},{"id":23264,"username":"jhonnmick","name":"Jhonn Mick","avatar_template":"/user_avatar/matsci.org/jhonnmick/{size}/25580_2.png","trust_level":0},{"id":19770,"username":"maaliso","name":"Maaliso","avatar_template":"/user_avatar/matsci.org/maaliso/{size}/22097_2.png","trust_level":1},{"id":18337,"username":"ywang","name":"","avatar_template":"/letter_avatar_proxy/v4/letter/y/bbce88/{size}.png","trust_level":1},{"id":10774,"username":"mostafa","name":"","avatar_template":"/letter_avatar_proxy/v4/letter/m/f0a364/{size}.png","trust_level":1},{"id":22093,"username":"aravind","name":null,"avatar_template":"/letter_avatar_proxy/v4/letter/a/e8c25b/{size}.png","trust_level":1}],"primary_groups":[],"flair_groups":[],"topic_list":{"can_create_topic":false,"more_topics_url":"/c/gulp/26?page=1","per_page":30,"top_tags":[{"id":7,"name":"howto","slug":"howto"},{"id":17,"name":"gulp","slug":"gulp"},{"id":30,"name":"error","slug":"error"},{"id":8,"name":"lammps","slug":"lammps"},{"id":65,"name":"optimization","slug":"optimization"},{"id":1,"name":"bug-report","slug":"bug-report"},{"id":75,"name":"forcefield","slug":"forcefield"},{"id":78,"name":"lattice","slug":"lattice"},{"id":40,"name":"stress","slug":"stress"}],"topics":[{"fancy_title":"A huge difference of MgB2 bulk parameters","id":65388,"title":"A huge difference of MgB2 bulk parameters","slug":"a-huge-difference-of-mgb2-bulk-parameters","posts_count":13,"reply_count":6,"highest_post_number":13,"image_url":null,"created_at":"2025-09-17T12:24:31.506Z","last_posted_at":"2026-08-12T16:55:44.180Z","bumped":true,"bumped_at":"2026-08-12T16:55:44.180Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello, \nNowadays, I am trying to compute the bulk properties of the MgB2. In literature, the bulk parameters a=b=3.06 A, c= 3.52 A. However, the resulting bulk parameters are a=b= 2.88 A (approximately 6  error, acceptab&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":7,"name":"howto","slug":"howto"},{"id":17,"name":"gulp","slug":"gulp"},{"id":78,"name":"lattice","slug":"lattice"}],"tags_descriptions":{},"views":175,"like_count":6,"has_summary":false,"last_poster_username":"Emre_Genc","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":23787,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":22114,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"CO2 angle distortions with pGFN-FF","id":67412,"title":"CO2 angle distortions with pGFN-FF","slug":"co2-angle-distortions-with-pgfn-ff","posts_count":3,"reply_count":1,"highest_post_number":3,"image_url":null,"created_at":"2026-08-10T04:10:59.595Z","last_posted_at":"2026-08-10T08:46:07.167Z","bumped":true,"bumped_at":"2026-08-10T08:46:07.167Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Prof. Gale, \nI am currently running MD simulations using the pGFN-FF model and have observed some unusual bending of CO2 molecules confined within a framework. Although not all molecules are affected, there are case&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":26,"like_count":1,"has_summary":false,"last_poster_username":"Alexa","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":16298,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"GCMC for guest adsorption in porous materials","id":67278,"title":"GCMC for guest adsorption in porous materials","slug":"gcmc-for-guest-adsorption-in-porous-materials","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2026-06-25T01:54:10.978Z","last_posted_at":"2026-07-22T23:47:43.227Z","bumped":true,"bumped_at":"2026-07-22T23:47:43.227Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Prof. Gale, I hope you are doing well. \nI was testing GCMC simulations in GULP to estimate guest uptake in a porous material. Out of curiosity, I also tried the pGFN-FF model and was able to set up an initial input.&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":38,"like_count":1,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":16298,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"built in potetntial libraries","id":67149,"title":"built in potetntial libraries","slug":"built-in-potetntial-libraries","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2026-05-07T16:44:34.865Z","last_posted_at":"2026-05-08T00:19:05.645Z","bumped":true,"bumped_at":"2026-05-08T00:19:05.645Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello, \nI’m new to GULP and by searching for literature for the buckingham parameters I stumbled upon the built in libraries. For example, i saw that there should be a oxide.lib or metal.lib but i couldn’t find them? \nIs&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":38,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":25280,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Examples of GCMC of gas adsorption onto a surface","id":67004,"title":"Examples of GCMC of gas adsorption onto a surface","slug":"examples-of-gcmc-of-gas-adsorption-onto-a-surface","posts_count":3,"reply_count":0,"highest_post_number":3,"image_url":null,"created_at":"2026-03-28T15:33:16.292Z","last_posted_at":"2026-03-30T19:46:58.130Z","bumped":true,"bumped_at":"2026-03-30T19:46:58.130Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear users \nAre there any example gin files for GCMC of gas adsorption? \nSpecifically I am interested in simulating calorimetry data, and calculating the isosteric heat of adsorption. \nI’m struggling to work out where to&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":17,"name":"gulp","slug":"gulp"}],"tags_descriptions":{},"views":60,"like_count":0,"has_summary":false,"last_poster_username":"IanShuttleworth","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":25067,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Custom Potentials?","id":66857,"title":"Custom Potentials?","slug":"custom-potentials","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2026-02-11T09:41:04.017Z","last_posted_at":"2026-02-11T10:29:51.742Z","bumped":true,"bumped_at":"2026-02-11T10:29:51.742Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hi, \nI’m a masters student at Oxford looking at disordered lattice dynamics in alpha-AgI and wondering if it is possible to define a custom potenital in GULP or if the library of potentials is the only possible option.","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":52,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":24852,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"How to select proper a-Si Stillinger Weber parameters","id":66321,"title":"How to select proper a-Si Stillinger Weber parameters","slug":"how-to-select-proper-a-si-stillinger-weber-parameters","posts_count":6,"reply_count":3,"highest_post_number":6,"image_url":null,"created_at":"2025-11-20T17:06:01.781Z","last_posted_at":"2026-01-23T02:09:30.551Z","bumped":true,"bumped_at":"2026-01-23T02:09:30.551Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello, I am new to GULP and molecular dynamics in general, but I have been working on creating a VDOS plot for amorphous silicon. \nFor context, I ran an MD simulation of 2744 a-Si atoms in LAMMPS, from which I exported t&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":141,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":24343,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":24486,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"fit core shell potential with variables charge and constraint","id":66772,"title":"fit core shell potential with variables charge and constraint","slug":"fit-core-shell-potential-with-variables-charge-and-constraint","posts_count":5,"reply_count":1,"highest_post_number":5,"image_url":null,"created_at":"2026-01-16T20:35:31.117Z","last_posted_at":"2026-01-20T09:12:59.039Z","bumped":true,"bumped_at":"2026-01-20T09:12:59.039Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear all \ni’m a new user of gulp and i would use gulp software to create new core-shell potential for ZnO. i woud include core and shell charges of cation in the fit process so i have used the species and variables secti&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":68,"like_count":0,"has_summary":false,"last_poster_username":"malahmer","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":24709,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"File structure of the *.frc file","id":66687,"title":"File structure of the *.frc file","slug":"file-structure-of-the-frc-file","posts_count":6,"reply_count":1,"highest_post_number":6,"image_url":null,"created_at":"2025-12-21T21:21:09.693Z","last_posted_at":"2025-12-29T00:37:20.696Z","bumped":true,"bumped_at":"2025-12-29T00:37:20.696Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello, \nI compute the force constants through the phonon calculations, but a very complex file is generated in which I can not extract the force constant of relevant bond. \nCould you please help me to understand the file&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":7,"name":"howto","slug":"howto"},{"id":17,"name":"gulp","slug":"gulp"}],"tags_descriptions":{},"views":92,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":22114,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"How to write the Tersoff_ZRL form in GULP","id":66587,"title":"How to write the Tersoff_ZRL form in GULP","slug":"how-to-write-the-tersoff-zrl-form-in-gulp","posts_count":4,"reply_count":0,"highest_post_number":4,"image_url":null,"created_at":"2025-12-12T00:21:00.702Z","last_posted_at":"2025-12-15T04:33:05.188Z","bumped":true,"bumped_at":"2025-12-15T04:33:05.188Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian and other users: \nI am a beginner in the GULP software. Currently, I am using GULP to study the BeO system, whose interatomic potentials are in the Tersoff_ZRL form. However, I am not sure how to convert the &hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":67,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":24503,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Dynamical matrices decimal point digits","id":47248,"title":"Dynamical matrices decimal point digits","slug":"dynamical-matrices-decimal-point-digits","posts_count":6,"reply_count":2,"highest_post_number":6,"image_url":null,"created_at":"2023-03-03T20:59:21.030Z","last_posted_at":"2025-11-04T01:15:04.849Z","bumped":true,"bumped_at":"2025-11-04T01:15:04.849Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hi, \nI am interested in doing some lattice dynamics calculations that require high precisions. \nWhen I looked at the dynamical matrix outputs from GULP, it seems to print out 6 decimal point digits. \nDoes this mean inter&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":318,"like_count":1,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":14000,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":20378,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Apparent bug handling alternative space group settings","id":64904,"title":"Apparent bug handling alternative space group settings","slug":"apparent-bug-handling-alternative-space-group-settings","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2025-08-07T10:48:03.706Z","last_posted_at":"2025-08-08T00:57:08.239Z","bumped":true,"bumped_at":"2025-08-08T00:57:08.239Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian and all, \nI believe I’ve encountered a small bug whereby non-standard settings of space group 38, Amm2, are interpreted as settings of space group 35, Cmm2. \nI think this can be fixed by a very small patch in&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":1,"name":"bug-report","slug":"bug-report"}],"tags_descriptions":{},"views":61,"like_count":1,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":13354,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Phonon Density of State (PDOS) using GULP","id":64733,"title":"Phonon Density of State (PDOS) using GULP","slug":"phonon-density-of-state-pdos-using-gulp","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":"https://matsci.org/uploads/default/original/3X/3/7/37db0aa2ec0d101b29d873502c0b5b8b0e16105f.jpeg","created_at":"2025-07-26T19:13:22.133Z","last_posted_at":"2025-08-01T00:10:49.001Z","bumped":true,"bumped_at":"2025-08-01T00:10:49.001Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Is it possible to calculate the in-plane and out-of-plane phonon DOS of a double-wall carbon nanotube (DWCNT) using GULP, as shown below? Please suggest some resources so that I can easily calculate it. I am new to GULP.&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":91,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":18099,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"How to plot phonon dispersion from GULP output files","id":64524,"title":"How to plot phonon dispersion from GULP output files","slug":"how-to-plot-phonon-dispersion-from-gulp-output-files","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2025-07-13T12:20:45.621Z","last_posted_at":"2025-07-14T00:44:29.090Z","bumped":true,"bumped_at":"2025-07-14T00:44:29.090Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear GULP users, \nI am new to using GULP and currently working on calculating the phonon dispersion of my system. I would also like to visualize the corresponding eigenvalues. However, I am not sure how to extract and pl&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":98,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":16560,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"How to enhance the accuracy of the Phonon calculations?","id":64476,"title":"How to enhance the accuracy of the Phonon calculations?","slug":"how-to-enhance-the-accuracy-of-the-phonon-calculations","posts_count":4,"reply_count":1,"highest_post_number":4,"image_url":null,"created_at":"2025-07-10T19:38:40.811Z","last_posted_at":"2025-07-12T17:34:18.524Z","bumped":true,"bumped_at":"2025-07-13T01:10:56.425Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello all, \nI am trying to understand the stability of the hypothetical MOF crystal. According to my analysis, there is only one imaginary frequency. In DFT calculations, this only one negative frequency can be mitigated&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":7,"name":"howto","slug":"howto"},{"id":17,"name":"gulp","slug":"gulp"}],"tags_descriptions":{},"views":186,"like_count":0,"has_summary":false,"last_poster_username":"jhonnmick","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":22114,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":23264,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Interaction energy with GFN-FF","id":64442,"title":"Interaction energy with GFN-FF","slug":"interaction-energy-with-gfn-ff","posts_count":3,"reply_count":1,"highest_post_number":3,"image_url":null,"created_at":"2025-07-09T08:16:26.503Z","last_posted_at":"2025-07-11T07:19:58.347Z","bumped":true,"bumped_at":"2025-07-11T07:19:58.347Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear prof. Gale, \nWhen calculating the interaction energy of an A–B complex using pGFN-FF (E_AB - E_A - E_B)​, I understand that computing the single-point energies of A and B separately causes the program to independent&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":86,"like_count":1,"has_summary":false,"last_poster_username":"Alexa","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":16298,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Pgfnf optimized molecule explode (sulphur polymorph)","id":63972,"title":"Pgfnf optimized molecule explode (sulphur polymorph)","slug":"pgfnf-optimized-molecule-explode-sulphur-polymorph","posts_count":9,"reply_count":7,"highest_post_number":9,"image_url":null,"created_at":"2025-06-11T14:11:54.486Z","last_posted_at":"2025-06-26T07:57:06.967Z","bumped":true,"bumped_at":"2025-06-26T19:07:34.222Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hi \nA sulphur modification explode after optimization in GULP by pgfnff. \nIt is not a charged system  = no ions, not the already discussed problem. \nWhen optimized by DFT a large cut off must be used to descipbe week int&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":144,"like_count":0,"has_summary":false,"last_poster_username":"Michal_Husak","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Calculating electrostatic potential in GULP to be able to add it to MACE energy","id":64169,"title":"Calculating electrostatic potential in GULP to be able to add it to MACE energy","slug":"calculating-electrostatic-potential-in-gulp-to-be-able-to-add-it-to-mace-energy","posts_count":5,"reply_count":3,"highest_post_number":5,"image_url":null,"created_at":"2025-06-21T11:05:53.122Z","last_posted_at":"2025-06-22T19:53:32.199Z","bumped":true,"bumped_at":"2025-06-22T19:53:32.199Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"I would like to use GULP for calculating electrostatic lattice energy as an external correction term to energy calculated by MACE in ASE environment. The work is targeted on molecular crystals - salts of organic acids. \n&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":126,"like_count":1,"has_summary":false,"last_poster_username":"Michal_Husak","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"GULP and MACE","id":64134,"title":"GULP and MACE","slug":"gulp-and-mace","posts_count":3,"reply_count":1,"highest_post_number":3,"image_url":null,"created_at":"2025-06-19T16:49:15.259Z","last_posted_at":"2025-06-20T05:50:12.151Z","bumped":true,"bumped_at":"2025-06-20T05:58:38.196Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Is there any chance to add to GULP support for MACE ML force field ? MACE is alredy integrated in ASE, Lamps, OpenMM. My test in ACE shows MACE outperforms any MM force field in quality of results. And atom typing is aut&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":125,"like_count":0,"has_summary":false,"last_poster_username":"Michal_Husak","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Charge Neutrality Problem","id":63928,"title":"Charge Neutrality Problem","slug":"charge-neutrality-problem","posts_count":8,"reply_count":3,"highest_post_number":8,"image_url":null,"created_at":"2025-06-09T13:13:23.986Z","last_posted_at":"2025-06-12T08:14:14.661Z","bumped":true,"bumped_at":"2025-06-12T08:14:14.661Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear all, \nI am a novice user of GULP and this is my first simulation to run with it. \nI am trying to use the core-shell model with the Buckingham potential implemented in GULP to optimize the unit cell of BaTiO3. I am f&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":17,"name":"gulp","slug":"gulp"},{"id":65,"name":"optimization","slug":"optimization"}],"tags_descriptions":{},"views":153,"like_count":3,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":19770,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"The result of GULP run depends on MPI parallelizm","id":63654,"title":"The result of GULP run depends on MPI parallelizm","slug":"the-result-of-gulp-run-depends-on-mpi-parallelizm","posts_count":5,"reply_count":2,"highest_post_number":5,"image_url":null,"created_at":"2025-05-23T11:40:36.713Z","last_posted_at":"2025-06-04T20:32:46.987Z","bumped":true,"bumped_at":"2025-06-04T20:32:46.987Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hi \nWe do test of GULP with gfnff … We use MPI parallelized GULP. All works OK e.g. for X23 test structures (e.g. CYHEXO) = different MPI parallelization gives identical number of cycles and identical results. \nFor a mor&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":110,"like_count":1,"has_summary":false,"last_poster_username":"Michal_Husak","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"WARNING : Breaking symmetry due to perfect biradical","id":63452,"title":"WARNING : Breaking symmetry due to perfect biradical","slug":"warning-breaking-symmetry-due-to-perfect-biradical","posts_count":5,"reply_count":2,"highest_post_number":5,"image_url":null,"created_at":"2025-05-12T14:55:59.980Z","last_posted_at":"2025-05-14T07:02:54.002Z","bumped":true,"bumped_at":"2025-05-14T07:02:54.002Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello Dear, \nI have received such an error in MOF crystal with P1 symmetry along with the charge calculation error. \nWhat might be the possible reason? \nThanks in advance. \nEmre \n\nINPUT: \nopti conv conse qok nomod pres s&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":30,"name":"error","slug":"error"}],"tags_descriptions":{},"views":90,"like_count":0,"has_summary":false,"last_poster_username":"Emre_Genc","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":22114,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Difference in Four-body potentials between GULP and LAMMPS","id":63319,"title":"Difference in Four-body potentials between GULP and LAMMPS","slug":"difference-in-four-body-potentials-between-gulp-and-lammps","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":"https://matsci.org/uploads/default/optimized/3X/8/2/82b2d9be0c757e55a69f745a8d45a63e0ff10f22_2_1024x171.png","created_at":"2025-05-03T16:58:02.192Z","last_posted_at":"2025-05-12T00:57:10.817Z","bumped":true,"bumped_at":"2025-05-12T00:57:10.817Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian and other users, \nI am using GULP to compute the single point energy and compare with lammps to verify that everything is set correctly in GULP. \nthe output shows everything is in good agreement with lammps (&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":49,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":18337,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Bond parameter fitting","id":62817,"title":"Bond parameter fitting","slug":"bond-parameter-fitting","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":null,"created_at":"2025-04-03T16:47:02.079Z","last_posted_at":"2025-04-04T02:15:50.116Z","bumped":true,"bumped_at":"2025-04-04T02:15:50.116Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian, \nI am confused by the fitting procedure of gulp to parameterize the bond terms (K and r0 in harmonic functional form). \nMy goal is to release the stress to 0 GPa and to keep the volume the same as the experi&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":55,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":18337,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Wrong bonding in output cif file","id":60132,"title":"Wrong bonding in output cif file","slug":"wrong-bonding-in-output-cif-file","posts_count":15,"reply_count":10,"highest_post_number":15,"image_url":"https://matsci.org/uploads/default/original/3X/7/7/7780d313a413ac44158cd9a9d1ff0dd837816ca3.png","created_at":"2025-01-09T18:47:13.977Z","last_posted_at":"2025-03-27T23:55:44.812Z","bumped":true,"bumped_at":"2025-03-28T14:27:24.115Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Users, \nMy crystal is a metal-organic framework, and I want to use GULP to perform a series of optimizations. Before that, I think it’s necessary to make sure the code describes the crystal correctly. \nthere are 6 a&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[{"id":7,"name":"howto","slug":"howto"}],"tags_descriptions":{},"views":279,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":18337,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":18864,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Writing the charges and bond order","id":62314,"title":"Writing the charges and bond order","slug":"writing-the-charges-and-bond-order","posts_count":5,"reply_count":1,"highest_post_number":5,"image_url":null,"created_at":"2025-03-25T14:12:31.009Z","last_posted_at":"2025-03-27T08:32:24.448Z","bumped":true,"bumped_at":"2025-03-27T08:32:24.448Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Hello, \nI am trying to output the charge and bond order during molecular dynamics simulation using ReaxFF.  The line in the input file related to this is: \n\noutput qbo append charge_BO \n\nWhen I run GULP (6.3.2 latest ver&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":84,"like_count":0,"has_summary":false,"last_poster_username":"mostafa","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":10774,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Deviation in Buckingham Potential Parameter Optimization","id":61968,"title":"Deviation in Buckingham Potential Parameter Optimization","slug":"deviation-in-buckingham-potential-parameter-optimization","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":"https://matsci.org/uploads/default/original/3X/0/6/06289b3b429e632df03eb1b2b370ae5f1835698f.png","created_at":"2025-03-17T16:12:08.334Z","last_posted_at":"2025-03-25T00:18:01.853Z","bumped":true,"bumped_at":"2025-03-25T00:18:01.853Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"I have obtained the Buckingham potential parameters (A, ρ, and C) for individual atomic interactions from the literature. I currently have only one observable (spring constant) to optimize these parameters. When I attemp&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":74,"like_count":0,"has_summary":false,"last_poster_username":"aravind","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest single","description":"Original Poster, Most Recent Poster","user_id":22093,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"How to obtain the internal stress of initial geometry","id":61894,"title":"How to obtain the internal stress of initial geometry","slug":"how-to-obtain-the-internal-stress-of-initial-geometry","posts_count":6,"reply_count":2,"highest_post_number":6,"image_url":null,"created_at":"2025-03-14T13:25:30.384Z","last_posted_at":"2025-03-21T00:20:10.666Z","bumped":true,"bumped_at":"2025-03-21T00:20:10.666Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian and all users, \nI am thinking about how to obtain the static stress/pressure of my system with a force field, and comparing it with the results from lammps. \nI saw keywords like ‘strain’ + ‘stress_out’ or ass&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":92,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":18337,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Difference in single point energy between GULP and LAMMPS","id":62002,"title":"Difference in single point energy between GULP and LAMMPS","slug":"difference-in-single-point-energy-between-gulp-and-lammps","posts_count":3,"reply_count":1,"highest_post_number":3,"image_url":"https://matsci.org/uploads/default/optimized/3X/e/7/e7e2a2790178b94cc3d1311a7b346974e18d0c4f_2_1024x223.png","created_at":"2025-03-18T18:06:45.334Z","last_posted_at":"2025-03-19T17:56:01.490Z","bumped":true,"bumped_at":"2025-03-19T17:56:01.490Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Dear Julian and users, \nI have computed the single point energy of my periodic MOF system in a cubic box with unit cell length of 26.3 Å. \nThe 1-4 electrostatic and van der Waals interactions in my force field are both s&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":85,"like_count":0,"has_summary":false,"last_poster_username":"ywang","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":"latest","description":"Original Poster, Most Recent Poster","user_id":18337,"primary_group_id":null,"flair_group_id":null},{"extras":null,"description":"Frequent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]},{"fancy_title":"Variation in Anion-Anion Core-Shell Charges Across Different Cases","id":61836,"title":"Variation in Anion-Anion Core-Shell Charges Across Different Cases","slug":"variation-in-anion-anion-core-shell-charges-across-different-cases","posts_count":2,"reply_count":0,"highest_post_number":2,"image_url":"https://matsci.org/uploads/default/original/3X/d/2/d2c7ff37374f070250387e4a2d668a4c20a26b78.png","created_at":"2025-03-12T16:25:47.452Z","last_posted_at":"2025-03-13T01:42:03.980Z","bumped":true,"bumped_at":"2025-03-13T01:42:03.980Z","archetype":"regular","unseen":false,"pinned":false,"unpinned":null,"excerpt":"Why do anion-anion core-shell charge distributions vary across different cases? \nWhat factors influence these charge variations in a core-shell model? \nAre there specific guidelines to ensure consistency in charge assign&hellip;","visible":true,"closed":false,"archived":false,"bookmarked":null,"liked":null,"tags":[],"tags_descriptions":{},"views":42,"like_count":0,"has_summary":false,"last_poster_username":"julian","category_id":26,"op_like_count":0,"pinned_globally":false,"featured_link":null,"has_accepted_answer":false,"can_vote":false,"posters":[{"extras":null,"description":"Original Poster","user_id":22093,"primary_group_id":null,"flair_group_id":null},{"extras":"latest","description":"Most Recent Poster","user_id":3299,"primary_group_id":null,"flair_group_id":null}]}]}}